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Laminar Compressible Boundary Laver 
Along a Curved Insulated Surface 


* 


I. FLUGGE-LOTZ?+ ann ARLO F. JOHNSON? 


Stanford University 


ABSTRACT 


Several approximate methods exist for computing the com 


pressible laminar boundary-layer flow along a curved surface 


with given pressure gradient and given wall temperature. The 
boundary-layer flow for given pressure gradient and insulated wall 
is an essentially different problem, because for Prandtl Number 
unequal to one the wall temperature is unknown and a rather 
accurate computation of this temperature is desirable. In this 
paper a method is developed for solving the laminar compressible 
boundary layer over insulated bodies of arbitrary shape via the 
integrated momentum and energy equations and the expression 
of velocity and temperature profiles by polynomials. Examples 
are worked out for a flat plate which is compared with Crocco’s 
solution, for an arbitrary pressure distribution which is com 
pared with the solution of Fliigge-Lotz and Eichelbrenner, for a 
biconvex airfoil, and for an inflected wall. The method gives 
results accurate to a few per cent so long as the flow is accelerating 
Like the Pohlhausen technique—its counterpart in incompres 
sible flow—it is less accurate in decelerating flows and predicts 
separation too far downstream. The time required to calculate 
the drag coefficient, boundary-layer thickness, velocity, and 
temperature profiles at any given section using this method is 
of the order of a few hours 

The solutions presented indicate that compressibility has an 
appreciable effect on velocity, temperature, drag, and boundary 
layer thickness. Prandtl Number variation appears to create 
a significant change only in the temperature profile 

A method that Gruschwitz developed to solve the same prob 
lem is briefly discussed at the end of this paper 


NOTATION 


= arc length along the wall, ft 
coordinate normal to the wall, ft 
velocity in the x direction, ft. per sec 
velocity in the y direction, ft. per sec. 


fluid density, slugs per ft 


Received April 7, 1954 

* This study is part of a research project at Stanford Uni 
versity sponsored by the Office of Naval Research [Contract 
N6-onr 251 Task Order 6 (NR-035-104) 

t Lecturer, Department of Mechanical Engineering, Stanford 
University 

t Formerly, Graduate Student at Stanford University. Now, 
Associate Professor, Department of Mechanical Engineering, 
University of Utah. 


fluid pressure, lb. per sq-ft 

absolute viscosity, slugs per ft. sc« 

u/p kinematic viscosity, sq.ft. per sec 

u(Ou/Oy) viscous shear stress, lb. per sq.ft 

specific heat at constant pressure, sq.ft. per 
sec.2 °R 

specific heat at constant volume, sq.ft. per 
sec 2 R 

Cp/Cr (1.4 for air) ratio of the specific heats 

absolute temperature, °R 

cyT (perfect gas) enthalpy, sq.ft. per sec 

ie | u*?/2) total specific energy, sq ft. per 
sec 2 

thermal conductivity, lb. per sec. °R 

uCp)/(Rr), Prandtl Number 

gas constant in the equation of state 

empirical constant in the Sutherland Equa 
tion, °R 

C/T), modified Sutherland constant 

speed of sound, ft per sec 

u/a, Mach Number 

(4 1)/2} M2 = u,?/2i, compressibility 
variable 

local drag coefficient 


total drag coefficient 


pu l : 
( l ) dy momentum thick 
Pelle Ue 


ness, ft 


a. 


in pu y 
(: ) dy displacement thick- 


Pelle 


0 


. 


0 
ness, ft 

U/ Ue, velocity ratio 

p-/p, density ratio 

7/7T,, temperature ratio 

normalized y coordinat¢ 

boundary-layer ‘‘thickness’’ used in de 
fining 7», ft 

coefficients in the velocity polynomial 

p-/Pw, basic parameter in the density poly 
nomial 

coefficients in the density polynomial |Eq 
15 

coefficients in the density polynomial [Eq 
(16) 


basic parameter in the velocity polynomial 
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X = (b'u,)/(bu’,), abbreviated notation for this This equates flow-wise (x) components of force to M: 
aenned momentum change. The equation for y components ) 
K = alternate parameter to \ pam vitae bythe 't : 
a generates to wis 
L = characteristic length, ft enemies oe ; 
Re = (uL)/vm, Reynolds Number, based on the 
; Op/oy = 0 ‘ | 
reference values - ) 
ti, fo, ete. = functions of d | f ti t thi f the | len 
cause oO > extreme 155 — 
B Re ee ee eECAUSE 1€ ex reme hinnes O he oundary layer ex] 
R = (b — 1)/m, abbreviated notation of the and the deductions that follow therefrom. Phi 
given ratio (3) The energy equation, expressed in alternate 
: > : ret 
AB = tabulated function of A forms: 
ky = constant associated with the isentropic en ne! 
ergy equation [ pu(O7, Ox) | + | pu(Oj oy)] = enc 
Subscripts and Superscripts (O oy) (ur) + (0 Oy) [k(OT oy (4a ave 
e = value in the flow just external to the the 
boundary layer = (1/Pr) (0/Oy) [u(07/dy)] - is t 
w = value at the wall (1 — Pr) Pr] (O/Ov) (ur (4b bou 
) = value at the reference point in the flow : 1 
(free-stream values have been taken here) It is derived under the usual assumption that there js 
sep = value at the separation point as eeuseniaee h : f os in li : Eqs 
. ; O % Cli > excNange OF energy - x “i 
stag = value at the stagnation point re Ce : —— Cirection. to 4 
, oo Sie Because yw and k, vary with temperature, while Pr js ” 2 
‘ . ; Sa : ; iné 
. = 02/dx? constant for a given fluid, Eq. (4b) proves preferable rye 
Viiuhtcs tn imanitaniees Bove for mathematical analysis. It reveals by inspection 
° ° ° /J0 
vy = x/L the particular solution, 7 = constant when Pr = |, 
/ 7 . ~ . °,° e 
ya = (¥V Re)/L which fits the physical boundary conditions of flow over 
2 me apes an insulated wall. 
q¢ = e/ lo _ . = ii 
tqui state a perfect gas: The 
i? ws a (4) Equation of state for a perfect gas: 
ta’ = (Lu,')/uo oe 
= Rip7 (5 
pa = pe/po P ip 
= / 
Va = VWe/ Vo ~ ye . ° 
an ») scosity-temperé  relz s: } 
62 = (6Re)/L? (5) Viscosity-temperature relation 
642’ = (0?’Re)/L — —— whi 
(u/VT) [1 + (C/T)] = const. 6) 
INTRODUCTION , : a roe ; 
where C = 203 for air, 7 in degrees Rankine. 
y i \HE PROBLEM ATTACKED HERE, the so-called ‘‘ther- 
mometer problem’ of compressible flow, had been PL = Pelle ( Rep 
previously solved with the help of finite differences by : ve , ae 
P ss z i ; ‘ Eq. (6) is the Sutherland empirical relation which the 
Fliigge-Lotz and Eichelbrenner,' but since their solu- . é ; *® 
‘ é ‘ ‘ nicely fits experimental data. Because of its awkward 
tion is laborious, the present work is offered as an : Ps : ae 
. : : : eee form, it will be used only to evaluate yu(7°) for the po- 
approximate solution that will give satisfactory accu- , ea os 
é ; as fis 5 tential flow just outside the boundary layer. Passing 
racy with a reasonable expenditure of time. The line aa 
: s “emia at constant x through the boundary layer, we shall -e 
of attack follows that conceived by von Karman. ; ; ——< rs ; sep 
; e Baa substitute the simpler Eq. (7). Crocco,* substantiated forn 
executed by Pohlhausen,” and refined by Holstein and ae ; : : 
ray oan se ; by Howarth,’ found that for the flat plate this substi- | 
Bohlen?’ for the corresponding incompressible fluid case, oe ; : 
ae é tution wrought no appreciable effect on temperatures 
Because of its length, the mathematical development me : . : 
: A and velocities, only a minor overestimating effect on 
can be presented here only in skeleton form. One may ‘ : ' 
; : : ; ee ie a drag. We assume the pressure gradient on our curved This 
find its detailed development and discussion in refer- : . , us 
surfaces will not magnify these errors appreciably. I 
ence 4. : plan 
n — 
BASIC RELATIONS IN THE BOUNDARY LAYER BOUNDARY CONDITIONS yo 
a ven iia , a : erab 
To solve the boundary-layer problem exactly, one [he boundary conditions for these differential equa : 
we ‘ aN 2 ; . - F subs 
must satisfy the following equations: tions are: that 
ae ; . a 
(1) Equation of conservation of mass: . a . . = 
(a) No slip at, or flow through, the wall—1.¢., Mw». = must 
[(0/Ox) (pu)] + [(0/Oy) (pv)] = O (1) Vs = 0 Al 
, _—— oe ae ; (b) Insulated wall, (07'/oy),, = 0 Eqs. 
Note that this consideration is limited to two-dimen- , Lie atas : 
, ma (c) Asymptotic transition of velocity and temper the , 
sional flow. These results can be extended to the three- se Mie ' 
; ; : ; ature from boundary layer to potential flow, u, > '; 
dimensional flow about axially symmetric bodies by foe . : 
; : Pee : T, => 1. 
the method of Mangler.* 
‘ ‘ It is assumed that the velocity u, (xv) is know! 
(2) Momentum equation: s cid ie . ; walk 
For a body of given shape, this function will be ca Ex 
[pu(Ou/Ox)] + [pv(Ou/dy)] = culated by the theory of frictionless fluid or else evalu- if we 
two | 


—(Op/Ox) + (O7/dOy) (2) ated experimentally. 
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BOUNDARY LAYER ALONG A Cl 


MATHEMATICAL DEVELOPMENT: SATISFYING AVERAGE 
MOMENTUM AND ENERGY RELATIONS 


The coupled, nonlinear differential equations [Eqs. 
1), (2), and (4)] force one to conclude that the prob- 
lem cannot be solved in exact form with a reasonable 
expenditure of calculation effort. The method of 
Fliigge-Lotz and Eichelbrenner, which is nearly exact, 
requires a calculation point by point, in both x and y 
increments, of the basic equations using finite differ- 
ences. If we take 


averages in the y direction and seek only to satisfy 


The work involved is prohibitive. 


the equations governing such averages, the y increment 
is thereby in effect extended to the entire width of the 
boundary layer and only x-wise steps need to be taken. 

The averaging is done by multiplying each term of 
Eqs. (2) and (4) by dy and integrating with respect 
to y through the boundary layer, using Eq. (1) to elim- 


inatev. The integration 
[ {Eq. (2)] dy yields 
7 0 


Tw / Pe = U0’ + u,.u,’0(2 — M.2 + (6*/0)] (8) 


The integration [Eq. (4)] dy yields 
70 


{° (i, — fi dy = 0 
rome | pu(je — j) dy = 


which for the insulating wall is 


| pu(j. — j) dy = 0 


0 


Replacing j by ¢ + (u?/2) and normalizing y through 


the substitution (first suggested by Dorodnitzyn) 


ba dy / ” dy ] 9 dy 
7 = =- (9) 
Jo pp) 0 p, 56/0 p, 


(1.e., p,6 dn = dy) we obtain the energy relation in the 


I 
- m [ u;x(1 — u,*) dn (10) 
/0 / V0 


form 
>] 
U,( p, 1) dn = 


This peculiar way of normalizing y requires both ex- 
planation and justification. The simpler substitution 


n = y/y. would give » = O when y = O, 7» > 1 when 
y > y., as does our substitution, and appears pref- 
erable because of its simplicity. However, such a 
substitution would yield approximate velocity profiles 
that do not change shape with Pr and M],, as they 
must (Fig. 1). 

Although Prandtl Number 


Eqs. (8) and (10), it does exert its influence through 


does not appear in 


the p, term, as will presently be seen. 


INTRODUCTION OF VELOCITY AND DENSITY 


POLYNOMIALS: ACCELERATED FLow 


Equations (8) and (10) become amenable to solution 
i! we can find expressions for u, and p, involving only 
For the accelerated flow 


two unknown parameters. 
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FIGURE | 
Effect of Mach Number and Prandtl Number 
on 
Flat Plate Velocity Profiles 
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i.e., from the stagnation point to the point of minimum 
pressure—let us assume 

ur, = Ao + a\n + den” + a3n?* + ayn! (11 
and apply the boundary conditions: 


At the wall satisfy (1) the condition of no slip (u = 
0) and (2) the differential momentum equation [Eq. 


(2) ]. 
At » = 1 simulate an asymptotic approach of u, to 
unity: (3) u, = 1; (4) Ou,/On = 0; and (5) 0°u,/On? = 


0. Eq. (11) then becomes 
u, = n(2 — 2n*? + n*) + (An/6) (1 — n)*-~= (12) 

where 
A = [(06?)/v, Ju,’ (13) 


It is the fundamental parameter of this profile. 


At Pr = 1, 
which yields the relation 


j = je ort + (u?/2) = 4. + (u,?/2), 

pp = 1+ m(1 — u,”) (14 
Thus we require an eighth-degree density polynomial 
to be consistent with u,. 


pr = b+ dn +... + den (15 


Since we do not wish to specify eight boundary condi- 
tions, we also write it in alternate form. 


+ ml — 4,7) + 
(1 — Pr) (d + din - 


p= 1 


+ dsn°’) (16) 


and apply the following boundary conditions: 

At the wall satisfy (1) the condition of no heat trans 
fer; (2) the differential energy equation [Eq. (4)]; and 
(3) 0/Oy [Eq. (4) ]. 

At 7 = 1 simulate an asymptotic approach of p, to 
unity: (4) p, = | and (5) Op,/On = 0. 

Eq. (15) becomes 
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(6b — 1) (1 — 5n* + 4n*) + 
PrX }[AC12 + A)]/36} n® (1 — 9)? 
Prm [(12 + A)/36] n? (1 — 9)? X 
(1211 + 29) + ACL — 2n)] + 
m} [n'(1 — »)?]/36f [86(1 + 2n — 9?) + 
12\(2 — n) (1 — 9) — A*X(6 — 40 + 7’) ] 


(17) 


The density (or temperature)* profile is seen to depend 
upon the following variables: Pr (constant for a given 
fluid); m = |(y — 1)/2]M,? (compressibility variable) ; 
\ (basic profile parameter); and X [X = (b’u,)/(bu,") |. 
This last variable, Y, exerts through }’ a controlling 
influence on the variation of temperature along the 
wall. When Pr = 1, X = 2m, showing X to be another 
manifestation of compressibility effect. 

At this point a few words need to be said about the 
choice of boundary conditions for the u, and p, poly- 
nomials. 

Experience in the comparable Pohlhausen method for 
incompressible fluid has shown that a different choice 
of boundary conditions (e.g., impose O*,/On* in addi- 
tion to 0°u,/O0n? = O, or in lieu of it at » = 1) may 
greatly change the predicted location of the separation 
point and, to a small degree, reshape the velocity pro- 


AERONAI 


STLcAaL 


SCIENCES POUL y; Vis 

wall (6b = py/p. = T./7T.) which appears essential but 
is not otherwise obtained. Note that no attempt is 
made to assign separate thicknesses to the thermal and 
kinetic boundary layers. In each case, it is assumed 
that temperature and velocity attain potential value 


at the same point, 7 = 1. 
three outer boundary conditions, is required to ap 


However, velocity, with 


proach its potential value more rapidly than temper 


ature, with only two outer boundary conditions 
After all, an artificiality must be introduced when 
finite ‘‘thickness’’ 


actually asymptotic in profile, and the above device 


is assigned to the layers that are 


is at least as reasonable as the carrying of separate 


thermal and kinetic boundary-layer ‘‘thicknesses, 


and much simpler mathematically. 
VELOCITY AND DENSITY POLYNOMIALS: DECELERATED 
FLOW 


On the basis of our studies of choice of boundary 
conditions it was concluded that a better prediction of 
separation point would result if for decelerating flow 
from the point of minimum pressure to the separation 


point—we change to a different velocity polynomial: 


files. Therefore careful attention should be given to 
this point, and boundary conditions should be chosen Uu, 
which lead to best agreement with those exact solu- 
A study by Mangler* on this 


= ant... t asn° (18 


One additional boundary condition is imposed: Satisfy 
0/Ooy [Eq. (2)] at the wall.f Putting in the boundary 


conditions, Eq. (18) becomes 


tions that are available. 
point augmented our own calculations and compari- 
sons and led to the choice indicated. The experience of 
Fliigge-Lotz and Eichelbrenner on this same problem 


‘ : : , ie : u, = (n/3) (5 — Jn? 3n*) + 
led to the inclusion of the third condition for the p, ”o Tt pa 4 (1 14 19 
° 7 = 3 ) —7)° (1! 
profile }0/Oy [Eq. (4)] satisfied at the wall}. It brings mn ie 1 
our variable Y, thence b’, into the solution and thereby 
t Imposing this additional boundary condition on the acceler . Soe ’ 


gives a control of temperature variation along the ; 
ating flow is found to decrease accuracy there, contrary to what 


* Since Op/Oy = 0, pr = Ty. one might expect. 





The 
TABLE 1 not cl 
1 y 
For Pr = 1 aye 
r K f fy + fs +f 2(0/6 Jay aa 
12 0.0948 —(0.095 1.043 0.711 ‘ 
11 0.0941 —0.091 1.038 0.709 
10 0.0919 —(. O80 1.023 0.7038 
9 0.0882 —(0.061 0). 997 0.693 
8 0.0831 —(), O34 0.963 0.679 
7.3 0.0819 —(),027 0.956 0.677 
7.6 0. O807 —(0 020 0.948 0.674 
7.4 0.0794 —().013 0.940 0.670 
= 0.0781 —0 0086 0.932 0.666 
rf 0.0767 0.002 0.923 0.662 
6.8 0.0752 0 O10 0.914 0.659 
6.6 0.0737 0.019 0.904 0.655 
6.4 0.0721 0.028 0.893 0.651 Eqs. 
6.2 0.0705 0.037 0.882 0.647 
6 0.0689 0.046 0.872 0.643 revert 
5 0.0599 0.098 0.816 0.620 b—] 
} 0 0497 0.158 0.754 0.595 
3 0.0385 0,225 0.688 0. 567 
2 0.0264 0.299 0.618 0.536 SIMPI 
1 0.01385 0.381 0.545 0.504 
0 0 0.470 0.470 0.470 
0 0 0.414 0.414 0.414 The 
= —() 0157 0.506 0.316 0.355 
= —().0316 0.604 0.225 0.293 substi 
—3 —().0469 0.704 0.148 0.229 
=~ —0 0610 0.804 0.072 0.165 de? 
—5 —0 0735 0.900 0.019 0.101 
—6§ —().0835 0.987 —(0.016 0.039 





il but 
Ipt is 
il and 
umed 
value 
with 
) ap 
nper- 
tions 
when 
t are 
levice 
arate 


sSes, 


ATED 


idary 
on of 
OW 

ation 


al: 


(18 


itisfy 


dary 


(19 


‘celer 


what 





BOUNDARY LAYER ALONG A CC 
TABLE 
For Pr = 
A kK } Q 
10 0.0919 —() OSD 0 3S4 
g 0 O882 —() 061 0.389 
8 0.0831 —O0O 034 0.395 
7.8 0.0819 —(). 027 0.396 
7.6 0. OS8O7 —() O20 0 298 
7.4 0.0794 —() 013 0 399 
42 0.0781 —() O06 0.400 
7 0.0767 0. OO2 0.401 
6.8 0.0752 0.010 0.403 
6.6 0.0737 0.019 0.404 
6.4 0.0721 0.028 0.406 
6.2 0 .O705 0 O37 0. 407 
6 0 0689 0.046 0.408 
5 0 0599 0.098 0.416 
it 0. 0497 0.158 0.425 
3 0. O385 0.225 0.434 
4 0. O2b64 0. 299 0.445 
l 0.0135 0.381 0.457 
0 0 0.470 0.470 
0 0 0.414 0.414 
— |] —() OLST 0.506 0.413 
—_ 2D —() 0315 0.604 0.414 
—3 —() 0468 0.704 0.417 
—3.5 —0.0540 
—4 —() 0610 O.S804 0.423 
—4 3 —() 0636 
—44 —(). 0662 
—4 6 —() 0687 
—+4.8 —() O7F11 
— 5 —() 0734 0.900 0. 431 
—5. 1} —(). 0745 
—5.2 —() O756 
—5.3 —(). 0767 
5 —() O778 
—5.5 —() O788 
—5.6 —(). 0798 
—5 7 —() O8OS8 
se —(0.O0818 
—5.9 —(). O827 
—6 —(). O836 0.987 0.441 
—6.05 —0._ 0840 
—6.1 —(). 0844 
—6.15 —(). 0849 
—6.2 —(0.0853 
—6.25 —().O857 


See Table 1 for 2(0/6)a 


The boundary conditions for the density profile are 
not changed, but, since it must now be a tenth-degree 
polynomial, several terms are altered, giving 


(6 — 1) (1 — 5m + 4°) + 
PrX [X(20 + 3A)/72] n3(1 — 9)? — 
Pr m[(20 + 3A)/12] n?(1 — n)* X 
[(5/3) (1 + 29) + (A/12) (8 — 2n)] + 
m }[n4(1 — »)2]/144} [—16(30 — 4y — 
Sn? — 12? + 9m) + 240(5 + Sn — 
3n? — 5n*® + 3n*) + 9A2(1 — 4n + 
n° + 2n* — n*)] 


p,-1= 


(20 


Eqs. (17) and (20) may be checked by noting that they 
revert to Eq. (14) when we let Pr = 1, for which case 
b-~1=NX/2=~m. 

INTEGRATED MOMENTUM AND 


SIMPLII YING THE 


ENERGY EQUATIONS 


The velocity and density polynomials may now be 
substituted into Eqs. (8) and (10) to yield 
de°/dx = [v,/(bu,)] [fi + (6 — Ife + 
m( fs + Pr f;)] 


(21) 


RYED INSULATED SURFACE $49 
0.725 

O.725f B AB 

0.626 10.32 67.7 

0.596 7 56 70 4 

0.558 t SU ia.a 

0.550 4.25 73.9 

0.541 3.70 74 5 

0.532 3.14 75.2 

().522 2.59 75.9 

0.512 2 OA 76.6 

0.501 1.48 i 

0 490 0.93 72.0 

0.479 0 38 78.7 

0.468 —() 18 79.4 

0.457 —() 73 8) 2 

0.395 —3.51 84.1 

0.326 —6.28 SS_6 

0.251 —9 06 93.7 

0.171 —11.85 99.4 

0. O87 — 14.64 106.0 

0 —17.48 113.7 

0 — 14.45 132.2 

—0.095 — 16.50 153.8 

—(0.189 — 18.59 IS1L.0 

—0.275 20.68 225.1 
57 

—0.350 —29 79 01.9 

324.7 

R51 5 

383.4 

22.1 

—0.412 — 24.92 $170.2 

{98 8 

531.3 

568.6 

611.7 

662.2 

722.2 

794.5 

883.6 

995.9 

—0.456 —27.07 1141.8 

1232.5 

1339.2 

1466.6 

1621.3 

1813.0 

db ‘dx = (bm/u,) (u,’/) [B + (1 — R)AB] (22) 


where f;(A), B(A, Pr), and AB(A, Pr) are known func- 
tions that have been calculated and tabulated in Tables 
1, 2, and 3. 

These equations will be of more general application if 


all terms are rendered dimensionless: 


6," /dxa = [vg (bua) | [fi + (b = 1) fs + 
m (fs + Pr fs)] (23) 


db/dxq = (bmua’/uad\) [B + (1 — R)AB] (24) 
In these equations ug, %q’, va, and m are known from 
Variables \, 6,, and 0} are to be eval- 
One additional relationship 


the potential flow. 
uated as functions of xq. 


is needed. It arises out of Eq. (13): 


02/52 = 02u,'b/v. = k (25) 
1 

Since 6/6 = f (u, — u,*) dy, the group [A(@?/6?)] is a 
0 

function of \ only. We designate it « for convenience. 

In dimensionless form, 


Kk = 047ta’b/ va (26) 
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TABLE 3 


For any Pr 


A J ds g 
10 0.590 0.049 45.51 
9 0.565 0.043 45.51 
8 0.531 0.087 45.52 
7.8 0.523 0.036 $5.53 
7.6 0.515 0.035 $5.53 
4 0.507 0.034 $5.53 
7-3 0.498 0.033 45.54 
‘4 0.489 0.032 45.54 
6.8 0.479 0.031 15.55 
6.6 0.469 0.030 45.55 
6.4 0.459 0.028 $5.56 
6.2 0.448 0.027 45.56 
6 0.437 0.026 45.57 
5 0.379 0.020 $5.61 
} 0.314 0.015 15.65 
} 0.243 0.011 S74 
2 0.166 0.007 415.78 
l 0.085 0.008 $5.85 
0 0 0 45.94 
0 0 0 38.05 
—] —(0).094 —(0. 002 36.23 
2 —(), 187 —(0.008 34.45 
—3 —().273 —(0. 003 32.69 
—4 —0.349 —(0.002 30.98 
a —0.411 —0.001 29 33 
—6 —0.456 —0. 000 > i 3 
B = [2X — (1 — Pr)gs|/Pr. 


AB (Pr = x) = (0.725/x+)AB (Pr = 0. 
See Table 2 for f;, fs, AB (Pr = 0.725). See Table 1 for 2(0/6)a,. 


725) 
pee 


METHOD OF SOLUTION 


The solution of a given boundary-layer problem may 
now be made by the following steps: 

(1) Determine u,(x) for the body in question (evalu- 
ate experimentally or calculate by the theory of fric- 
tionless flow). 

(2) Evaluate reference values 7) and um at some se- 
lected reference point in the flow, such as the undis- 
turbed free stream. 

(3) Express the velocity “,(x) in dimensionless form, 
Ua(Xq). 

(4) Calculate the constants 

Ci = 203 T» 
ky = 2cyTo/Uo? = 1/mo 


where c, is expressed in (ft. Ib.) (slug)~! (°R.)~'. The 
value 5,920 is suggested for air. 
(5) Calculate the potential flow quantities 


m = ua?/(1 + ky — UQ”) 
Ta = (1 + ky 
va = (1/Ta) [A + C)/(Ta + Ci) ] 


Ua") ‘Ry 


(6) Evaluate the initial values for 5, by’, », K, 64” 
6°’. If the calculations are initiated at the stagnation 
point, these will be found to be (for details see refer- 
ence 4, p. 84-86) 


b = | 

by’ = @ 

d = 7.05 - 
« = 0.0770 wre 
6a” = 0.0770 (v¢/ua’) 

647’ = —0.0653 [(vaua")/(ua’)?] 
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where the values of vag, ua’ = Og/OXa, and nu," = 
O°u4/OX," are those at the stagnation point. 

If the calculations are initiated at some other point, 
it may be necessary to have—presumable from exper; 
ment—velocity and density profiles “(y) and p(y) at 
two points closely adjacent (a distance Ax apart) from 
which to calculate these initial values. Use the de. 


fining equations 


b = pe/ pu 
ba’ = OD OXg 


bé*ue’ ; a 
r where 6 = dy 
p 0 p 


\(0/6)? = A[(37/315) — (A/945)- 98) 
(A2/9072) |? 


P Uy f° pu (: u ) ; 
d = — ay 
\ wLJ0 pelt, uu, 


16,7 /dXq 


II 


II 


> 
| 


For supersonic flows around airfoils with sharp nose 
there will, of course, be no stagnation point. It will 
probably be satisfactory in such cases to use the initial 
values) = 1: 6, = O0:A = « = 6, = 6,” = 0. 

(7) Choose an arbitrary small increment Ax, and cal- 
culate 9,” and 6 at the station x, (initial) + Ax, by 
finite-difference equations of the form 

s(w + Aw) = s(w) + 2’(w)- Aw 

Greater accuracy will be obtained if the average value 
of by’ (and 6,2’) over the increment Ax, is used instead 
of the initial value as indicated in the above equation. 
This, however, involves some trial-and-error calcula- 
tions since },’ and 6,7’ at x, + Ax, are not known until 
calculations for the given increment of Ax, are com- 
pleted. 

(S) Calculate at x,y (initial) + Ax, the remaining 
variables x, 6,°’, and b,’ from the equations 


K = O47Ua’b/ v4 
64°’ = [va/(bua)] [fi + (6 — 1)fe + m(fs + Pr fi)] 


where /; _4 are functions of « or A, tabulated in Tables 1, 


2, OF 33 
ba’ = [(bmug’)/(ugd)] [B + (1 — R)AB] 


where B and AB are tabulated in the same tables. 
(9) Calculate the wall temperature 7), from 


= TaT vb (29a 


(10) Calculate the local skin-friction drag coefficient 


from 
CaV Re = 2(0/5)ai(va/Oa)UaT a (29b 


where 2(6/65)a; is a function of «x found in Table 1. 

All quantities are now known at station xy + 1%. 
Choose additional increments in x, and repeat steps 
7-10 for each. See the end of this paper for examples 
of these calculations. It will be observed that at the 
point of minimum pressure (x, = 1.67 in the Fligge 
Lotz example) there occur discontinuities in the func- 
tions f, fz, AB, and B, causing slight irregularities in 
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the values of 04°’, ba’, and Cav/ Re. A change in the 
form of the velocity and density polynomial at the 
minimum pressure point causes these discontinuities. 
The computer should terminate a Ax, increment at the 
minimum pressure point so that he will not need to 
interpolate tabulated values of fi, fo, ete., over an in- 
terval including these discontinuities. Otherwise this 
point requires no special treatment. 

(11) Locate the separation point at the value of x, 
for which « = —Q0.0890. Calculation difficulties at 
the separation point make it expedient to terminate our 
calculations short of the separation point and then 
locate X¢ «wep. by extrapolation of the curve x¢(«). 

Since only the wall temperature is appreciably af- 
fected by the value of the Prandtl Number (see p. 453) 
if this temperature is not needed the calculations may 
For this 


be made under the assumption that Pr = 1. 
mathematically simpler case, b = 1 + m and can be 
calculated a priori, obviating the calculations of b,’ and 
bin steps 6, 7, and 8. Step 9 is eliminated completely 
since 7, = ZT» (1 + mo) = const. for Pr = 1. 


EXAMPLES AND SUMMARY 


The thinness of the boundary layer makes it almost 
impossible to obtain experimental data sufficiently 
wecurate and definitive to afford a means of checking 
the theory. Velocity profiles can be measured, but 
it is extremely difficult to measure temperatures accu 
rately. We can only test our results against those of 
more mathematically exact theory insofar as it exists. 
The flat plate case has been satisfactorily solved by a 
number of workers. Crocco has done extensive work 
m this case. For the curved surface there is the 
Fliigge-Lotz, Eichelbrenner report. Due to the lengthy 
computations entailed, they have calculated essentially 
one example, for the case Pr = 1, although there are a 
few results for Pr = 0.7. We shall compare results 
with these two cases. Beyond that our calculations 
must stand alone until other exact solutions appear for 


further check. 


The Flat Plate 


The flat plate is assumed to be extremely thin and 
set at zero angle of attack so that the pressure gradient 
It follows that ug’ = \ = « = 0, uw = 
Tl; = v = 1,m = const. Evaluate Eq. (23) at \ = 0 
noting from Tables 2 and 3 the values f; = fs = 0.470. 

= /,= 0: 


will be zero. 


d6_7/dx~a = 0.470 

Integrating with respect to xy we obtain 

64° = 0.470 xa + const. 
Where the constant is zero since 6, = 0 when x = 0, 
measuring x, from the beginning of boundary-layer 
lormation—i.e., from the leading edge of the plate. 
Evaluate Eq. (24) at \ = 0, using Tables 2 and 3: 
db/dxg = bm(ug’ Ugr) | } [45.94(Pr — 1)]/Pr} + 


11 — [(b — 1) m]! [(0.725/Pr) 113.7]) (30 
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The ratio u#q@’/A = 0/0 is evaluated through its alter 


nate expression [cef. Eq. (25) } 
up'/X = (n/X) [v, (A°b 


or in dimensionless form, using also Eq. (28), 
Ua’ /X = (K/A) [ve/(Oq2b)] = (37/315)? [1/(0.470x4d) | 
31 


Eq. (30) is a simple differential equation in the form 
db/dxa = (1/Xqa) (bke + ks 
for which the only finite solution is 
b = — ks/ko = 1 + m (0.4428 + 0.5572 Pr 


Observe that b = 1+ mfor Pr = 1, asit should. The 
wall temperature comes from this relation since 


b = Pp, p : ye ‘i 
At Pr = 0.725 (air), 
T, = T, (1 + 0.8468 m 
Crocco finds the exact relation to be 
Tl, = T, i + 0.8513 


Our approximate theory is here in error by 0.2 per cent 
for 7, = 2,orm = 0.8. If A/, is lower and, or Prandtl 
Number is nearer unity, the error will be still less. 

The local drag coefficient is calculated, with the aid 


of Eq. (29b), to be 
Ci Re = 0.6854 V1 
The total drag coefficient is then 


| an 1.371 
C/o Co VW Re dx 
Ya JO V Na 


Crocco gets CpV/ Re = 1.328/ VV xq. It is interesting to 
note that this value is the same as that obtained by 
Blasius for the incompressible case. In other words, 
the drag of the insulated flat plate at zero angle of 
attack is unaffected by Mach and Prandtl Numbers. 
Our approximate theory overestimates drag by 3.2 per 
cent. 

To compare velocity and temperature profiles, use 
the following relation which derives from Eqs. (9) and 
(13), 

Ava i, 
% = | dy (32 
0 


bu," 
Using Eas. (31) (17), and (32) we obtain 
g hq 


va = 5.84 [n + m((0.4428 + 0.5572 Pr) X 


tn — o + (2/3)n°} — 4 Pr} (n°/3) — [(3/5)n*] + 
(n°/3)} + (n° 5) — [(4/7)n7] + 
(9°/2) — (9° 9))| 


Combined with Eqs. (12) and (17), this evaluates 
u;(y¢) and p,(y¢). Calculations have been made for 
air, Pr = 0.725, and the results compared with Crocco’s 


profiles in Fig. 2. 
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Flugge-Lotz Example 
In the Fliigge-Lotz-Eichelbrenner report, one ex- 
ample—for Pr = 1 The authors 


chose a velocity distribution that would incorporate a 
flows over 


had been completed. 


number of the characteristic features of 
yrodynamic bodies. The velocity starts from a stag- 
nation point, has a region of accelerated flow, reaches 
, maximum velocity, and then has a region of deceler- 
ited flow which persists until separat on of the bound- 
ary layer occurs. The shape of body over which this 
velocity distribution would occur was not ascertained. 
In that report the reference values are those at the stag- 
nation point, with reference velocity m taken as the 
speed of sound in the stagnant fluid. 


We shall make a number of calculations based on this 
example and variations of it. It is thought convenient 
to reference the example to the free-stream conditions. 


he velocity is then given by the analytic expression 


Ua = [Xal(xa + 5)]/(x%a + 1)? 
and the reference values are % = 997 ft. per sec.; 
Ty) = 435°R.; and » = 2.091 X 10~‘ sq.ft. per sec. 


By means of our calculations we shall show for this 
example (1) accuracy of our approximate method; 
2) effect of compressibility; (3) effect of Prandtl Num- 
ber; and (4) effect of velocity gradient u,’. 

(1) Accuracy of Method. 
first calculation is made for conditions identical with 
those taken in the calculation. Since 
the latter may be considered to be exact, a comparison 


Our Approximate Our 


Fliigge-Lotz 
of the results serves to measure the accuracy of our 


In Fig. 3 the drag and momen- 
Also sepa- 


ipproximate method. 
tum thickness are compared (Pr = 1). 
ration is located at the point where the drag drops to 
zero. The momentum thickness is about 10 per cent 
low. The separation point is 31 per cent too far aft, 
ind this is largely responsible for an overestimation 
of drag over the major portion of the curve. 

A comparison of velocity and density (or temper- 
ature) profiles was also made at approximately the 
mid-region of accelerated flow and of decelerated flow. 
Figs. 4 and 5 portray these results. There is also 
added the velocity profile for incompressible flow, cal- 
culated for the same velocity distribution. Temper- 
iture is, of course, constant in the incompressible 
flow, and the (7, — 1) line would coincide with the y 
axis. It is seen that an extremely satisfactory accu- 
racy is obtained in the accelerated flow. The approxi- 
mation is somewhat rough for decelerated flow, al- 
though our method is set up so that the wall value of 7, 
is given exactly [at b = (7,)~ = 1 + m] for the case 
Pr = 1. It will be almost exact for other values of 
Prandtl Number. 

Local drag per unit area, tr, = [u(Ou/dy) Jo, will 
depend on (7,),,, which fixes y,,, and on (Ou Ov)». The 
good reproduction of (7°), means that if there is any 
maccuracy of drag prediction it is due to the error in 
the velocity profile slope at the wall. 


IRVED INSULATED SU 


RFACE 153 


(2) Effect of Compressibility.—Some effect of com- 
pressibility may already be seen in Figs. 4 and 5. The 
velocity profile of the compressible flow is changed 
appreciably from that of the incompressible flow. 
The boundary layer becomes thicker in the studied 
cases. In order to study further the influence of com- 
pressibility, new examples were added and the vari- 
ables that characterize the boundary layer (@,, Cp, and 
Two = T,,/To) were plotted. 

First the basic example for Pr = 0.725 was calculated. 
In addition, the comparable incompressible case was 
calculated after setting m = 0, 7, = vg = 1. A third 
case was formed by increasing the reference velocity 
by approximately 17 per cent. Dimensionless ve- 
locity uq(xq) remained unchanged in the three cases, as 
did reference temperature 7). Fig. 6 presents the re- 
sults. It is seen that every variable compared feels a 
major compressibility influence. We recognize, of 
course, that this is not a comparison of boundary layers 
over the same body traveling at different speeds, since 
the form of the velocity variation ,(xqg) would be 
altered by a change of speed of travel. However, the 
comparison is certainly of qualitative value. 

It may be added that increasing the free-stream 
Mach Number by 17 per cent—from 0.982 to 1.15 
while holding w4(xa) the same was accompanied by an 
increase in maximum local Mach Number from 1.80 
to 2.27, a 26 per cent increase. 

(3) Effect of Prandtl Number. 
amples were recalculated at Pr = | (except for the in- 
compressible case where Prandtl Number has no sig- 
The re- 


The preceding ex- 


nificance and does not enter the calculations). 
sults of Crocco show that greater Prandtl Number ef- 
fect will be encountered at the higher speeds. There- 
fore, we have compared the cases, in Fig. 7, at the 
higher Mach Number to magnify any Prandtl Number 
effects. It is seen that there is no appreciable alter- 
ation of momentum thickness, drag, or separation point 
by going from Pr = 1 to the value for air, Pr = 0.725. 
In addition, two completely supersonic flows (not shown 
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in this paper) have also been calculated—shock-wave 
effects have been ignored—for Pr = 1 and 0.725. They 
further support the conclusion that the Prandtl Number 
affects only the temperature to an appreciable degree. 
This is a pleasant discovery, since it means that for 
most purposes we may make our calculations at Pr = 
|, which is a much simpler and less tedious calculation 
than for Pr ¥ 1. Prandtl Number does alter wall 
temperature ratio 7’,,. At Pr = 1 the ratio remains 
constant at 7) = 1 + m. At Pr < 1 this ratio has 
the same value at the stagnation and separation points 
but dips to lower values in between. We must not 
conclude from the small percentage wise variation of 
7... that here again the Prandtl Number effect will be 
of no practical importance. In the example shown the 
minimum value of 7. is only 7.4 per cent below its 
stagnation value, but a 7 per cent decrease could 
amount to 70° for AJ, = 4, say. This would be of 
considerable importance in calculating the temperature 
of a pilot’s canopy. 

(4) Effect of Velocity Gradient uy’.—-It is readily 
seen that increasing the reference velocity while hold- 
ing velocity ratio 4, = u,/u) unchanged results in an 
increase in “u, at a given value of x and therefore cor- 
responds to a numerical increase in u,’ at all points 
except the minimum pressure point where it is zero. 
Thus the method employed above in ascertaining com- 
pressibility effects is as much a measure of the effect of 
the slope of the u,(x) curve as it is of compressibility. 

A more clear-cut manifestation of slope effect will be 
found in holding the reference values constant while 
increasing the value of u,. Partial examples were made 
for the two cases: 


ug = 1.15 [xu(5 i Xa) |/ (Xa = 1)?! 


for 0 < x; < 1:2. 


Ug = 1.25 } [xal(5 + xa) )/ (va + 1)? — 0.390 


for 1.67 “9 Xa SJ (Xa)sep.- 


of an increase in u,’ in the accelerating part of the 
boundary layer. The second reveals the influence of 


The first shows the effect 


an increased deceleration. To give a smooth transition 


from the basic w(x) curve, this latter alteration was 
begun at the point of zero slope u,’ with the same values 
of uq and 6, which were found for the basic case. Fig, § 
shows the results and reveals an extremely strong slope 
influence in the region of accelerated flow. The sepa 
ration point is materially advanced by increased ad 


verse pressures. 


LIMITATIONS 


(1) No attempt has been made to account for bound 
ary-layer shock-wave interaction. 

(2) Under conditions of high positive velocity gra 
dients, or sustained gradients of moderate magnitude 
the proposed method may break down both mathe 
matically and physically. The formal calculations 
may require a value of «x greater than 0.0948. But, 
reaches a peak value at 0.0948 (A = 12) and further 
computations cannot be made. Physically, values 
of X > 12 are unsupportable because they will give 
boundary-layer velocities greater than the exterior 
velocity at the same cross section of flow. Neverthe 
less, while some practical cases may not be amenabk 
to solution by the above method, there are many more 
which can be solved with good accuracy and reason 
able expenditure of time. 

COMPARISON WITH THE METHOD OF (GRUSCHWITZ 

Gruschwitz® studied the same problem at approxi 
mately the same time but published it earlier. He 
assuines different thicknesses for the velocity and the 
temperature layers, and, by an ingenious method, he 
avoids the trouble that usually is caused by employing 
two different thicknesses. He takes as his polynomials 
u, = filn) and p,(1 — u,) = fo(m). The latter poly 
nomial is required to satisfy the outer boundary con 
ditions, fo = Of2/O0n = O2f/2/0n? = 0. But it satisfies 
these by virtue of the character of u, and places no re 
quirement upon p, to reach unit value when w, does 
Thus he can make allowance for the difference in ther 
mal and kinetic thicknesses while carrying just one 
thickness through his calculations. There is a serious 
disadvantage to this method, however: He loses con 
trol of the ratio p, and has no assurance it will reach 
unity or its slope will become zero at any point. Cal- 
culations of his density profiles were made for several 
reasonable cases, and it was found that p, is not well 
approximated; it can even become negative near the 
outer edge in some cases. The fact that in his second 
polynomial the region of objectionable values of p, 0¢ 
curs where 1 — uw, is extremely small minimizes its 
harmful effect on the rest of his results. In fact, he 
achieves in the region of accelerated flow better agree- 
ment with the exact Fliigge-Lotz example than do 
the present authors. Considering the renegade action 
of p,, one cannot help feeling that the agreement might 
be only a coincidence peculiar to this example. Cer 
tainly its radical departure from physical reality forms 


(Continued on page 490) 
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On the Stability of the Laminar Boundary 
Laver in a Compressible Fluid’ 


D. W. DUNNT ann C. C. LIN} 
Massachusetts Institute of Technology 


ABSTRACI 


A more complete theory is developed for the stability of the 
laminar boundary layer in a compressible fluid, removing some 
limitations and modifying certain general conclusions in the 
earlier treatment of Lees and Lin. The main changes are the 
inclusion of three-dimensional disturbances and the demonstra 
tion that, in many cases, the stability characteristics depend 
m temperature fluctuations (contrary to previous conclusions) 
The former change is made possible by a more careful analysis 
iid simplification of the basic equations of the disturbances; 
the latter modification is found through improved solutions of 
these equations 

At high supersonic free-stream velocities, the relation for the 
haracteristic values is found to be dependent on temperature 
fluctuations and the boundary conditions imposed on them, ex 
ept for three-dimensional disturbances traveling very obliquely to 
the direction of the free stream. For subsonic and slightly super 
sonic velocities, where the wave speed of the disturbances is 
usually not large, the earlier conclusion of Lees and Lin is still 
valid—i.e., temperature fluctuations are unimportant tor the 
stability problem 

For the boundary-layer flow over a flat surface, three-dimen 
sional disturbances are found to be of little importance at low 
subsonic free-stream Mach Numbers, since they usually have a 
higher minimum critical Reynolds Number than two-dimensional 
ones, except possibly under conditions of extreme surface cooling 
As the Mach Number increases, three-dimensional disturbances 
become significant under conditions that are less and less ex- 
treme, until finally, at a Mach Number between one and two, 
they begin to play the leading role in many cases of practical 
interest 

At supersonic free-stream Mach Numbers the boundary layer 
in never be completely stabilized by cooling with respect to all 
three-dimensional disturbances. There is always present a class 
of unstable three-dimensional disturbances no matter how low 
the temperature of the solid surface becomes. However, in 
many cases, surface cooling is still a very effective means of 
stabilizing the boundary layer 
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Numerical results illustrating these conclusions are presented 


for the case of the boundary layer on a flat plate 


(1) INTRODUCTION 


I RECENT YEARS there has been considerable interest 
in the theory of boundary-layer stability for com- 
pressible fluids. Most of the recent work on the sub- 
ject has been based on the original theoretical investi- 
gations of Lees and Lin.! As has been pointed out 
on several occasions, there are certain inadequacies in 
this basic theory which limit its significance and which 
in some cases make unreliable the results of numerical 
calculations based on it. The purpose of the present 
investigation is to clarify some of these problems and, 
in general, to gain further insight into the stability 
theory for the laminar boundary layer in a compress- 
ible fluid. 

One of the chief limitations of the theory up to now 
has been consideration of only two-dimensional disturb- 
ances by analogy with the procedure for an incompress- 
ible fluid. 
mation of the equations of the problem that for an 


Squire* has shown by means of a transfor 


incompressible fluid every three-dimensional disturb- 
ance is equivalent to a two-dimensional one at a lower 
Reynolds Number. Therefore, for an incompressible 
fluid the minimum critical Reynolds Number for two 
dimensional disturbances is the minimum critical Rey- 
nolds Number for all disturbances, so that the restric- 
tion to two-dimensional disturbances is justified. If 
an attempt is made to prove a similar result for a com- 
pressible fluid, it is found that the equations obtained 
by making the usual approximations (see reference 1) 
are too complicated to permit a transformation of the 
type used by Squire. In the present work an order- 
of-magnitude analysis of the terms in the complete 
equations for a three-dimensional disturbance is car- 
ried out. A consistent approximation scheme based 
on this now yields simplified equations of the same 
order of accuracy as the usual ones, but with fewer 
terms, so that a transformation to equations of the 
disturbance is 


form as for a two-dimensional 


However, these transformed two-dimensional] 


same 
possible. 
equations are not the equations of a proper two-dimen- 
sional disturbance, so that no direct conclusion regard- 
ing the importance of three-dimensional disturbances 
can be made. For a compressible fluid the properties 
of three-dimensional disturbances have to be found by 
means of a more elaborate investigation than for an 


incompressible fluid. The transformed two-dimensional 
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equations lead to a considerable simplification math 
ematically, however. 

It is found in the present investigation that for sub- 
sonic free-stream Mach Numbers three-dimensional 
disturbances are not of much importance, but for 
supersonic Mach Numbers they may play a leading role 
in certain circumstances. In particular, for supersonic 
sonic Mach Numbers there is always present a class of 
unstable three-dimensional disturbances no matter 
how low the temperature of the solid surface is, so that 
complete stabilization of the boundary layer in such 
cases, as predicted in reference 2, is impossible in the 
strict sense. 

It is known that at high Mach Numbers the usual 
form of the boundary-layer theory becomes inaccu- 
rate," ° so there is some question as to the proper form 
of the stability theory in such cases. An investigation 
of this problem® has shown that, as long as the proper 
velocity and temperature profiles for the mean flow 
are used, the simplified stability equations derived in 
this paper are still valid as a first approximation at high 
Mach Numbers. 
of a first approximation may not be great, so that higher 


However, in such cases the accuracy 


approximations and thus more complicated equations 
may have to be considered. t 

Various limitations to some of the approximations 
made in the theoretical developments of reference | 
have been brought out. In fact, as pointed out in refer- 
ence 7, the original theory of reference | has sometimes 
been used in situations in which it was not intended to 
be valid. This theory was formulated with the restric- 
tion to small values of the wave speed c, so that, strictly 
speaking, its range of reliability is limited to values of 
the Mach Number not too much larger than unity. In 
particular, the viscous solutions of reference | are accu- 
rate first approximations only in the immediate neigh- 
borhood of the point where w = c. However, at high 
Mach Numbers, c is often large in the stability calcula- 
tions, so that, in order to satisfy conditions at the solid 
boundary, the viscous solutions have to be accurate far 
away from this point. For an incompressible fluid, 
Tollmien* obtained viscous solutions that are uni- 
formly accurate all across the boundary layer. This 
work has been extended to the case of a compressible 
fluid, and a new mathematical formulation of the sta- 
bility problem in terms of the new viscous solutions 
has been obtained. In the present paper the modified 
formulas for the calculation of the characteristic values 
are given. The detailed mathematical discussion will 
be presented elsewhere. 

The conclusion reached in reference |—-that the rela- 
tion for the characteristic values is independent of the 
temperature fluctuations and the boundary conditions 


imposed on them—is now found to be restricted to sub- 


t As an example of the increase in complication, in a second 
approximation terms must be retained whose coefficients involve 
the v-component of the basic flow velocity, so that to this order of 
approximation the flow cannot be considered to be parallel 
This conclusion was reached independently in reference 6 and in 


an investigation by Cheng. '” 
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sonic and slightly supersonic Mach Numbers in gen 
eral. At high Mach Numbers it is true only for three- 
dimensional disturbances whose direction of propaga- 
tion is at a large angle with the free-stream direction 
In all other cases the form of the relation for the char. 
acteristic values depends on the boundary condition for 
the temperature fluctuations at the solid boundan 
However, for very small values of the wave number, 
the relation again takes on a simple form much 
like that of the earlier theory. This situation, which 
corresponds to very low frequency fluctuations, occurs 
in the limiting case of complete stabilization of the 
boundary layer by cooling, for example. Thus, it 
is still possible to obtain the critical temperature ratios 
for complete stabilization by fairly simple numerical 
calculations. In all cases at high Mach Numbers when 
the wave number is not small, the relation for the char 


acteristic values is more complicated. 


(2) GENERAL FORMULATION OF PROBLEM 


The basic differential equations for the motion of a 
compressible fluid, in Cartesian tensor notation, are as 
follows: 

Equation of motion 


(Ou; Of*) + [u(Ou; Ov;)| = 


(1 p*) (O7;;/Ox, 2.1 


Equation of continuity 


(Op*/ot*) + (0/Ox,) (p*u,) = 0 


Equation of energy 


p*C,|(O7T*,/ot*) + [uj(O7* Ox,)]) = 
(0/Ox;) [k*(OT*/Ox,)] + rier, (2.3 


Equation of state 
p* = | a a ling (2.4 


In these equations e;; is the rate-of-strain tensor 


and 7;, is the stress tensor, defined in the following way 


ei, = (1/2) [(Ou,/Ox;) + (Ou,/Ox;)] (2.0 


j 


Ti; = 2u*e;, + [(2 3) (A* — we ex; — p* \6:, (2.6 
All quantities are dimensional, as defined in the List of 
Symbols (p. 477). The coefficients of viscosity 4* and 
\*, the coefficient of heat conductivity k*, the specific 
heats C, and C,, and the Prandtl Number o are con- 
sidered to be functions of the temperature alone. The 
gravitational force terms are omitted, since they call 
be shown to be negligible for the large values of the 
Froude Number to be expected in most practical prob- 
lems (see reference 1). 

We consider a time-independent basic flow and a 
small time-dependent disturbance. Thus, any quat- 
tity Q(x;, (*) may be written as the sum of a part 
O(x;) and a part Q’(x,, f*): 

(2.7 


Q(x, *) = Oxi) + O' (xs 


We are interested in determining whether the disturb- 
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ances represented by Q’ are stable or unstable—that 
is, whether they decrease or increase with time. Both 
the quantities, Q in the disturbed motion and Q in the 
steady motion, satisfy Eqs. (2.1) to (2.4). Thus, if 
expressions of the form of Eq. (2.7) for all the variables 
are inserted into these equations and quadratic terms 
in the small disturbance quantities are neglected, the 
following system of linearized equations for the fluctu- 


tions Q' is obtained: 


Om, Ox,)p*’ + p*i(Ou,’ Of*) + [(OM;/Ox,)u,"] + 
[a ;(Ou;’/Ox;)]} = O7;;'/Ox; (2.8) 
Op*’ Ot*) + (0/On;) (4;p*’ + u,'p* Q (2.9) 
m .. aT’ oF? o7*’ 
p*’ + p* + u;, +t 
On ot* Ox Ov 


1 Oo | Coa o7*’ (cary er 
Le Ox... & Ox, . 


The quantities e;;’, 7’, (C,u* o)’, (LC 


u*’ are given by 


€:;' = (1/2) [(Ou,’/Ox,) + (Ou,'/Ox,)] (2.12) 


and 
Cu*/o)! = (d/dT*) (Cyu*/o)T* 
O/C) ed/dT) {i/ea 
1” = Ga /di” a 


g®’ <= (de*/dT*)T’ 


(2.14) 


Equations (2.8) to (2.11) form a system of linear 
partial differential equations for the fluctuations. The 
coeflicients are known functions given by the steady 


flow. 


3) STABILITY OF PARALLEL AND NEARLY PARALLEL 
FLOWS 


When the fluid is incompressible, the basic steady 
flow can be strictly parallel, as in the case of two- 
dimensional flow through a channel. The formation 
of the stability problem is then simple. For compress- 
ible fluids, however, there is no such simple basic flow. 
Flows satisfying the boundary-layer approximation are 
the closest to parallel flows. It will be shown that for 
such flows the stability problem is essentially the same 
as that for a parallel flow which agrees locally with the 
basic flow of the boundary-layer type. 

For proving the theorem in Section (4) and for appli- 
cation to boundary layers in three-dimensional flow, 
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the development will be made in a somewhat more 
general manner than usual. We consider a semi 
infinite flat plate in an oblique incoming stream. The 
free stream has components of velocity (2) and l’ 
parallel to the plate, and the flow inside the boundary 
layer is parallel to the free stream at each location 
(11, 3). Indeed, the boundary layer can be described 
in terms of strips parallel to the direction of the result- 
ant free stream. Each strip behaves like a flat plate 
with the leading edge normal to the resulting stream. 
The location of the leading edge coincides with that of 
the particular strip in question. Thus, if s is the dis- 
tance of a point from the leading edge in the direction 
of the free stream, the boundary-layer thickness 6 at 


that point is given by 


5/s = O(R OCR” 3.1) 

where 
R, = (sv U,? + U;*)/in* 3.2) 
R = R, = (6V U,? + U;3*)/n%,* (3.3) 


and 7,;* is the kinematic viscosity in the free stream. 
As is well known, the quantities Q associated with the 
basic flow satisfy the following order-of-magnitude rela- 
tions (if the free-stream Mach Number is not too 


large): 
(00 Oxz): [(0Q/Ox1) or (OQ Ox3)] = 
Ok, *) = OR) 34) 
(%, OF 3): to = O(R,’?) = O(R (3.5a) 
T* = 0(7,*), p* = O(7:*), a@* = O(f;*) (3.5b) 


We shall now simplify the equations for the small 
fluctuations by means of analogous considerations re- 
garding the fluctuations and their derivatives. From 
physical considerations and also from the mathematical 
nature of the boundary-value problem associated with 
the system of disturbance equations [Eqs. (2.8) to 
(2.11)], it may be expected that the fluctuations have 
the typical boundary-layer type behavior in a region 
adjacent to the solid boundary. That is, near to the 
boundary, viscous effects and therefore the highest 
order derivatives in the direction normal to the solid 
surface are important, and changes in the normal direc- 
tion are much more rapid than in any direction parallel 
to the surface. 

It should first be recognized that the fluctuations have 
a general length scale of the order of the boundary- 
layer thickness, so that 


1 Q) (0Q dx) O(R) (3.6) 


(1/Q’) (00’/dx;) = 


In the viscous region next to the surface, the character- 
istic length in the x2 direction (normal direction) is much 
smaller and is obtained by considerations analogous to 
the usual boundary-layer arguments for the basic flow. 
It is clear from Eq. (2.8) that the equation of motion 
for wu,’ contains a term of the form p*#,0u;'/Ox; from 
the inertial terms and a term of the form f*07u;'/Ox,? 


from the viscous terms. Thus, if these terms are to be 
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of the same order, .2-derivatives must be larger than 


x)-derivatives. Now, 
p*i,(On,’ Ox)) = O[(p1.* Uy’) 6} (3.7) 
pa* (07u,'/Oxe"”) = O[ft*(u'/52)R?”"] (3.8) 


where R” is a factor indicating the difference of order 
of magnitude between x; and x. derivatives: 


(0Q0’ Oxz): [(00"/Ox,) or (0Q’/Ox3)] = OCR") (3.9) 


A comparison of Eqs. (5.7) and (3.8) gives n = 1 2. 
The relation (3.9) is now comparable to Eq. (3.4) with 
R, replaced by R;. The relation analogous to Eq. 
(3.5a) is 

u»’ = 0(R’”’) 


(u,’ or u3’): (3.10) 


An examination of the equations also gives the follow- 
ing relations [which are used in Section (4) ]: 


| ad 7" = O( uy’ U)), Poa pi* = O(1,’ U;) (3.11) 


Therefore we have a boundary-layer phenomenon asso- 
ciated with the fluctuations within the boundary layer for 
In a thin layer next to the solid surface 
the fluctuations 


the mean flow. 
with thickness of the order of 6R 
have the viscous behavior characterized by rapid vari- 
ation in the x2-direction, while outside this layer they are 
usually independent of viscosity since the viscous 
terms in the equations are usually unimportant there. 
In general, the order-of-magnitude relations (3.9) 
and (3.10) are valid only in the viscous layer at the 
solid surface, while outside it the velocity fluctuations 
are all of the same order of magnitude and the char- 
acteristic length in all directions is the same-—that is, 
of the order of the boundary-layer thickness 6. 

In all of these order-of-magnitude estimates O(R ’*) 
may denote a quantity not actually of the order of 
magnitude of R“* in the strict sense, but smaller. 
Such a situation occurs in certain layers of the fluid 
where it is shown by a more complete analysis (see 
reference 6) that 0(R”*) is replaced by 0(R'’). In 
fact, when the disturbance is analyzed into periodic 
components [see Section (5)] and the ordinary differ- 
ential equations for the amplitudes of these components 
examined, it is found that such a layer (‘‘critical layer’’) 
always occurs at the x.-location in the boundary layer 
where the wave speed in the x,-direction of the par- 
ticular component considered is equal to the component 
of the velocity of the basic flow in that direction. There 
is the possibility that in certain cases this critical layer 
(where viscosity is also always important) may occur 
outside the viscous layer at the solid surface described 
above. In such cases, relation (3.9) with 0(R’’) re- 
placed by 0(R’”’) will also be valid in the critical layer, 
so that variations through this layer in the x2-direction 
will also be extremely rapid, just as in the layer at the 
wall. The order-of-magnitude analysis of the disturb- 
ance equations carries through in just the same way, 
so that it is not necessary for present purposes to con- 
sider the more complicated properties of the solutions 
of these equations. 
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(4) SIMPLIFIED EQUATIONS FOR SMALL DISTURBANCRs 
IN FLOWS OF THE BOUNDARY-LAYER Typ] 

By making use of the boundary-layer approximations 
outlined above [corresponding to relations (3.4), (3.5). 
(3.9), (3.10), and (3.11)] it can be shown that the equa 
tions for the small disturbances can be simplified to the 


following: 


_¢({Mm . Ou’ On,’ , Ot; 
a + UU; tT ie = 


T ut A 2 
or* On) Ox: OX» 
Op*’ O71,’ 
- + j* f.] 
On; On» 
i Ons’ _ Oy’ Ou»’\ _ Op*’ _, Our’ 
p . ity + Us =— + fi 
Ot Ox) OX: Ox» Ov,” 
1.2 
‘ (“ _ Ou;' Ou,’ On ) 
p ee + 1s + Us 
or Ox On: OX» 
weg os 
OX; One" " 
Op*’ Op*’ Op*’ Op* 
+ if + i + uy’ = —p*A’ | 
or* Ox; ON; Oe os 
ef fa ,or* 
BC. - + i + is + Ue» 
Ol On) OX: OX 
_ , iid 
— pr’ +k (4.5 
X's 
p*'/p* = (p*’/p*) + (T*'/T*) (4.6 
where 


A’ = (Ou;'/0x) + (Ome! /Ox2) + (Ou3'/Ox3) (4.7 
Equation (4.5) can also be written in the form 


P os we. _« a) 
p ty + u + Us = 


p = 113 
ot* Ox) Ox; Ox» 
1 /op*’ Op*’ e) rer 
: + u + t3 + (4.8 
C; ( ot* Ox; Ox; a Ox” 
where o = C,j* k*. These equations are correct ex- 


cept for terms O0(R™ “*). 

A few typical steps in the process of simplification 
are discussed below: 

Example 1.—-In the inertial terms of the equations 
of motion we have terms of the kind 74,(Om,'/Ox;) + 
2(Ou,’/Ox2). Only the first term is included in Eq. 
(4.1). This is justified by the fact that 


\(Ou,’/Ox,) = 
(te M1) [(Oz,’ Ox) (Ou,’ On)) | (4.9 


ho(OU,’/Ox2): 


which is 0(R~) outside the viscous layer and 0(R 
inside the viscous layer, according to Eqs. (3.4), (3.9), 
(3.9), and (3.10). 

Example 2.—None of the terms arising from the dissi- 
pation function (re + p*du;,/dx,) are included in 
iq. (4.5). This is justified as follows. The typical 
terms of the largest order which could contribute to Eq. 
(4.5) are 
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u*’(Oi,/Ox2)” + 2f*(OM)/Ox2) (Ou;'/Ox2) + 
Their orders of magnitude can be estimated by com- 
paring them with @*(O#,/Ox2)*, which is the largest 
term in the dissipation function for the mean flow. 


Thus, 
w*"(Oit,/Ox2)? = (u*"/M*)A*(OM)/Ox2)? — (4.10) 
7*(Oim, ON2) (Ow)’/ Ox2) = 
2[(Ouy’/ Ox») (Oily Owe) ]A* (Oi, xz)? (4.11 


On the other hand, the term p*7,C,07*’ Ox, retained 
on the left-hand side of Eq. (4.5) can be estimated as 
follows 

5* it, C OT *’/Ox) = 


* iC, [(OT*’ Ox) /(OT*/Ox,)] (OT*/Ox,) (4.12) 


Now, from the energy equation for the basic flow, 
5*ii,C,(O7*, Ox,) is of the same order as the dissipation 
term f*(O7,;/Ox2)*. Consequently, the various quan- 
tities given by Eqs. (4.10)—(4.12) are in the ratios 

|: R outside the viscous layer and in the ratios 1:R 

R within the viscous layer. We are therefore justified 
in neglecting the first two quantities, which are the 
largest contribution of the dissipation function to Eq. 
1.5). 
In addition to the approximations just described, we 
can also neglect the x; and x; dependence of the quan- 
tities in the basic flow, which determine the coefficients 
in Eqs. (4.1) to (4.6). This approximation, which has 
been repeatedly discussed by various writers on the 
stability problem, depends on the boundary-layer na- 
ture of the mean flow and can be shown to be equiva- 
lent to neglecting terms of the order of R~! (see refer- 
ences 6, 8S, and 10). 

The system of Eqs. (4.1) to (4.6) is essentially of the 
form expected from a strictly parallel flow. One im- 
portant characteristic of the system is the absence of 
dissipative terms in the energy equation. This makes 
it possible, as we shall see in Section (5), to reduce the 
problem for a three-dimensional disturbance to a two- 
dimensional problem. In previous investigations such 
a reduction was not possible because of certain terms 
retained in the equations for the disturbances that 
prevent a proper transformation. Such terms are now 
shown to be negligible. 

The discussion up to now has been restricted to the 
case of moderate Mach Numbers, for which relations 
o1), (3.4), and (3.5) are valid. For very high 
supersonic Mach Numbers, the order-of-magnitude re- 
lations in the boundary layer corresponding to relations 
3.1), (5.4), and (3.5) are different (see reference 5), 
and, consequently, those for the disturbances corre- 
sponding to Eqs. (3.9), (3.10), and (3.11) are different. 
An order-of-magnitude analysis (see reference 6) similar 
to that described for the case of low Mach Numbers 
shows that the disturbance equations are again, in the 
first approximation, of just the same form as Eqs. (4.1) to 

1.6). The proper velocity and temperature distribu- 
tions given by a modified boundary-layer theory (see 


reference 4+) would have to be used, of course, at 
very high Mach Numbers. The error estimate for 
the neglected terms in the disturbance equations is 
now different, however, being expressed not in terms of 


R but in terms of a Reynolds Number 
Ry = (v:*/%0*)R 113 


based on fluid properties at some representative tem 


perature 7)* inside the boundary layer. Since 7)* is 
roughly proportional to .1/;° at high Mach Numbers, 
Ry, is much less than R in such cases. Consequently, 
the errors committed in the simplified equations [which 
are O(R, “*)| might not be very small at high Mach 
Numbers. 

The present analysis is carried out for any boundary 
layer over a flat surface, whether there is a pressure 
gradient or not. Thus, the results are valid, not only 
for the zero pressure-gradient flow over a flat plate, but 
also for the flows over wedge-shaped profiles. Although 
a pressure gradient is associated with the latter flows 
at subsonic Mach Numbers, terms due to pressure gra 
dient do not appear explicitly in the disturbance equa- 
tions to the accepted order of approximation. T 

Although the simplified equations | Eqs. (4.1) to (4.6) | 
are obtained by rather crude order-of-magnitude argu 
ments, an a posteriori check verifies that the approxi 
mations made are valid in the stability investigation of 


the boundary layer over a solid surface. 


(5) PERIODIC SOLUTIONS OF THE EQUATIONS FOR SMALL 
DISTURBANCES 


The Eqs. (4.1) to (4.7) contain x), x;, and ¢* only 
through the partial derivatives, so that solutions of the 


form 

Q’ = g(xz) exp [i(aw, + ayx, — ayc*t*)] (5.1 
are possible. It may be expected that a disturbance 
of a fairly general nature can be analyzed into normal 
+ 


modes of this type, with a; and a; real and c* in general 


complex. A normal mode is damped if c;* < 0, neutral 
if c;* = O, and self-excited if ¢;* > 0. 

Evidently solutions of the form (5.1) may be regarded 
as waves propagating in a direction in the (1), x;) plane 
specified by the numbers (a, a;) (with changing ampli- 
tude in time if c* is complex). 

Consider in particular the case a; = 0. 
terms involving the x3-derivatives drop out of the Eqs. 
(4.1) to (4.7), and the equations for ", 12’, p*’, p*’, T*’ 
are exactly the same as if the component #; of the basic 


Then, the 


+ A more careful examination is necessary when there is sur 
face curvature. Preliminary results indicate that for most 
three-dimensional disturbances curvature effects are negligible 
when the curvature is small, but that for those disturbances with 
propagation direction nearly normal to the free-stream direction 
see Section (5)] additional terms have to be retained in Eq 
(4.2) 

t Free boundary layers (such as wakes and jets) are excluded 
from all of the above discussion. The effect of viscosity on the 
stability of such flows is quite different. Further investigation 
of these cases is being carried out 
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flow were zero. Thus, for waves propagating in the 


x;-direction only the component #, of the basic flow 
(besides the distributions of the scalar quantities p* and 
Now, a rotation of the co- 


7*) has any influence. 
ordinate system in the (1), x3) plane will not change 
the system of Eqs. (4.1) to (4.7). Hence, we may con- 
clude that for waves propagating in any direction only 
the component of the basic flow in that direction has 
to be considered, in addition to the distributions of the 
scalar quantities. very periodic three-dimensional 
disturbance may be treated in terms of a corresponding 
two-dimenstonal problem. One need only take the 
proper component of velocity in the basic flow to carry 
out the simplification. 

If we reduce the above equations to dimensionless 
forms, it 1s clear that the waves propagating in the 
direction of the resultant free-stream velocity are 
associated with the Reynolds Number R and the Mach 
Number .J/, of the basic flow. For a wave propagating 
at an angle 0 relative to this direction, these numbers 


are reduced to 
R = Reos0, M, = M, cos 0 (5:2) 


A direct transformation verifying this statement will 
be given in the next section. 

The advantage of the present treatment lies in its 
generality and simplicity. Whether the incoming 
stream is oblique or normal to the leading edge of the 
flat plate makes no difference in the discussion. When 
specialized to the incompressible case, the present re- 
sults cover both the original work of Squire* and the 
case of an oblique incoming flow discussed by Kuethe.’ 
For two-dimensional parallel flows, Squire obtained a 
set of transformation formulas leading to a conclusion 
equivalent to Eq. (5.2). Supplementary discussions 
are required in Squire's method. Kuethe, dealing with 
the case of the oblique incoming stream by Squire's 
inethod, noted the variation of ‘‘effective’’ Reynolds 
Number with different directions. He did not, how- 
ever, arrive at the conclusion reached here that the 
maximum “‘effective’’ Reynolds Number is reached by 
waves propagating in the direction of the free stream. 
Consequently, these waves first give rise to instability as 
the Reynolds Number of the flow is increased. 


(6) EQUATIONS FOR THE AMPLITUDE FUNCTIONS 


Consider now the stability of a two-dimensional 
boundary layer over a flat plate coinciding with the 
(x1, x3) plane and introduce the dimensionless vari- 
ables given in the List of Symbols (p. 477). The free- 
stream velocity is now Ll), since U,; = 0. The free- 
stream reference quantities 6, 14, pi*, T,*, p:*, w1* are 
evaluated at a fixed point 1; = x,"—that is, at a definite 
distance from the leading edge of the plate——in the part 
of the boundary layer that is of interest. 

The dimensionless velocity and temperature dis- 
tributions, w = a,/U, and 7 = 7T*/T,*, depend on a 
number of factors—in particular, the boundary condi- 
tions at the wall and the free-stream Mach Number ./;. 
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We consider a small disturbance of the boundary layer 
in the manner discussed previously, with solutions of the 
form 


QO'(x, y, 2, 1) = gly) exp [t(ax + Bz — act)] (6.1 


for the disturbance quantities, as indicated in the List 
of Symbols. When relations of the type [Eq. (6.1 

are inserted into Eqs. (4.1) to (4.7), the following equa- 
tions for the amplitude functions f, ¢, h, 7, 7, 6 are 


obtained: 


plia(w c)f T w’ad| = — |[(1a) (yAT,7) |x 


(u/R)f" (62 
plia’(w — co] = — [1 (yM,? |x’ + (un Riad” 6.3 
plia(w — c)h] —[(78) (yM7)|a + (u/ RK” (64 


la(w — co)jr + p'ad t p(taf + 78h + ag’) 0 6.5 


plta(w — c)@+ T’ad] = 


=“ = 1) (taf Tt ibh . ag’) ." [ (yu) (oR a” 
(6.6 

Tv r @ 
=-+-— (6.0 

Pp p 1 


where a prime now refers to differentiation with respect 
toy. 

The combination of terms af + 8h suggests its re 
Indeed, the follow- 


placement by a new variable a f. 
ing transformations convert Eqs. (6.2) to (6.7) to those 
of an equivalent two-dimensional problem: 


af+ ph = a 

ap = ap 

an = an 

a? = ar 

ad = a6 (6.8 
«2 


Va't+tBP =a 

aR = aR 

al, = aM, 
The two-dimensional equations are: 


plia(w aa c)f + wad | = — [(7a) (yM,2) |z tT 
(u Rk)" (6.9 


~/ l 


plta?(w = cjg] = — [1 (yM,”) Ji 7 te R) ad" (6.10 
ta(w — c)F + p’ad + p(iaf + ad’) = 0 (6.11) 
plta(w — c)8 + 7’ ad] = 


—(y — 1) (af + ad’) + [(yu) (oR) 0" (6.12 
t/p = (F/p) + (6/T) (6.13 


Evidently the transformations [Eqs. (6.8)] are 
equivalent to Eq. (5.2) with cos 9 = a V a® + 8°, 90 
that 6 does in fact give the direction of propagation of 
the waves, Eq. (6.1). The dimensionless wave speed 
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cis not changed, which indicates that the phase velocity 
in the x-direction is actually sec 0 times the phase ve- 
locity in the direction of the wave normal. 

To complete the analysis, it is necessary to verify 
the equivalence of the boundary conditions. In the 
stability problem the boundary conditions are homo- 
geneous, so that there is no difficulty in this respect. 
The homogeneous boundary conditions of the three 
dimensional problem clearly transform to the proper 
homogeneous boundary conditions of the equivalent 
two-dimensional problem. For forced oscillations, how- 
ever, it is obvious that the equivalence would have to 
be more carefully examined. 

We have now shown that the boundary value problem 
associated with the system of Eqs. (6.2) to (6.7) for a 
three-dimensional disturbance is equivalent to that 
associated with the system of two-dimensional Eqs. 
6.9) to (6.13). Thus the stability problem for three- 
dimensional disturbances can always be formulated 
as a two-dimensional problem which is mathematically 
the same as that for two-dimensional disturbances. It 
should be realized, however, that the equivalence is 
purely mathematical since Eqs. (6.9) to (6.13) do not 
represent a proper two-dimensional disturbance. There 
is a decisive difference between the present results and 
those of Squire for an incompressible fluid. In the 
latter case, because the Mach Number does not enter, 
it is possible to transform the equations for a three- 
dimensional disturbance to equations that actually 
do represent a real two-dimensional disturbance at a 
lower Reynolds Number. It is then possible to conclude 
that three-dimensional disturbances can be ignored 
when the minimum critical Reynolds Number is being 
investigated. For a compressible fluid this is not true 
because of the necessity of transforming the free-stream 
Mach Number. It is not even true that every three- 
dimensional disturbance is equivalent to a two-dimen- 
sional one with a modified Mach Number as well as 
Reynolds Number, because with a reduced Mach Num- 
ber the velocity and temperature distributions would 
In the equivalence dis- 


also have to be modified. 


cussed above these distributions must remain the 
same. 

To apply the results of a two-dimensional analysis 
to the three-dimensional problem we may use the trans- 
formations (6.8), with the quantities g associated with 
the two-dimensional problem. It is convenient, there- 
think of three-dimensional disturbances as 
waves propagating in various directions. For each 
value of M,—that is, for each value of the direction 
angle of the waves-the stability problem is investi- 
gated by means of Eqs. (6.9) to (6.13) with the usual 
Thus, for 


fc re, to 


methods of the two-dimensional theory. 
each value of M, there is a neutral curve in the (a, R) 
plane which transforms to one in the (a, R) plane by 
transformation (6.8). The family of 
all such neutral curves, for all possible values of Mj in 
the range 0 << ,< M, gives the stability character- 
istics of three-dimensional disturbances. 


means of the 
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(7) FORMULATION OF THE CHARACTERISTIC-VALUE 
PROBLEM 


Having established a two-dimensional formulation of 
the problem for three-dimensional disturbances, we 
may make use of the theory of Lees and Lin! for further 
investigations. However, some modifications are nec- 
essary when the wave speed of the disturbance is not 
small compared with the free-stream speed. For super- 
sonic free-stream velocities this is usually the most 
important case. Detailed 
of the necessary modifications will be given in a sepa- 
rate paper. Here we shall merely give a summary of the 


mathematical discussions 


theoretical arguments and indicate how the final for- 
mulas differ from the existing ones. For a brief outline 
of the mathematical theory see the Appendix. 

To determine the analytical nature of the system of 
Eqs. (6.9) to (6.13), it is convenient to introduce the 


variables 


Zi, =] 
Zs = f° 
Z = We - 
Zi = =/Mi | ~ 
Z,=06 
en ¥ 


Then the equations transform to a system of six equa 
tions of the first order. The latter, when reduced to 
normal form, have coefficients analytic in the inde 
pendent variable y and the parameters 1°, @R, 4’, 
and c. 

Six homogeneous boundary conditions must be satis 


fied for a natural oscillation. These are given by 


f,@ 6 bounded as y > « | 
f (0) = ¢(0) =0 (7.2 
ao’(0) + 56(0) = O | 


where a and 6 are constants depending upon a, c, 0 and 
the physical properties of the gas and the adjacent solid 
wall (see Appendix). 

The condition of boundedness for y — © can be put 
into an analytical form with the aid of the differential 
equations [Eqs. (6.9) to (6.13) ]. 
effects become negligible, so that the terms of order 
1 R can be neglected, and the equations reduce to 


For large y, viscosity 





(d/dy) \|(w — c)d’ — w'o)/[T — Mi2X(w — c)?]} - 


{[~9 Tl 2 79 
, [ar(w — « (9 =O (6.5 
For y — © this is easily seen to become 


Ww" 


o” — Bo = 0 / 


~J 
res 


8? = [1 — M21 — )?)\ 


Equation (7.4) has two solutions exp(+y), where 
(for real c) 6 denotes the positive square root 


B= aV1 — M1 — c)? (7.5 
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Thus, for real c the condition of boundedness at in- 
finity is automatically satisfied if 


c<1— (1/M) (7.6) 
In this case, the characteristic values are continuous, 
and the solutions represent disturbances whose phase 
velocity in the (x, z) plane is supersonic relative to the 
component of the free-stream velocity in the direction 
of propagation of the disturbances in this plane. These 
supersonic disturbances correspond to sound waves 
propagating into or out of the boundary layer (see refer- 
ence 1). 
For subsonic disturbances with 
c>1— (1/M) (7.7) 
the solution of Eq. (7.3) which is bounded at infinity 
is exp(—@y). Therefore the boundary condition for 
yo © is 


o'¢—> —Basy—> ~ (7.8) 


In this case the characteristic values are discrete and 
are determined from a secular equation which can be 
reduced to the form 


E(a&, c, M,?) = F(z) }1+ A} (7.9) 
In this equation, F(z) is the function of Tietjens defined 
by 
o—¢g » 
F(z) = 1+ | eee, SY (er) ** tae 
& i 
. Zz 2 ; \ 
3 | OH, | — (ig) | dee (7.10) 
; 3 
where 


yp 


5 3 °c c¢— @ . = 
S = tan)’ < \ dy (4.22) 


and the function /:(a, c, M;") is given by 


C= & 5 "0 [-— 2 ] E ; vi _ 
ay “(@, €, di . — 
v 2 \ v , 7 


dix Ou’ + Bou 


Pox doy’ + Bd 


(Tie r Mw, ‘COiw) 8 = M,°*c*) | ou’ 7 Bi 
[(Tdew’ + Mi? Ww'Ch2w) (Te — Mi?c?)] bx’ + Boer 
(7.12) 


The quantities ¢),, ¢2. are the wall values and 1, ¢21 
the free-stream values of two linearly independent 
solutions of Eq. (7.3). 

The function A = A(a, c, M2, a, 6, 2, 7), whose ana- 
lytical definition is given in the Appendix, depends on 
both the inviscid and viscous solutions of the dis- 
turbance equations and on the thermal boundary con- 
dition in Eq. (7.2). 

The principal difference between the present formu- 
lation and that given by Lees and Lin!’ is in the presence 
of the function A. The appearance of this new term 
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in the secular equation is due essentially to the fact that 
we are now placing no restriction on the value of , 
In fact, for finite values of 7°,,7 A can always be written 


as 


where A; is another function depending on the same 
Thus, it is clear 
On the 
other hand, at high Mach Numbers, when c is usually 


quantities as A which is 0(1) at most. 
that A is always negligible when c is small. 


large, A is in general 0(1) and cannot be neglected, since 
the value of A, depends on the thermal boundary con 
dition and is therefore small only in very special cases 

In reference 1 it was concluded that the character- 
istic value problem is independent of the temperature 
fluctuations and the boundary conditions imposed on 
them. The 
conclusion is of general validity only when c is small, 


present considerations show that. this 
since it is only then that A can be neglected without 
any restriction being placed on the value of A). For 
practical purposes, this limits the range of validity of 
the conclusion to subsonic and _ slightly supersonic 
Mach Numbers, when c is not large in most cases of 
interest. For example, the error corresponding to the 
neglect of A for an insulated wall at J/, 1.6 can be 
estimated from Eq. (7.13) to be about 10 per cent (on 
My). 
considering the usual order of accuracy of the asymp 
totic theory. 
Numbers the error increases rapidly. 


taking c ~ 1 — 1 This is still not very large, 


However, for somewhat higher Mach 
The preceding conclusion is true in an extremely re 


Mach 
those disturbances whose direction of propagation 1s 


stricted sense at high Numbers—namely, for 
at a large angle with the free-stream direction so that 
Mis small. For these, it is clear that A is again negli- 
gible without any restriction being placed on the value 
of A). 

In the cases just discussed the secular equation re- 


duces to 
E(a, c, MW?) = F(z) (7.14 


which has a form much like that in the theory of refer 


ence 1. In all other cases, the form of this relation 
depends on the thermal boundary condition in Eq. 
(7.2) 


A detailed analysis of the thermal boundary condition 
leads to the important conclusion that when the wave 
number a& 1s very small the relation for the character 
istic values again reduces to Eq. (7.14), since A, then be- 
Such a 
fluctuations, 
An important example 


corresponds to 
time- 


comes negligible. situation 


very low frequency since the 
frequency is given by ac cos 0. 
is the limiting case of complete stabilization of the 
boundary layer by cooling, for which all values of @ on 
the neutral curve must approach zero in the limit [see 
Section (10)]. Thus, the critical temperature ratios 
for complete stabilization can still be obtained by fairly 
simple calculations with the use of Eq. (7.14). 


+ The exceptional case 7, = 0 is discussed in the Appendix 
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In all cases at high Mach Numbers when the wave 
number is not small, the relation for the characteristic 
We shall 


not consider these more complicated situations any 


yalue is more complicated than Eq. (7.14). 


further here and shall limit all subsequent calculations 
in this paper to situations in which Eq. (7.14) is a valid 
approximation. The detailed presentation of the ex- 
tended theory will be given in future communications. 
\n examination of Eq. (7.14) shows that it reduces to 
the corresponding relation given by Lees and Lin for 
small c but that it is somewhat different for large « 
because of the changes in the definition Eq. (7.11) for 
zand the first factor in Eq. (7.12). The results of the 
previous theory are reproduced by assuming c to be 
small and keeping only the first term in each of the 
power series expansions of (w — c) and vy about y,. In 
the supersonic case, when ¢ is not small, such a pro- 
cedure does not represent a good approximation for our 
purposes.t For a discussion of the differences between 
results obtained by the present and older theories, 
see Section (10). 
8) DETERMINATION OF THE CURVES OF NEUTRAL 
STABILITY FOR OBLIQUE WAVES 


Consider oblique waves with direction angle are 
cos M, 14, where M, is some value in the range 0 < 
M,< ™M,. 
values (a, R, c) can be transformed to 


The condition (7.14) for the characteristic 
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u+iv = 1+ (w,'c/T,) ) [ba’ + 


aV1 — M1 — c)*u]/ [on’ 4 
aV1l — MA(1 — c)? gn}} (8.2c 
and $11, ou, ¢1', gx’ are as defined by Eqs. (S8.6a 
(S.6d). 
For some purposes it is more convenient to solve Eq 


(8.1) for u + iv, so that 


u+w= WV,(z,c) + Viz, « (8.3 
where 
V,(z, (1 + A)®, — ACO? + O,*)] 
{1 + AC — ®,)]? + A2@,7} (Sida 
V(z,c) = [(1 + A)®,]/, [1 + A — &,)]? + APb?} 


(S.4b 


Before discussing the determination of the neutral 
curve we shall write down the required formulas and 
make certain transformations that are useful in the 
calculations. The quantities ¢;, and @»; are the free 
stream values of two linearly independent solutions of 
the inviscid equation 
(d/dy) \|(w — c)d’ — w’g] [T M,*(w — c)?}} 

&i[(w—c) The =0 (85 


These two solutions may be obtained as power series in 


a’, as described in reference 1, so that 


(zs) = (1+ A) | (w+ w)/[1 + A(u + w)]f (8.1) 
- ; édu = (1 — oc) > a*"FAb, (S.6a 
with b(z) = 1/[1 — F(z)], where z, \, and u + iv are n=0 
given by 
gn = (1 -— oc) > # a*" Ko, (8.6b 
. (3 l" c—w ) , seals 
~ > 3 ‘ 
aR = 3° dy (8.2a) 
oJ, \ v x , 1 — M*(1 — ©)*] . 
Pu = (1 — © > a » Hp 1 
P ' (1 — c) . Pee 
3 Wy Vy ve lc — Ww , 
A(c) = 7, \ dy — 1 (8.2b) (S.6¢ 
= m v : 
1— M21 — c)? 
+ The approximation is worse in » than in (w — ¢ én = (1 — ¢) (l —<¢ 2X ill 
— Hn 0 
(S.6d 
where 
~ "1 (w — c)° ~ - r ine M7 (w —c)-* ‘iw «€ 
m= 1, Ay = } — dy, H.= / dy [ dy, 
Ji) 7 e io == ig" a t 7 
Pp (S.ia 
~ "l(w— c)° y 7 = VW, ( — c)- Y(w— c)- 
a | dy | dy [ dy 
a 7 J | iw- = J0 7 
> ~ m1 i = M?(w — c)° ~ “lL (% ¢} ) a a Mw — ¢)° 
Ko=1, Ky = | dy, Kez= | dy [ dy, 
Jo (w— Cc) Jt Jo (w— ¢ : 
(S.7b 


For values of c that are not too large, the integrals 
‘./a) and (8.7b) can be put into forms that are more 


suitable for calculation purposes. For example, in 


Ko we use the relation 
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~7T — M,*(w — c)? 1T — M,2(w — c)? 
} = dy = z dy — 


(w— c)* (w — c)? 
'T — M,2(w — c)? 
dy (8.8) 
ae w— Cc)? 
to transform A» into 
Ry = AR, = Ne (8.9) 


where 


~, liw — e)* 17 — M,2(w —c)’ 
No = { es dy | - dy 
0 ? ie (w — c)? ° 


Thus .\2 has an integrand 


(w—c)? ('T — Mw — c)? 
| dy (8.10) 


ti , (w — c)? 

which is not numerically large for any value of y, so 
that N;2 may be expected to be numerically smaller 
than A». Similar transformations are made for the 


other integrals as follows: 


K,, = R.A, ;=—= N,, n = o 


’ 
‘ 


No = (8.11) 


H, = AH, -» —_ M,, n > ES 


In general, M,,, and N,, will be numerically smaller 
than the corresponding H,, and K,, for small c. The 
preceding transformations were first used in the 
stability theory for incompressible fluids,’’ where c is 
usually fairly small, and then extended to the case of 
compressible fluids. Obviously they are not useful 
for reasonably large values of c, in which case other 
transformations may be more appropriate. This situ- 
ation occurs for large enough Mach Numbers—for 
M, greater than about two, perhaps. When J/, be- 
comes very large, in fact, c is often very close to unity. 

The relation (8.2c) may now be written as 

ee Pi 
u+Ww= = IR, + [T7,,/(wy’c)] + (1/a)G(a, c)' 
T,, | f 
(8.12) 
where 
- | * 1-— M,(1 — c)? 
x 


G(a@, c) = (1 — Ara?) - 
vee aie (1 — c)? 


(aa — > a Tess) | (8.13) 


l 
If G(a, c) is expanded as a power series in a, the real 
and imaginary parts of Eq. (8.12) are given by the 
following expressions in the case of real c: 
u = (Wy't/Tw) \ Bi, + [Te/(wy'c)] + (k/a&) [1 - 
kPa + (kA? — 2Ae,)a? — +]} (8.14) 


v= (Wy 'C Es) (Ay; + (k a) (—2H, a+ — 


where 
k= [1 — MP1 — c)?7] 7/01 — 


Such series expressions are expected to be good ap- 
proximations when ka is small. This is true for al] 
cases at subsonic and slightly supersonic Mach Num 
bers, but at high supersonic Mach Numbers it appears 
to hold only for cases of extreme cooling since otherwise 
k may be very large. 

In the present work we shall neglect all terms of 
order a in Eqs. (8.14) and (8.15) for the sake of sim- 
plicity, so that 
u = (wy'c/Ty,) \RKi, + [Te/(uy’c)] + 

(ka) (1 — kha ' S.l6a 
v = (c) = (w,'c/T,) Ky (S.16b 

In the calculations of reference 2 the same approxi- 
mation [Eq. (8.16b)] for v is used, but terms of order 
a* are retained in the relation (8.14) for u. The results 
show, however, that often the error corresponding to 
the neglect of terms of order a” in Eq. (8.16a) is no larger 
than that in Eq. (8.16b). 
mations (8.16a) and (8.16b) appear to be reasonably 


In any event, the approxi- 


adequate for all subsonic and slightly supersonic Mach 
Numbers, with errors within the limits of error of 
the asymptotic theory. Moreover, they are very 
accurate at any Mach Number when @& is very small 
(which occurs, for example, when cooling at the 
solid surface is sufficient). [See Section (10)]. For 
Mach Numbers greater than about two, however, this 
is true only for cases of extreme cooling. The present 
approximations, therefore, are probably not accurate 
at such high Mach Numbers in most cases of practical 
importance, and the method of numerical calculation 
may have to be revised. t 

Now let us consider the procedure for finding the 
characteristic values (a, R, c) on the neutral curve. 
The various steps can be illustrated most clearly by 
geometrical considerations, with the aid of the condi- 
tion for the characteristic values in the form of Eq. 
(8.5) rather than Eq. (8.1). 

For the values of c that occur in the calculations, 
A(c) is always small, so that W,(z, c) and V,,(z, c) are ap- 
proximated fairly closely by ®,(z) and ®;,(z). In fact, 
forc = 0, ¥,(z, 0) = ®,(z) and V,(z, 0) = ©,(z). There- 
fore, for each value of c, the graphs of V, and VW; against 
z have much the same shapes as those of ®, and ®;, re- 
spectively (see Fig. 1). The maximum Y,,,(c) on the 
graph of W;(z, c) against z is generally close to the maxt- 


) 


mum value ®; = 0.580 at s = 3.22 on the ®; curve. 


7 All existing calculations suffer from limitations of this kind 

; : : ' ; ; 
In some investigations, numerical calculations were carried out 
even though 


by using methods adequate only for small ¢ and a, 
Such 


the resulting values of these quantities were not small 
results must be treated with reserve 
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FiG. 1 rhe Functions ®,{2) and #,(2 


For Mach Numbers that are not greater than about 


two, the graph of v(c) against c has one of the typical 
forms shown in Fig. 2a (which is plotted from the data 
given in reference 2). In general, the greater the 
cooling that is, the lower the wall temperature— the 
higher the % curve is in the (c, 7%) plane. For larger 
Mach Numbers the v curves may have more com- 
plicated shapes, as indicated in Fig. 2b. For the pres- 
ent, we shall consider only Mach Numbers for which 
the curves have the simple behavior shown in Fig. 2a. 
The pairs of values (a, R) are to be found for a series 
of real values of c. The minimum value of ¢ for sub- 
some Mach Numbers is ¢ = 0, and for supersonic Mach 
Numbers it is either c = Oorc = 1— 1 Mif > 1. 
For each value of c, relation (8.3) is solved for z and 
u with the aid of Eq. (S.16b). First, 2 is obtained by 
solving w(c) = V,(s, c), which is equivalent to finding 
the value of s corresponding to the point of intersection 
of the curve W,(s, c) against s with a straight line 
parallel to the z-axis at a distance 7,(c) from it. From 
the behavior of the graph of W,(z, c) [which is similar 
to the graph of @,(z) in Fig. 1], we see that if m(c) < 0 
there is one intersection point (one value of z), and if 
w(c) > O there are generally two intersection points 
two values of s) if m(c) is not too large. For some 
value of c large enough it is clear that the straight line 
will be just tangent to the corresponding graph of 
¥(3, ¢) at its maximum point, and for larger values of 
¢ there will be no intersection points at all. This 


maximum value of c, therefore, is given by the inter 


section of the graph of W;,,,(c) with that of w(c) (see 
Figs. 2). Corresponding to it there is just one value of 
zs. With the value of z determined, u = W,(z, c) is ob 
tained from the graph of W,(s, « 
Then @ follows by soly 


[which is similar to 
the graph of ®,(z) in Fig. 1]. 
ing Eq. (S.16a), since all the quantities in this equation 
and finally Kk 


The corresponding values 


besides @ are now known functions of ¢, 
is obtained from Eq. (S.2a 
(a, R) are then determined by means of the transforma 
tions (6.8). The points (a, RX) for all possible values of 
c give the curve of neutral stability in the (a, R) plane. 

Evidently 2(c) and u(c) could be obtained graphically 
as just described by drawing graphs of V,(s, c) and W;,(z, 
c) against z for each value of c. However, for the pur 
pose of the numerical calculations it is more convenient 
and more accurate to solve the two real relations ob- 
tained from Eq. (8.1) for z and w by an iteration pro 
cedure with the aid of the tabulated values of ®,(z) and 
®,(z) (see Table 1), as described in reference 2. The 
maximum value of c and the corresponding values of 
z and uw can also be found accurately by an iteration 
method (see reference 13).7 

Present calculations are carried out for the boundary 
layer on a flat plate at a Mach Number of J/; = 1.6 


+ A detailed discussion of the procedure for the numerical cal 
culations described in this and following sections and complete 
tables of various functions used in the calculations are given in 


reference 13 











TABLE 1 
The Functions ,(=), ® ,®? ,®? 
' ' 

z £ (z) ¥, (2) t z) [, 2 
1 fy 2c rie =-2 557 
| 
2 l 77012 om 2252 2 71 zZ 72 
ed 2.54 
Lod 2. 2683 -1. 71669 1.505 2.937 
i.s 850 2.63 
1.¢ 2.44985 -1, 138600 4245 2.384 
rT 139C 2.127 
1.8 2.48104 - 75892 - 06428 1.907 
1.9 - .2140 1.726 
2.0 2.43927 - .41253 - .3337 1.573 
23 - .4377 1.443 
2.2 2. 3519¢ - ,1234 - 5333 329 
2.3 7 242 1.213 
ye 2.22724 + 4498 - .7108 1.1007 
2.5 - 7912 9833 
at 2.06929 31558 - .3625 Sod 
b 7 209 7254 
2.8 1, 8856¢ - 46043 - .9627 - 5853 
2.9 - 985 4414 
» 1.48932 54272 - ,985 2972 
pe - .964 159 
3.2 1,.4972¢ 58082 - .9223 2 
3.3 - .8629 
3.4 1,3251¢ 56401 - .7893 - 
1 Pe - .7054 - 
3.¢ 1.12429 51074 - 151 - 

7 - .5222 - 
3.8 1.07982 435¢€ - .4302 - 
309 - .3422 - 
4. 1.01118 §220 - ,2609 - 
4.1 - .1865 - 
4.2 97361 27133 - .1172 - 
4.3 - .0625 - 
4.4 96056 20038 - .01470 - 
4.5 - ,0061 
4.€ 95989 13601 = ,0034 
Ny + ,OL7€E 
4.8 976-59 9503 + 179% 
4.7 
5.0 99582 72 
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The approximation is now made that the minimum wher 
| critical Reynolds Number is very close to the valu 
{ of R corresponding to the maximum value of c. This ae 
i was found to be true in previous investigations for an = 
' incompressible fluid’? and for a compressible fluid 
G The neutral curves in Fig. 5 show that in the present 
case it is also an excellent approximation, with only 
a slight error. The minimum critical Reynolds Num [Ther 
ber R,, approximated in this way is plotted against the the 
; direction angle of the oblique waves in Fig. 4 and tabu- 8.24 
lated in Table +. As might have been expected, R, 
becomes infinite at direction angles of O° and 90°, R 
The minimum value of R,, occurs at J, = 1.030 (cor 
sponding to a direction angle of 49.9°) and is given 
by R,, = 1.64 X 10° This value R,. = 1.64 X 108 
therefore, is the minimum critical Reynolds Number 
for all three-dimensional disturbances in this particular an 
case. vaiue 
boun 
(9) THE DEPENDENCE OF THE MINIMUM CRITICAL w(y 
REYNOLDS NUMBER ON THE DIRECTION O} and 
PROPAGATION OF THE WAVES mine 
For those cases in which the approximations (S.16a F BN 
and (8.16b) are sufficiently accurate, it is possible to “I 
derive fairly simple analytical results for the directional 
dependence of the minimum critical Reynolds Number 
that is, dependence on the parameter Wj. By solv- i 
ing Eq. (S.16a) for @ we have “ee 


a= a (1+ kya ) 9.1 





Fic. 2a. The function w(c) versus ¢ for various ratios 7/7; 
of wall to free-stream temperature at a free-stream Mach Number Puen 
of M, = 0.7 w bs 


and ratio of wall to free-stream temperature of 7), 7) = 
1.073. This temperature ratio is the critical ratio be- 
low which the boundary layer is completely stable 
with respect to two-dimensional disturbances [see Sec- 
tion (10)]. The Prandtl Number is given the value 
0.75, and, for simplicity of calculation, the viscosity 
coefficient is taken to be directly proportional to the 
absolute temperature, with the constant of proportion- " 
ality chosen so as to give the correct result at the wall. 





This procedure for the viscosity coefficient is known 
to be a good approximation in the case of the boundary 
layer itself when the Mach Number is not too large. 
Although it probably is not as accurate in the stability 
calculations, estimates show that it does not signifi- 
cantly alter the general trend of the results. 

In Table 3, a, a,, as, R, R,, R, are tabulated as func- 
tions of c for several values of 17, in the range 0 < I, < 
1.6, where a, and R, are wave number and Reynolds 
Number based on a reference length x, (the distance 
from the leading edge of the plate) and a, and R, are 
those based on the momentum thickness as reference - 
length. oe se? ~¥ ; " 

In Fig. 3 neutral curves in the (a,, R,) plane are shown : . . 
for the three cases J, = 1-4, 1.0, 0.6 corresponding to 1yC. The fumetion o4¢) vere cfr varios Numb 
direction angles of 29.0°, 51.3°, 68.0°, respectively. of I, = 4 ‘ 
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[w/c V1 — M1 — =) 

i aad : )| (Wy ‘t Vs R 

u ———(Ki, + — = : (9.2) 
L le Wy C} | u— L 


Then the minimum value of R on the neutral curve in 
the (a, R) plane is given, with the aid of Eqs. (6.8) and 


Mm, MM ./3 7" /e—w ; 
R. = V7, R. = m,” \2 / \ : dy Xx 
1 + kha 


a 


where ¢ is approximated by the maximum permissible 
value of c described previously. Once the particular 
boundary-layer flow is specified—i.e., the functions 
w(y) and 7(y)—-this value of c and therefore 2(c), u(c), 
and all the other functions of ¢ in Eq. (9.3) are deter- 
mined. For small A, previous discussions show that c 
is given approximately by w(c) ~ 0.580, and also that 
3(c) = 3.22, u(c) 


1.48 (see Fig. 1). 


We now have 
R, = K} [1 (Ma) | ca (kA M,)} (9.4) 


where A depends on the given boundary-layer flow but 


is independent of M,. In order to determine the de- 


TABLE 2(a 


The Function vo(c) for Various Ratios of Wall to Free-Stream 
Temperature at a Free-Stream Mach Number of ().7 





TABLE 2(b) 


The Function v(c) for Various Ratios of Wall to Free-Stream 


Temperature at a Free-Stream Mach Number of 4.0 


The function v (c); M = 4, iE ee PuPy oO = 0.75 











Fic. 3. Curves of neutral stability at a free-stream Mach 
Number of 14, = 1.6 and a ratio of wall to free-stream tempera 
ture of 7/7; = 1.0738 
pendence on M, of the other quantities entering into 

R., we note that 
me ‘(ww — c)* F 
Ay = — dy = f[, (9.5 
wi 7 
is independent of M,, and, with the use of Eqs. (8.7b) 
and (9.2), that 


u—-L=u—-—L — (w,'c/T.~) (M;? — M,’) 9.6 


(9.6a 


Then, with the aid of the definition for k, Eq. (9.4) can 
be written as 

R, = (A + BM,?)/[IhV1 — M20 — 02] (9.7 
where 
(w,p’c)] (« — L) — My} + 


[77,/(1 — c)?]) (9.8a 


A=K((1 —c)* (7, 


B= K[(1 —c)? — Ih] (9.8b 


The quantities A, B, and c are independent of M, and 
are determined from the distributions w(y) and 7(y) 
of velocity and temperature, respectively, in the 
boundary layer. 


From Eq. (9.7) we find 


dR./dM, = | M°[2A(1 — c)? + B] — Aj {Me X 
[1 — M21 —c)?)"} (9.9 


Therefore 
dR,/dM, = 0 
when 
M, = 1/V2(1 — c)? + (B/A (9.10 
and 


dR. dM, <0 (9.11) 


when 
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a 
TABLE 3 
Wave speeds, wave numbers, and Reynolds Numbers for Neutral 
Subsonic Disturbances at a Free-Stream Mach Number of 1.6 
and Ratio of Wall to Free-Stream Temperature of 1.073 
(a) HM +06, M =1.6, T/T = 1,073 
1 b wl 
. 107" a Raxl07? R_x10~ 
¢ a ~ --) ay oa © 
- Lo) 
0 @) oo 0 oo oo aed 53 
35 .0515 40.7 -00534 42.1 512 4240 x 
okt -0642 21.3 00665 22.1 334 1110 = 
oe .0779 12.2 00807 12.7 23 3€5 >. 
25 110 4.82 0114 5.00 129 56.7 
oatS ola 3.23 0134 be 102 25.5 
325 178 1,63 0184 1.69 70.8 6.52 
- 3625 .231 1.09 0239 pe 61.7 2.92 
od t2 . 268 1,05 0278 1,09 68.8 2.70 
3625 .278 1.40 0288 1.45 94.7 4. 7¢ 
325 253 2.63 -0262 2.73 162 16.9 
eate - 208 6.23 0216 6.45 317 94.8 
25 . 186 9.89 0193 10.3 450 239 
2 144 28.2 -0149 29.2 988 1938 
0 0 eo 0 co aad aad 
i0 20 30 40 50 6 ? a 
@ in degrees 
able 2 (Continued Fic. 4. Minimum critical Reynolds Number versus direction 
% “ r/t of propagation of oblique waves 
or <x10~ a, 9x f a, x10™ 
o ao oo Pa, 
BAD 498) 290 5160 60 
ya a a 4 i 
325 2 23 2 7 : 
625 1295 0306 el? 
+ 6 fh TABLE 4 
35 | 35 15 Minimum Critical Reynolds Numbers and Direction Angles fot 
2 318 163 . , . . , : . 
3 284 227 Oblique Waves at a Free-Stream Mach Number of 1.6 and Ratio 
4 = es : of Wall to Free-Stream Temperature of 1.075 
ts 1099 109 O10: tf 28900 : 
£53 748 : 7x1c? M = 1.6, T/T 1,073 
3 718 2 x 9x08 i 
oo oo © 
olY 
Direction angle = cos M,/M, 
~ al~ -6 
} cos ~M,/M R, x10 
2S Xo 
Table 3 (Continued) 
0 90° Co 
(c) H = 14,4) 1.6, T/T, = 1.073 — 86.4 > Pe | 
2 82,8 19.2 
. e e 
C a xi0~* A, ¥gx10 Oy Ryxl0~ 4 Tae? 5.23 
6 68.0 2.70 
a7 021 7, 8 0022 ‘ 80 2 ed 8 60,0 1,89 
2875 0219 by PB ~00227 17,80 92.0 72 
292 10421 8.63 "00436 8.94 88.8 182 20 SL. 3 1.65 
3 0662 S37 . 00686 5.36 83.7 65.3 
25 hae 2.33 -0129 2.41 70.9 Jue  e 030 49.9 1.64 
35 183 1.41 0189 1. 46 63.3 4.89 
3625 .21¢ 1,17 .0224 1,21 61.8 3.35 1.2 41.4 — 
375 +263 1,07 .0272 1.22 68.8 2.81 
3625 . 250 1.55 0259 1.61 94.9 5.87 1.4 29.0 E 3 
35 .225 2.09 10233 2.17 115 10.7 1.5 20.4 §.0 
25 - 164 4.01 0172 4.15 163 39.3 
3 0936 9.92 "00969 10.28 227 241 1.576 10 19.6 
292 0605 17.05 .00627 17.67 252 711 
2875 0318 34.5 100329 35.7 268 2904 1.594 ) 76.7 
2857 «0 oo i) oo 00 0 0o 
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M, < 1°V2(1 — c)? + (B/A) 


M, < A, it follows that 


Since OU 


WW, V20 — cc)? + (B/A) < 1 (9.12) 
is the condition that all three-dimensional disturbances 
have a higher minimum critical Reynolds Number than 
two-dimensional ones. Thus, if it can be shown that 
the parameter C = J1/,V 2(1 — c)* + (B/A) for a 
specified boundary-layer flow is less than unity, then 
three-dimensional disturbances are more stable and are 
not very important in the stability problem. 

In Table 5 and Fig. 5 are given the results of calcula- 
tions of C as a function of 1/7, for an insulated flat plate 
fom WV, = 0 to VW, = 2.4 and of C as a function of 
T, T,; at a Mach Number of M, = 0.7 from 7,,/7; = 
0.70 to 7, 7 = 1.25. For the insulated plate, it is 
clear that three-dimensional disturbances are more 
stable than two up to a Mach Number of about 1.8 
The fact 


that C decreases and again becomes less than unity when 


when they become slightly more unstable. 


JJ, increases beyond 1.8 may not be of much signifi- 
cance because of the large errors at such Mach Num- 
bers due to the approximations used in the deductions. f 


10) THe INFLUENCE OF SURFACE COOLING ON 
BOUNDARY-LAYER STABILITY 


The problem of boundary-layer stability at super- 
sonic speeds is different from that at subsonic speeds 
in one important respect—-there is always a class of 
oblique waves (those with 1 < M, < M,) for which 
the values of ¢ are bounded away from zero (i.e., ¢ > 
1—1/M,>0). 


in the case of surface cooling. 


This leads to interesting consequences 


Let us first consider strictly two-dimensional dis- 
turbances (for which MM, = .1/,), since the situation is 
then somewhat simpler. The minimum value of ¢ on 
the neutral curve is 1 — 1 1/4), and the maximum value 
is given by w(c) = W;,,(c)—that is, it is the value of c 


corresponding to the intersection point of the graph 


of w(c) with that of W;,,(c) (see Figs. 2). For suffi- 
ciently small 7,,/7;, w(c) > 0 for c = 1 — 1/M, 
and v%(¢) is monotonically increasing for c > 1 — 1/M,. 


Also, further decreases in 7/7; raise the whole of the 

curve higher in the (c, %) plane. This behavior of 
the % curve is typical for Mach Numbers up to 1/, = 
‘48.[ For Mach Numbers greater than 4, = 7.48, 
no amount of cooling is sufficient to raise the v curve 


sufliciently to make w(c) positive atc = 1 — 1/M,. 


T It should be emphasized that, because of the use of the ap- 
proximate relations (8.16a) and (8 16b), the results of this sec- 
“ion are in most practical cases valid only for Mach Numbers 
hot too much larger than unity. Their implications in regard to 
the influence of surface cooling at low Mach Numbers are dis- 
cussed at the end of the next section 

+ The actual value of this upper limit to the Mach Number 
depends on the value of the Prandtl Number and the particular 
Viscosity law used. The value 1/, = 7.48 corresponds to a Prandtl 
Number of 0.75 and a linear viscosity-temperature relationship. 
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We see, therefore, that in this range of Mach Num 
bers, for 7,7), small enough, the range of values of ¢ 


on the neutral curve is bounded below by | 1/M, 
and above by an upper limit given by w(c) = W;,,(c), 
which can be made as close to 1 — 1. J; as desired 
by making 7,,/7; smaller. Thus, there is a critical 
temperature ratio (7), 7)), given by w(1 1/M,) = 
V,,(1 — 1/4/,) for which the whole neutral curve 


This can be seen by solving Eq. (S.16a) 
for a and noting that a—~ Oasc— | 1/\/;. There- 
7, — (7, 7)),, all values of c on the 


disappears. 


fore, when 7), 


TABLE 5 


The Function C = 1vV 2(1 c)? + B/A for Various Mach 
Numbers and Rates of Heat Transfer 


Te 7 « 
Table 5(a 





The function C = My 2(l-c)© + B/A for zero heat transfer; 


- n . al? - e 
The function C = M V 2(l-c)* + B/A for various temperature 
4 














ratios at M = 0.7; O-l, om=l 
7 / 0.70 Rg 9 on 5 
Ty qT 0.7 80 0.90 1.2 
C 0.871 797 0.722 0.544 
1.0 a 
0.8 + + + 
0.6 + + 
g 
04 + + + + 
02 + + + + 
1°) | 
04 0.8 1.2 1.6 2.0 24 
M, 


THE FUNCTION C=M 2(i-c® +8 FOR ZERO HEAT TRANSFER AT 


VAR 
10US MACH NUMBERS FIGURE 5(a) 
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THE FUNCTION C=My/2(i-cl?+2 FOR VARIOUS TEMPERATURE RATIOS 


AT A FREE-STREAM MACH NUMBER OF M =07 


FIGURE 5(b) 
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TABLE 6 
Critical Temperature Ratios for Complete Stability with Respect 
to Two-Dimensional Disturbances for Free-Stream Mach 
Numbers from 1 to 8 


o= 0.75, pPPp* Pu) 
T/T 


582 
. 766 
827 
1,018 
1.262 
1.453 
1.590 
1.678 
1.700 
1.641 
1.548 
1.146 
842 
~ 453 
- 240 
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neutral curve approach | — 1/1/,, so that all values 
of a Oand all values of R> o. 

It is to be noticed that in this limiting situation rela- 
tions (8.16a) and (8.16b) become exact since the errors 
are ()(a@"), so that the relation 


v1 — (1 M,)| = V;,[1 — (1 M,)] (10.1) 


for the critical temperature ratio (7°,/7)), is exact. 
When 7,,/7, < (T,/T;), evidently w(c) > W;i,,(e 
for c > 1 — 1/M,, so that no solution is possible, at least 
for small a. The inequality 7,,/7; < (7,/7)), is 
usually interpreted as the condition for the complete 
stabilization of the boundary layer with respect to 
two-dimensional disturbances.t The preceding argu- 
ments are essentially the same as those used by Lees.’ 


Calculations have been carried out to determine the 
variation of (7),/7)),. with J/, for the boundary layer on 
a flat plate with Prandtl Number equal to 0.75 and vis- 
cosity coefficient directly proportional to absolute 
temperature. The results are tabulated in Table 6, 
and the curve of critical temperature ratio against 
Mach Number for two-dimensional disturbances is 
shown in Fig. 6. Since the critical temperature ratio 
approaches zero at 1/,; = 7.48, no amount of cooling is 
sufficient to completely stabilize the boundary layer 
with respect to two-dimensional disturbances when 
M, = 748. The present numerical results do not 
differ appreciably from those of Van Driest,'! even 
though corrected viscous solutions [which enter into 
the determination of \(c) and thus of W;,, (c) | are used. 
The explanation is probably traceable to the fact that 


t This interpretation will be discussed later. 


w(c) is so extremely sensitive to variations in 7), 7 
that small errors in W;,, in Eq. (10.1) do not affect the 
solution 7), 7) much.f 

The conclusions depend, of course, strongly on the 
value of the Prandtl Number and the particular vis 
cosity-temperature relationship used.'! However, for 
the same reason as given above, it is not likely that th 
use of the present viscous solutions would significantly 
alter the results for these other cases obtained in refer 
ence 11. 

Now consider oblique waves with direction angle ar 
cos M,/\M,, where 1 M, < M,. Just as in the case 
of strictly two-dimensional disturbances, there is again 
the possibility of completely stabilizing the boundary 
layer with respect to oblique waves propagating in this 
direction for small enough 7), 7). The arguments 
carry through much as before except for one complica 
tion due to the fact that when the Mach Number is 
large the curve %(c) against c may have an inflection 


point and thus a maximum and minimum in the inter 


val0 <c 1 (see Fig. 2b). Thus, even when 7°,./7; is 
small enough for the relation 
%(1 — (1/M)] = ¥i,[1 — /M)] (102 


to hold, the corresponding w(c) may not be monotoni- 
cally increasing for c > 1 — 1 M4, when M, is small 


t Van Driest’s computations, in addition to the use of the in 
accurate A(c) of reference 1, also involve what is essentially an 
approximation for ¥j;,,(c) in terms of A(c). It is interesting to 
observe that if the exact procedure of calculation is used, in com 
bination with the inaccurate value of \(c) given by reference |, 
several critical cooling ratios will be found for each value of 
c= 1 1/.\,. This has been pointed out by Bloom.'* How 
ever, with the present function \(c) whose influence is small, such 


difficulties do not arise 
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enough.+ There could still be values of ¢ > 1 —- 
VW, for which w(c) < W;j,,(c) (see Fig. 2b), so that a 


neutral curve might still exist. We shall make the 


tentative assumption} that the critical temperature 


ratio, below which the boundary layer is completely 
stable with respect to oblique waves with direction 
are cos M, ‘A, is the largest 7), 7 for which 


ingle 
( > 8 ( ad l | M 10.3 


[his value of 7, 7) is evidently the one for which the 
corresponding v%(c) Curve is just tangent near its mini- 
mum point to the W;,,(¢) curve (see Fig. 2b). 

Obviously if 4, < 1 the minimum value of ¢ on the 
neutral curve is justc = 0. Therefore, 1t is not possible 
by any amount of cooling to completely stabilize the 
boundary layer with respect to oblique waves with 
direction angles greater than are cos | .\/;.. The limit- 
ing directions are thus the normals to the Mach lines in 
a plane parallel to the solid surface. 

For the boundary layer on a flat plate at a Mach 
Number .1/; 


ind a linear viscosity temperature relation, the critical 


} with Prandtl Number equal to 0.75 
temperature ratio 7), 7) has been calculated as a func 
tion of direction angle are cos M, \/,; and plotted as a 
full-line curve in Fig. 7. The calculations are carried 
out by first finding 7, 7; as a function of , from 
Eq. (10.2). The resulting values of 7), 7; when 
plotted against are cos M, Ml, give a curve with a 
maximum and a minimum for JM, in the range 0 - M, - 
,. The proper curve for the critical temperature 
ratio is then obtained by drawing a horizontal tangent 
to this curve at its minimum point (thus cutting off 
the dashed part around the maximum point in Fig. 


—\ kx 


i). The numerical results are listed in Table 7. 
Figure 7 shows that for a temperature ratio 7,, 7) = 
400 all two-dimensional disturbances are completely 
stabilized, and, for a somewhat smaller ratio 7, 7) = 
74, the three-dimensional disturbances corresponding 
to oblique waves with direction angles from 0° to 
(4.3° are completely stabilized. The critical tem- 
perature ratio 7), 7, decreases to zero rapidly, how- 
ever, when the direction angle increases from 74.3 
to are cos | M, = Fae 


When the temperature ratio 7), 7 is as low as 1.474, 
therefore, all oblique waves are stabilized except those 
with direction angles greater than 74.3 Estimates 
of the minimum critical Reynolds Number of the latter 
all of which have extremely small values of the wave 
speed c) show that R. ~ 10° or R,. ~ 10.'* Since such 
large values of the Reynolds Number are often far 
beyond the range of practical importance, it follows 
+ This difficulty does not occur when M, is as large as Jj, 
east for the boundary layer on a flat plate 
+ This assumption seems the most reasonable one to make, but 
vhether or not it is always justified will depend on results of 
calculation of neutral curves for large values of @ See com- 
ments further on in this section 
Chis is equivalent to the use of the condition mentioned in 
onnection with relation (10.3) 


dashed par 


The exact significance of the 


t of the curve is not yet entirely understood 
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that if the wall temperature ratio can be reduced to as 
low a value as 1.474 the boundary layer will be highly 
stabilized from the practical point of view. However, 
since the boundary layer is sensitive to changes in 7), 7 
even a small increase above 7), 7) 1.474 might make 
the boundary layer relatively unstable 

For each value of .1/; up to J/ 7.48 it 1s probabk 
that the nature of the oblique waves is about the same 
as at AM, 


ature ratio 77, 


{; that is, for a sufficiently small temper 
iy, t & possible to completely Sta 
bilize oblique waves with direction angles less than are 
cos | .J/,; but not those with larger direction angles. 
At the present time the state of knowledge regarding 
the general problem of boundary-layer stability at high 
Mach Numbers is mainly restricted to such results as 
the preceding regarding complete stabilization by 


cooling. It is only in such limiting cases, for which a 
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Fic. 7. Critical temperature ratios for complete stability with 
respect to three-dimensional disturbances at a free-stream Mach 
Number of J/,; = 4 


ABLE 7 
Critical Temperature Ratios for Complete Stability with Respect 
to Three-Dimensional Disturbances at a Free-Stream Mach 
Number of 4 
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TABLE 8 
The Functions Mc), Yin(c) for pp = pupw, ¢ = 0.75 and Zero 
Pressure Gradient 

c A(c) 94,,(¢) c A(c) 95,6) 

O 0 ” 58114 Pe 0273 ~61307 
el . 11520x107 . 58216 pa ~OL62 63625 
“a -9237x10-3 . 58251 8 0769 267622 
3 3, 136x10-3 - 58529 9 .1298 «75297 
4 0068 - 58891 95 shtto 83165 
3 0152 - 59868 1.0 2320 edo 


approaches zero in the limit, that the present calcula 
tion procedures are satisfactory. 

It should be emphasized that for high Mach Numbers 
the proper interpretation of the critical temperature 
ratios calculated by present methods is not completely 
certain. All that is known in any particular case is that 
for a temperature ratio s/ightly above the critical value 
a neutral curve exists at all points of which a is 
very small, and, when the temperature ratio is re- 
duced to the critical value, this neutral curve recedes 
to infinite values of RX with all values of @ on it approach- 
ing zero at the same time. It is possible that if we start 
with a temperature ratio considerably above the critical 
value and allow it to approach the latter the neutral 
curve might split into two branches, one of which dis- 
appears in the limit as described above, while the other 
continues to exist at finite values of & and R. In order 
to determine whether or not this is so, a method would 
have to be devised to allow the calculation at high 
Mach Numbers of neutral curves for values of @ which 
are not very small. 

For Mach Numbers not too much larger than unity, 
the situation is clearer, since neutral curves can be cal- 
culated by existing methods for finite values of @ and 
thus for all temperature ratios. It appears that in 
these cases the possibility suggested above does not 
occur and the entire neutral curve disappears in the limit 
so that the critical temperature ratios actually do cor- 
respond to complete stabilization of the boundary layer 
with respect to the particular class of disturbances under 
consideration. 

For low Mach Numbers (up to about 1/7, = 2), the 
theoretical results of Section (9), which are then satis- 
factory approximations, also provide considerable use- 
ful information. The numerical calculations indicate 
that the parameter C increases as 7,,/7 decreases (see 
Fig. 5). Thus, the effect of cooling a flat plate at a 
given Mach Number is evidently to make three-dimen- 
sional disturbances relatively more unstable, so that 
with sufficient cooling some of them can actually be- 
come more unstable than two-dimensional disturbances. 
Since C for an insulated plate increases from C = 0 at 
JJ; = 0 to a value close to unity when J1/, =~ 2, 
it appears that the amount of cooling required to 
make any of the three-dimensional disturbances more 
unstable than two-dimensional ones is extreme at low 
subsonic Mach Numbers but gradually becomes less 
and less extreme as J/, increases, until, at a Mach 
Number of the order of two, little is required. In par- 
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ticular, at /, = 0.7, the wall temperature ratio must 
be reduced below 7), 7, = 0.58 (see Fig. 5). 

It is known that at subsonic Mach Numbers cooling 
of the solid surface is effective in stabilizing the bound- 
ary layer if only two-dimensional disturbances are con- 
sidered.” The preceding remarks indicate that prob- 
ably this conclusion is still valid even if three-dimen- 
latter 
Reynolds 


sional disturbances are considered, since the 
usually have a higher minimum critical 
Number than two-dimensional disturbances 


under conditions of extreme cooling). 


except 


Even at the moderate supersonic Mach Numbers 
when three-dimensional disturbances begin to play the 
leading role, cooling still has an important stabilizing 
effect, although complete stabilization at supersonic 
Mach Numbers as predicted in reference 2 is impossible 
in the strict sense, according to the present theory. For 
example, estimates show that for an insulated flat plate 
at J, = 1.6, R,, is of the order of 10°, whereas the cal 
culated results show that for a cooled plate at the same 
Mach Number 1/; = 1.6, with a temperature ratio 
T/T, = 1.073 which is just small enough for two- 
dimensional disturbances to be completely stabilized, 
R,. is of the order of 10°. 


(11) CONCLUSIONS 


(a) An order of magnitude analysis of the complete 
linearized equations for a three-dimensional disturbance 
superimposed on a two-dimensional boundary layer 
shows that they can be reduced to much simpler equa- 
tions in the first approximation. These simplified 
equations are still valid at high free-stream Mach 
Numbers if correct velocity and temperature profiles 
for the boundary layer in such cases are used. At 
very high Mach Numbers, however, the accuracy 
of a first approximation may not be great, so that higher 
order approximations (and therefore more complicated 
equations with more terms retained) may have to be 
considered. 

(b) The present equations for a three-dimensional 
disturbance transform to equations that have the same 
form mathematically as those for a two-dimensional 
disturbance but which do not, however, correspond 
to a real two-dimensional disturbance. Although no 
direct conclusion regarding the importance of three- 
dimensional disturbances is now possible, the use of 
these transformed two-dimensional equations leads 
to a considerable mathematical simplification of the 
theory. 

(c) Corrected viscous solutions have been obtained 
which lead to a new formulation of the stability theory. 
The conclusion reached by Lees and Lin that the 
relation for the characteristic values is independent 0! 
the temperature fluctuations and the boundary condi 
tions imposed on them—is now found to be restricted 
in general validity to subsonic and slightly supersonic 
Mach Numbers and is true only in a restricted sense 
namely, for very obliquely traveling waves at high Mach 


Numbers. At high Mach Numbers the relation for the 
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characteristic values depends in general on the thermal 
boundary conditions. It is found, however, that for 
very small values of the wave number the relation 
again has the same simple form as at low Mach Num- 
bers. In all cases at high Mach Numbers when the 
wave number is not small, the relation is more complhi- 
cated. The detailed mathematical development of this 
improved theory is to be presented elsewhere. 

d) Consideration of the present analytical and nu- 
merical results suggests that probably three-dimen- 
sional disturbances are of little importance at subsonic 
free-stream Mach Numbers, since they usually have a 
higher minimum critical Reynolds Number than two- 
ones, possibly under conditions 


dimensional except 


of extreme surface cooling. However, further numeri- 
cal calculations for oblique waves in this range of Mach 
Numbers would be useful. 

e) As the Mach Number increases, three-dimensional 
disturbances enter the picture under conditions which 
become less and less extreme, until finally for free- 
stream Mach Number .J/; between one and two they be- 
gin to play the leading role in many problems of prac- 
tical interest. 

(f) At supersonic free-stream Mach Numbers, the 
boundary layer can never be completely stabilized with 
respect to a// three-dimensional disturbances. Al- 
though oblique waves whose directions of propagation 
are at an angle less than are cos | .J/; with the free- 
stream direction can be completely stabilized by suth- 
cient cooling when .J/; is not too large, those whose 
directions of propagation are at a greater angle cannot. 

(g) The graph of critical temperature ratio (below 
which strictly two-dimensional disturbances are com- 
pletely stabilized) versus Mach Number (Fig. 6) is not 
significantly different from that found by Van Driest,"! 
in spite of the corrections to the viscous solutions. 
Although these corrections are very important for 
the accurate determination of the curve of neutral 
stability at supersonic Mach Numbers, they are not as 
important for the determination of the critical tem- 
perature ratios. 

(h) For Mach Numbers up to about two, cooling of 
the solid surface is effective in stabilizing the boundary 
layer. For higher Mach Numbers, reliable calculations 
of the minimum critical Reynolds Number over wide 
ranges of the wall temperature ratio have not yet been 
made; but this general conclusion would probably 
remain unchanged for Mach Numbers up to six, say. 

(i) Little more can be said about the problem of 
boundary-layer stability at Mach Numbers greater 
than two, since the present methods of numerical cal- 
culation for the curve of neutral stability (and thus 
Number) do not 
appear to be adequate in such cases. When ¢ and 4 
are both large, there are difficulties in working with 


the minimum critical Reynolds 


the more complicated relation for the characteristic 
values and also in evaluating the inviscid solutions of 
the disturbance equations. - These problems are being 


Studied further. 
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Numbers there is an 
even than the 
namely, that the present asymptotic theory, which is 


(j) At Mach 
more serious difficulty 


very high 


preceding 


essentially a first-order theory, may in itself be inac 
More 


information about this matter will be provided by nu 


curate [as pointed out in paragraph (a) above | 


merical results (including estimates of certain param 
eters) for high Mach Numbers, which can be obtained 
as soon as the difficulties mentioned above in paragraph 
(1) are resolved. 

(k) It is believed that the approximate method of 
calculation of the neutral curve described in Section (S 
is adequate up to a Mach Number of the order of two, 
especially in cases involving cooling, when c and @ are 
Additional neutral curves for the oblique 
Maeh Numbers or perhaps 


both small. 
waves in this range of 
further calculations of the type described in Sections 
(9) and (10) would be of considerable value. It might 
also be of interest to determine in a few typical cases 
the effect of using a more accurate approximation than 
Eq. (S.16a) by keeping terms of higher order in a, as 
Lees has done.” 

(1) For practical purposes, the transition Reynolds 
Number is much more important than the minimum 
critical Reynolds Number. Thus, information con 
cerning the rate of amplification of the boundary-layer 
disturbances after they enter the unstable region would 
be extremely useful. Further investigations are being 
directed toward developing the amplification theory 
for three-dimensional disturbances, which may _ be 
somewhat different from that for strictly two-dimen 
sional disturbances. 

(m) Applications of any of the present theory to cases 
involving finite wall curvature should be made with 
caution. Although there are no restrictions regarding 
pressure gradient, the basic equations derived here are 
strictly correct only for flow over a flat surface. Al- 
though curvature effects are usually extremely small 
for moderate values of the direction angle of the dis- 
turbances (see references 6 and 15), they may become 
important for values of this angle approaching 90 
For example, whereas the present theory predicts an 
infinite minimum critical Reynolds Number for dis- 
turbances with direction angle of 90°, it might be ex- 
pected from the work of Gortler’® that for small but 
finite concave curvature such disturbances would be 


unstable at finite Reynolds Numbers. 


Appendix 


DERIVATION OF THE RELATION FOR THI 
CHARACTERISTIC VALUES 


The basic relation (7.9) is obtained as the practical 
form of the secular equation on which all subsequent 
deductions and calculations are based. A complete 
derivation involves a detailed study of the asymptotic 
solutions of the system of Eqs. (6.9) to (6.13) and will be 
presented elsewhere. Here we shall merely outline 
the major steps and point out the nature of the approxi- 
mations involved. 
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For large values of aR and finite values of (y — y,), 
formal expansions give rise to solutions of the following 
three types: 

(a) Two solutions for which each of the functions 
has an asymptotic form 
=Z, (y 


Zily + [1/(aR)]Z;""(y) +... (Al 


(b) Two solutions for which each of the functions has 
an asymptotic form 
Pj} |Z m a } 


(4m)~ "Zz," 4 al | 


P(y ac - ¢)/v]'" dy | 


Zi(¥ lexp [+ (aR) 


(c) Two solutions for which each of the functions has 
an asymptotic form 


pexp | (aR) O}} [Z a a 
(aR) / 


Zily | 
Be Hed 
Qy) = | [o(w — c)/v] “dy | 


x 
« 


In the above, y, is the point where w = c. The ap- 
proximate solution Z;"’(y) is the inviscid approximation 
and has the exponential behavior exp | * By| for large 


( 


y, as discussed in Section (7). 
Consider the case of real c. It is clear that, if the 
solution is to remain bounded at infinity, only the 
negative signs are allowed in a// the exponentials. 
Thus, we have one permissible solution from each of the 


three types, (a), (b), and (c). Let us denote these 
solutions by 

(a) (f, @ 8) = (F, ®, 0) 

(b) (f, &, 8) = (fs, 3, 43) 

(c) (f, ¢, 0) = (fis Psy 95) 


To satisfy the conditions at the wall, i.e., 
f=¢=0, a’+bdb6=0 (A4) 
where a and / are given constants, we must have 
F(O) f3(0) f5(0) 
+ hb (0) o3(0) 6;(0) = O 
O(0) @;(0) 45(0) 


F(O) f3(0) 5(0) 
a P(0) o3(0) $5(0) 
0’(0) 63’(0) 6;’(0) 


(A5) 


Actual calculations reveal that the general orders of 
magnitudes of the functions are given by [see Eq. (A9) | 


CF @;,0) oo (1,1, 1) 


(fa, bs, 03) ~ (C1, €, O) (A6) 
(fs, D5, 45) ~ (0, é, 1) 
to the order of « = (a&R)~' Hence, we obtain 
(0) F(O) = [¢@3(0) f:(0)] (1 + A (AZ 
where 
ado’(0) + b0(0 4;(0) ;(0) f3(0) 
= ; (AS 


F(O) aé;'(0) + b0;(0) @;(0) o3(0 
We shall consider A in more detail later and just ob 
serve at this point that it is always 0(1) at most. 

From Eq. (A6), it follows that ¢;(0)/f;(0) = O(e), 
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which shows that in the limit « > 0 Eq. (A7) reduces t 
the inviscid formulation (0 0, as it should 

On the basis of the general nature of the solutions 
Al), (A2), and (A3) 
nitude of the solutions as shown in Eq. (A6), we may 


and because of the orders of mag 


refer to these solutions as (a) the inviscid solution, (b 
the viscous solution of velocity variation, and (c the 
viscous solution of thermal variation. 

For practical purposes it is found to be desirable t 
use asymptotic forms of the viscous solutions which 


remain accurate even for (y ¥.) > 0. For exampk 
for the viscous solution of velocity variation, detailed 


studies give 


AQ 
where 
Biy i(“ - ‘) E (‘ " ‘) dy ly ’ 
v 2S pv om 
A10 
and 


Al 


This solution agrees asymptotically with Eq. (A2) wher 
¢ is large but is also valid for small ¢. Furthermor 
if each of the integrals (A10) and (A11) is evaluated in 
the neighbourhood of y, by retaining only the first tern 
in the power series expansion of the integrand, the solu 
tion given in reference 1 is obtained. An improved 
version of the viscous solution of thermal variation may 
also be developed. * 

By substituting Eq. (AQ) into Eq. (A7 
The only additional stey 


, the relation 
(7.9) can be easily obtained. 
is to express ®(y) as a linear combination of two linearl\ 
independent inviscid solutions ¢(y) and @»(y)—-ie 
two such solutions of Eq. (7.3)—-and then to requir 
that this linear combination shall be proportional t 
exp [—8y] for y > 1, where y = 1 corresponds to the 
edge of the boundary layer. 

Now let us return to Eq. (AS) for A and investigate 
the conditions under which it may be neglected and 
those under which it must be retained. It may be 
shown that, for 7,, bounded away from zero,t Eq 
(AS) becomes 
To) ((al(Z 0’ — Tu@y’)/Tol} 4 
b0:(0) |i x(s, @ Al2 


A = [(y — 1) 
bc) } {[0;(0) | /[aé@;’(O 


where x(2, o) is a function depending on the viscous 
solutions which is always 0(1) (z is defined by Eq 
* The improved solutions are obtained by extending a metho 
of Tollmien in the incompressible case.8 Their use does not 
alter any of the order-of-magnitude estimates in the present dis 
cussion 


7 The exceptional case 7, — 0 is to be considered latet 





whi 


whi 


prece 
of th 
negle 

In 

l4 
pend 
tions 
is ob: 
of Ec 


guish 


In « 
Vac h 
neglec 
eter qd 

It n 


the f Tr 


where 
basis ( 
been e 


(A13), 


Ices to 


utuons 
| mag 
e¢ may 
mn, (b 
C) the 


ble to 
which 
imple, 


tailed 


AQ 


A10 


when 
more 
red in 
term 
S( lu 
roved 


| may 


ation 
| step 
early 

1.€. 
quire 
al to 
o the 


igate 
and 
vy be 


Eq 





STABILITY OF THE 


7.11) and o is the Prandtl Number). It might be 
mentioned that the factor c in the first term is due to 
¢,(0) being proportional to c. 

In reference | the conclusion was reached that the 
characteristic-value problem is independent of the 
temperature fluctuations and the boundary conditions 
imposed on them, to the order of accuracy of the theory. 
From Eqs. (A7) and (A12) it is clear that this will be 
true when A can be neglected without any restrictions 
being placed on the values of a and 6. Since the factors 
multiplying (y — 1) M,2c 7), in Eq. (Al2) are all0(1 
at most, it follows that the conclusion is valid in two 


ain cases 


c< 1 (A13 
which was considered in Reference 1, and 
M? <1 (Al4 


which corresponds to oblique waves with direction 
angles close to 90 

Case (A13) covers most practical situations at sub 
sonic and slightly supersonic Mach Numbers, when ¢ 
is usually not large. Case (Al4) extends the conclu- 
sion to very obliquely traveling disturbances at high 
Mach Numbers. 

It might be remarked that since ¢;(0) is proportional 
to c and since the theory of reference | was restricted 
to values of c of order ¢, it is understandable that the 
preceding conclusion should be a natural consequence 
of the asymptotic theory [in which quantities O(«) are 
neglected |, as developed in reference 1. 

In all situations at high Mach Numbers when Eq. 

\l4) does not hold, the exact form of Eq. (A7) de- 
pends on a knowledge of the thermal boundary condi- 
1.e., a knowledge of the constants a and ). If it 
4;(0) = O(1 €), an examination 


tions 
is observed that @;’(0) 
of Eq Al12 


guished 


shows that three cases may be distin- 


b/a = O01); A = Ole Ald 


A 1 1)M,2c?/T.,| ; 14;(0) | [(a/6)65’(0) + 


4;(0) t x(Z, o A16 


4 ™ I(4 oy M2 ' ey (3, o (A17 


In case (A135), A may again be neglected, so that Eq. 
A7) again reduces to the same simple form as at low 
Mach Numbers. In the other two cases, A cannot be 
neglected, and in Eq. (A16), moreover, a new param 
eter a 6 enters into the theory. 

It may be noted that when A is not 0(e), it is always of 


the form 
M,°c" fai A, A18 


Where A, = O(1 


basis of small c is somewhat smaller than might have 


Thus, the error in neglecting A on the 


been expected from the discussion regarding condition 
{(A13 
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To determine the proper form of the thermal bound- 
ary condition (A4) and thus the values of the constants 
a and 6 in practical situations, an analysis of the heat 
conduction problem in the medium below the interface 
between the gas and solid wall is necessary. Since 
this is straightforward but fairly lengthy, we shall just 
give the final conclusions here for a few typical cases. 
Consider a boundary-layer flow of a gas extending 
from the solid surface at y = 0 to y = ~, and suppose 
that from y = 0 to y = —h, isa thin metal wall, from 
: —ha channel through which a cooling 
liquid is flowing, and finally from y = —h to y = : 
Although the problem is con 


y= —-htoy = 


an insulating material. 
siderably idealized, it illustrates some of the essential 
features in certain cases of interest. For this problem, 
it can be shown that the condition for the temperature 
fluctuations in the boundary layer is given by 


6’(0) — jw — [(1-k) (dk dT) (OT Oy)),{0(0) = 0 
Al9 
where 
k, Ww 
an at f(r 4 52H) 
k { k Ww 
(ste V [le ste) 
I1-—-S e 1 + + 
Rn w a Rw 
kw 4 
(1 5 Je a A20 
RW, f 


in which the various quantities appearing are defined as 


follows: 
| kw | ( se) ww (I ) 
_= ' 1 + e + 
A k WwW) ie Rw, 
(NY (i+ tem 
Rw, J Rw 


(A2] 


—1ar ( PinCn R, (A22 


Wm" —_ 
The 
symbols p,,, ¢,, and k,, denote density, specific heat, 
The sub 


with corresponding definitions for w, and «a. 


and thermal conductivity, respectively. 


scripts 0, m, /, and 7 refer to gas at y = O, metal, 
The quantity ac 


the fluctu 


liquid, and insulator, respectively. 
is the nondimensional time frequency of 
ations {see Eq. (6.1) ]. 

Relation (A20 
reduces to 


is appropriate for a cooled wall. It 


WD = @, (R,, k (A23 
for high frequencies, and to 
ao = w(R k (A24 


for low frequencies. The case of low frequencies is of 


most interest here. In Section (10) it is shown that if 
the simple eigenvalue relation (7.14) obtained by 
neglecting A is used, then in the limiting case of com- 
plete stabilization a@ — 0 everywhere on the neutral 


curve. Since the frequency is given by ac = ac cos @, 
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where 0 is the direction angle of the waves, it follows 
that w — 0 in the limit and thus from Eq. (A15) that 
A is indeed negligible. Thus, we see that the use of the 
simple relation (7.14) in calculating the critical temper- 
ature ratios is justified. 

The definition of w for an insulated wall is easily ob- 
tained from Eq. (A20) by letting 4; = h. For high fre- 
quencies this also reduces to Eq. (A253). If there is no 
metal layer so that the insulator extends right to the 


surface, h = 0, and w is given by 
a = @ (k ky) (A25 


for all frequencies. 

In all cases at high Mach Numbers when a neutral 
curve exists, a = O(1), c = O(1), are reasonable esti- 
mates. Then, on taking cos 0 = O(1), we find that Eqs. 
(A23) and (A25) are accurate approximations and that 
w is extremely large. This is true of Eq. (A23) even 
when & = 0(10~*), say. Thus, in most cases of interest 
at high Mach Numbers when cooling is not sufficient 
to give complete stabilization, w is extremely large, and 
the relation for the characteristic values is that corre- 
sponding to Eq. (A17). 

The exceptional case 7), > 0 has not been considered 
so far. This occurs in the calculations of the critical 
temperature ratios for Mj, ~ 1 and for 1, ~ 748 
(see Figs. 6 and 7). We may justify the procedure of 
using the simple relation (7.14) for these calculations 
a posteriori—that is, we assume Eq. (7.14) to hold in 
these limiting cases and then check the consequences to 
see if they are consistent with the assumption. In the 
case MW, > 1,c = 1 — 1 MM, > O, and an examination 
of Eqs. (10.2) and (8.16b) shows that 7),«c. Then 
Eq. (A12) shows that A is now O(c), so that the pro- 
cedure is justified. In the case 1/,; ~ 7.48, c remains 
finite, and relations (10.1) and (8.16b) imply that 
K,, = 0 -i.e., that [d dy(w’ T)]. = 0. This indicates 
that the solution obtained by the use of Eqs. (7.14) 
or (S.1) reduces to that obtained from an inviscid for- 
mulation. On the other hand, the more accurate rela- 
tion (7.9), when put into the form (8.1), has ®(z) re- 
placed by [1 — F(z) (1 + A)]~'. At the same time, Eq. 
(A12) shows that A = O(1/7,,). It can be shown 
that for A with such a behavior the more accurate rela- 
tion also leads to the inviscid solution in the limit 
7, 0. Thus, the results obtained by the use of Eq. 
(7.14) are again justified. 

Summarizing our main conclusions, we have found 
that for very low frequencies the boundary condition 
for the temperature fluctuations is 


0’(0) + [(1/k) (dk dT) (OT Oy) ], 0(0) = 0 (A26) 
and the corresponding secular ec,uation 1s 
E(a, ¢, M,*) = F(e) (A27) 


while for very high frequencies, the boundary con- 
dition is 


a0) = O (A28) 


and the corresponding secular equation is 


E(a, C, M,?) = F(z)} 1 + [(¥ — 1)M,?c?,T 
x(2, a); A29 


The boundary conditions (A26) and (A2S) are 
actually just what would be expected physically, as 
the following considerations will show. In the sta- 
bility theory, the temperature fluctuations are essen 
tially natural oscillations originating in the boundary 
layer. Although these fluctuations can penetrat 
into the adjacent medium (gas, liquid, solid, etc.), they 
must go to zero at large enough distances, For 
very low frequencies, it is clear that such fluctuations 
will go far into the surrounding medium before dying 
out. Thus, for low enough frequencies, the change in 
amplitude will be negligible over distances in the solid 
of the order of the boundary-layer thickness, so that the 
fluctuating part of the heat flux will be negligible 
Therefore, since the left-hand side of Eq. (A26) is pro 
portional to the latter quantity expressed in terms of 
boundary-layer properties at the wall, relation (A26 
will hold closely. 

On the other hand, for very high frequencies 
the thermal inertia of the medium will prevent the 
fluctuations at any given location from reaching the 
maximum amplitude attained near the origin of the 
disturbance. Thus, the amplitude will die away to 
zero within a short distance inside the solid wall, for 
example.* Ultimately, for high enough frequencies, 
even the amplitude at the surface will become negli- 
gible, so that the condition 6(0) = 0 will be a very 


accurate approximation. 
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* This ‘‘surface-layer effect’’ or ‘‘skin effect’’ is well known ™ 


other problems involving temperature fluctuations 
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List of Svmbols 


Reference 


Dimensional quantities Nondimensional quantities quantities 


Positional co-ordinates 


1) 4 x rs) 

2) x y 6 

3) x e 6 
Time 

$) 6/U; 
Velocity components 

5) m = tm + my wy) + f(y) exp [i(ax + Bs — act l 

6) ue = te + Uy’ oy) + ad(y) exp [i(ax + Bz — act)] U; 

7) u = i h(y exp li(ax + Bs — act 7 
Density 

8) p* = p* + p'* p(y) + rly) exp |[i(ax + Bz — act ~~ 
Pressure 

* 


9) p* = p* + p”™ Ply) + rly) exp [i(ax + Bs — act p 
Temperature . 
10) 7* = 7* +7" T(y) + Oy) exp [t(ax + Bs — act)] 7,* 


Viscosity coefficients 


11) w* = g* + pp" u(y) + (du/dT) 0(y) exp [i(ax + Bz — act ju * 

12) A* = v* +A” My) + (dd/dT) 0 y) exp [i(ax + Bz — act . 
Thermal conductivity 

13) R* = k* + k’* u(yv)/o(y)] + (d/dT) (u/o)0(y) exp [i(ax + 82 — act Crin* 
Wave numbers of the disturbance 

a ") 
l4) a ag 6* 
a x) 

15) a 3 fr) 
Wave velocity of the disturbance 

16) c* ( U7; 
Specific heat at constant volume 

17) C l C, 
Specific heat at constant pressure 

18) C, ¥ i 
Gas constant per gram 

19) R* y¥-1 C. 
Acceleration due to gravity 

20) g 1/F2 1.2/8 
Momentum thickness of boundary layer 

21) 0* = fo” pw(l — w) dxe 


Reynolds Number 
22) Based on boundary layer thickness 6 
23) Based on momentum thickness 6* Re = pi*U,0*/ii 


24) Based on distance x; from leading edge R, = ~i*U\x/ii* 
Mach Number 
or ~~ 
i) My = Ui/V y,R*1 
Froude Number 
26 F= U,/W 26 
Prandtl Number 
97 o¢ = Gu b* 
a) A bar | denotes mean value, a prime denotes fluctuation. The subscript 1 refers to the value at the edge of the boundary 


layer at the particular x-location considered, and the subscript w refers to wall value. The subscripts r and 7 denote the real and imagi 
ary parts of a quantity, respectively. 

b) The nondimensional mean quantities 1;(x1, x2)/Ui, ue(x1, x2)/U1, p(x, y), p(x, ¥), are represented as functions of y alone—that 
Is, as w(y), v(y), p(y), p(y), , respectively. This corresponds to an approximation depending on the boundary-layer nature of the mean 
flow [see Section (4)]. 

c) The quantities C,, C,, u, \, k, o, y are all considered to be functions of 7 only. 
by Q’ = dQ/dTT’. 

d) For moderate Mach Numbers and rates of heat transfer C,, C,, y, ¢ vary only slightly in the boundary layer 
regarded as constants equal to their free-stream values in the numerical calculations. 

e) The reference length 6 is the boundary-layer thickness at the particular value of x considered. Its definition is somewhat arbitrary 
and may be taken in any convenient manner. In the present numerical computations it is taken to be the value of x2 at which 4, /U, = 
0.999. 


The fluctuation Q’ of any such quantity Q is given 


They are therefore 








The Unsteady Forces Due to Viscous Wakes 
in Turbomachines 


NELSON H. KEMP* anv W. R. SEARS* 


Cornell University 


ABSTRACT 


The configuration of viscous wakes of cascade blades is ap 
proximated from single airfoil experiments. The unsteady force 
and moment on a downstream blade passing through such wakes 
is then calculated on the basis of the theory of isolated thin air- 
foils in nonuniform motion. The results indicate that the force is 
nearly proportional to the profile-drag coefficient of the upstream 
blades. For typical values of this coefficient and conventional 


cascade the unsteady forces arising from passage 


through viscous wakes are of about the same size as those due to 


geometry 


aerodynamic interference between the moving blade rows, pre 


viously estimated 


SYMBOLS 


b = axial distance between centerlines of stator and rotor 
rows (Fig. 2) 

b’ = axial distance between trailing edge of stator row 
and leading edge of rotor row 

b’ = b — ¢, cos a, — ¢, cos (B — a 

( = airfoil half-chord 

d = pitch of blades in direction of rotation of turboma 
chine 

q = dynamic pressure, (1/2)pV? 

t = time 

u = perturbation velocity component in the stream direc- 
tion 

v = nonsteady upwash induced at airfoil 

% = coordinate along airfoil, measured from mid-chord, 
positive toward trailing edge 

y = coordinate perpendicular to x, positive upward 

v’, y’ = oblique coordinate system (Fig. 2) 

Cp = profile drag coefficient 

C, = lift coefficient, L/(1/2)pV?2 

Cy = moment coefficient, 17/(1/2)pV?(2c)? 

G = coefficient of nonsteady upwash iaduced at rotor 
blade, Eqs. (28) and (29) 

K = dimensionless half-width of viscous wake Eq. (21) 

K,, = modified Bessel function of the second kind 

L = lift, positive upward 

M = moment about airfoil center, positive in the stalling 
direction 

R = real part 

S(w) = {tw[Ko(iw) + Ki(iw)}}~ 

U = rotational speed of rotor, positive downward 

V = speed of stream relative to airfoil 

) = half-width of viscous wake 


a = angle of stagger of blades (Fig. 2) 


8 = angle between stator and rotor blades 
v = circular frequency, radians per sec. 

p = gas density 

o = solidity of cascade, 2c/d 

w = reduced frequency, vc/V 
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Subscripts 


stator 
r = rotor 
7 = total effect due to all viscous wakes 
( center of viscous wake 
free stream 
Q = taken constant along rotor blade 


Superscripts 
rv = effects on rotor blade due to viscous stator wakes 


= stator 


INTRODUCTION 


- AN EARLIER INVESTIGATION! an approximate calcu 
lation was made of the forces and moments that act 
on the blades of a turbomachine due to their passage 
through the disturbed fields of flow of adjacent blade 
rows. The extreme complexity of the problem made it 
necessary to introduce several important approxima- 
tions. First, the calculation was restricted to the case of 
an elementary stage consisting of a stator and rotor. 
Second, the familiar approximations of unsteady thin- 
airfoil theory were adopted. Finally, the theory 
adopted was a first-order one based on the assumption, 
later confirmed, that the fluctuation of blade circulation 
is always small compared to the mean value. 

Within these approximations, calculations were car- 
ried to the fluctuations of the lift, 
moment, and circulation due to passage through the 
potential-flow field of the adjacent row of blades. For 
the downstream (‘‘rotor’’) row, this included the effects 


out determine 


of passage through the sheet-vortex wakes of the up- 
stream (‘‘stator’’) blades, although this is, to be sure, a 
‘‘second-order’”’ effect in the spirit of the approximations 
lift, in particu- 


mentioned above. The fluctuations of 


lar, were worked out for typical stage geometries and 
were found to be appreciable. 

In the work just described, no account was taken of 
the effects of viscous wakes. This is, of course, an ob- 
vious shortcoming of the investigation, since such 
wakes are always present in a real fluid, and one may 
even imagine that the viscous-wake effects are larger 
than the purely potential-flow effects that were con- 
The authors’ purpose, however, was to con- 
The present paper 
Their 


sidered. 
sider these phenomena separately. 
reports on a study of the viscous-wake effects. 
evaluation involves additional approximations, as will 
be explained, although the mathematical details are 
actually somewhat simpler than in the earlier task. 
In outline, the present work involves the following 


steps. 
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[he flow downstream of a cascade of blades with 
viscous Wakes will be represented by an inviscid shear 
fow. The velocity profile of the wake behind each 
blade is made to resemble the results of measurements 
behind isolated airfoils. 

2) The theory of thin airfoils in unsteady flow is ap 
nlied to calculate the force and moment resulting from 


passage through such a series of wakes—-1.e., gusts. As 
is often the case, the presence of vorticity in the oncom 
ing stream is neglected in making this application 

In view of the small-perturbation type of approxima 
tion that is inherent in thin-airfoil theory, it is consistent 
to add together the fluctuating forces, moments, and 
irculations calculated here and in the earlier paper.! 
Naturally, to be precise one would add these various 
terms in proper phase—e.g., superimpose them as 
complex numbers before taking real parts. Asan upper 
limit, however, one might simply add their absolute 
values. The absolute values are easier to compute and 
wre the numerical values presented in the graphs and 
tables of this report and the earlier one. 


DESCRIPTION OF THE VISCOUS WAKES 


The model adopted here for the viscous wakes of 
cascade blades is based upon the experimental results 
of Silverstein, Katzoff, and Bullivant,? who measured 
the wakes of airfoils and generalized their results. It 
will be assumed that the wake behind each blade is the 
same as the wake behind an isolated airfoil. Since the 
wakes are relatively narrow, compared to the blade 
spacing, in the region of the adjacent downstream row, 
this assumption seems justified. 

As a further simplification, the wake will be con 
sidered to consist only of an inviscid, symmetrical, shear 
perturbation of the undisturbed stream. It will be re- 
called that all other perturbations due to the upstream 
airfoils have been included in the previous study;' thus, 
ina small-perturbation theory, only the introduction of 
the wake disturbance into the free stream is required 
here. 

Silverstein, Katzoff, and Bullivant have shown that 
the half-width of the wake, Y, may be calculated from 


the following formula 


Y = O68V 3C, “clx/c — 0.7) (1) 


where c is the airfoil half-chord, x is measured along 
the wake axis (free-stream direction) rearward from the 
center of the airfoil (see Fig. 1), and Cp is the airfoil 
If g, is the dynamic pressure at 
that in the free stream, 


profile-drag coefficient. 
the center of the wake and q 
reference 2 also shows that 


l - (G-/q,.) = (4.84C p \/(x/co — 0.4 (2) 

Finally, the wake velocity profile is also given in the 
same reference in terms of the dynamic pressure g as 

l q/q ; T y ; . 

cos Gy (3) 


qo 2 Y 


IIA 


l q 


where y is the coordinate perpendicular to the stream 
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Fic. | Viscous wake of an airfoil according to reference 2 


(see Fig. 1). Eqs. (1), (2), and (3) completely define the 
wake. 

In order to convert the wake equations from dynamic 
pressure to velocity, the small perturbation assumption 
is used. Let uw denote the perturbation velocity in the 
x-direction, which is the only velocity perturbation 
considered, as has been mentioned above. Since the 
free-stream velocity is also in the x-direction, one has, 
assuming u< /, 


(q/q 1— [((V + 4)*/V 2(u/V) (4 


In view of Eq. (4), Eqs. (2) and (3) become relations for 
the perturbation velocity at the center of the wake, %,, 
and the wake velocity profile: 


ue/V = —(2.42Cp ”*)/(x/e 0.4) (5)4 


u/Ue = cos” [(r/2)(v/Y y 


In order to simplify the mathemat‘cal analysis, the 
cosine wake profile of Eq. (6) will be replaced by a pro 
file in the form of a Gaussian error curve that has the 
flow in the 


same momentum flow. The momentum 


wake is proportional to the square of the velocity, 


(Vv + u)*? = V*+ 2u. The momentum flow across any 
section therefore depends on fu dy taken across the 
wake, and it is this integral that will be made the same 
for the cosine profile and the assumed Gaussian profile. 


For the former the integral is 


a} 
Je 9 7 
| cos* -~}dy = } (7 
ow =% > 3 

written u/u. = 


The Gaussian profile may be exp 


—(ky/Y)*], where the constant k is to be determined 


by the momentum flow condition. The integral is 


+} 
( exp [—(ky/Y)?] dy = 
} 


0 


») 


y Var Vf 2 
e dy = 
b ~ b Vrs 


The expression in parentheses is the well-known error 
function denoted by erf (k), and comparison of Eqs. (7) 


reader will easily verify that the velocity perturbation 


+ The 
u- according to this formula may be so large 10 per cent of V 
just behind the airfoil as to cast doubt on the validity of the small 


perturbation approximation However, th« perturbation at 


tenuates downstream and rapidly in the cross-stream direction 
Eq. (6)], and it is therefore hoped that the approximation is 
reasonably well justified. The small disturbance theory appears 


to be the only theory available for this type of investigation 
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and (S) shows that for the momentum flows to be equal, 
k must satisfy 


k/V 4 erf(k) (9) 
Now, erf (Wa) = 0.9878, so a satisfactory approxi- 
mate solution to Eq. (9) isk = YW, and the new ve- 


locity profile to be used in place of Eq. (6) will be taken 
as 

u/u; = exp [—7(v/Y)?] (10 
This expression closely approximates Eq. (6) in the 
range —Y < y S Y, and, although it extends to in 
finity, its values for |v > JY are negligibly small. 


For convenience a new coordinate «* along the wake 
axis is introduced in Eqs. (1) and (5). It is defined by 
x* =x — 0.7¢ (11) 


and in terms of x* the wake half-width and velocity at 
the center become 

Y = 0.68V 2¢c(Cpx*/c) (12) 

ue/V = —(2.42Cp“*)/(x*/ce + 0.3) (13) 

Eqs. (10), (12), and (13) completely define the wake 

to be used in this analysis of the nonsteady effects on a 

rotor blade moving through the viscous wakes of the up- 


stream row of stator blades. Before this analysis is be- 
gun, it should be noted that the ratio of wake half- 














ROTOR 


Stator Vs 


Fic. 2. Diagram denoting stator and rotor cascades and stator 
blade wakes. 
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width to airfoil half-chord is small. For a value of ¢ 


typical of compressor blades, say 0.02, Eq. (12) shows 


that at x = 1:26, 

w/eo= 05, Y/t 0.1 
and at 4 aC, 

uv /¢ 25.. Yi 0.2 


Thus, if the airfoil is a turbomachine blade, its wak, 
is still relatively narrow compared to its chord and 
presumably, compared to the circumferential blac 
spacing, in the region of the next blade row 


APPLICATION TO TURBOMACHINI 


The viscous-wake model described in the previous 


section will now be applied to find approximately th 
nonsteady effects experienced by a blade in the roto: 
row of a turbomachine as it passes through the viscous 
wakes shed by the blades of the upstream stator row 
First, an approximate expression for the nonsteady up 
wash experienced by the rotor blade is obtained. Then 
the concepts and results of nonsteady thin-airfoil theory 
are used to find the lift and other nonsteady quantities 
associated with this upwash 

The problem may be visualized with the aid of Fig. 2 
This represents a stage of a typical turbomachine as an 
infinite cascade of two-dimensional stator blades fol 
lowed by a similar rotor row moving downward relative 
to the stator with velocity U’. The stator blades have 
chord 2c¢,, spacing d,, and stagger a, which represents 
the swirl of the incoming flow whose velocity is J’ 
The rotor blades have chord 2c, and have mean chord 
lines in the direction of the relative flow, I’, which 
makes an angle 6 with the velocity I’,. The distance 
between the centerlines of the stator and rotor rows is / 
It is assumed that each stator blade sheds a viscous 
wake of the form described in Eqs. (10), (12), and (13), 
where c, V and x* now have the subscript s, and Cp is 
to be read C,’, the stator-blade profile drag coefficient 
The wakes are symmetrical about the direction of 1’,, as 
shown. 


Velocity Induced at Rotor 


Since the rotor blade moves along a line oblique to the 
x, (or x,*) axis, it is convenient to introduce oblique 
coordinates x’, y’ as shown in Fig. 2, where the y’ axis 
is parallel to the direction of rotor motion. The relation 
between x,*, vy and x’, y’ is given by 
* , 


x,.* =x’ — y’sina, vy = vy’ cosa ld 


This simply means that the coordinate x’ is always 
measured from a vertical line (the y’ axis) in Fig. 2. 
But, as has been seen above, the wake is narrow in the 
region of interest, so that y’/x’ is small in the wake it 
self, and one may write, approximately, 


oF ae" 


y = y’ cos a, 
This approximation simplifies the analysis because 1 


retains the simple dependence of the wake half-width 
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nd centerline velocity on only one coordinate—4.e., 


from Eqs. (12) and (13), 
Y= 0.68V 2¢.(Cp'x'/e 16 
—~(2.42V Cp’)/(x'/c, + 0.3 17 


[he velocity profile from Eq. (10) is now 


u/u, = exp } —2[(y’ cos a,)/¥]*5 IS 
Eq. (18) gives the profile of the wake shed from the 
center stator blade in Fig. 2 (the zero-th blade). Each 
stator blade sheds such a wake, so that the pattern is 
repeated infinitely many times in the y’-direction. The 
profile of the wake shed from the mth stator blade may 
be found from Eq. (18) by replacing y’ by y’ — nd,. To 
obtain the total velocity field of all the wakes, one then 
sums over 2 from —© to ~. The result is the total 
perturbation velocity field due to all the wakes, wr. 
U7 [= 


COS @ : 
: ae exp | —» ( - ) (y nd 2 | (19) 
u Y } 


The assumption that the configuration of wakes be- 
hind a cascade can be described by superposition of 
isolated-airfoil wakes is, of course, approximate. As is 
often the case, the simplicity afforded by this approxi- 
mation is attended by a certain complication namely, 
Eq. (19), as it stands, does not satisfy the principle of 
continuity of the axial flow through the stage. 
of the periodic nature of the flow in the y’-direction and 


Jecause 


the fact that a shear flow with only one velocity pertur- 
bation component has been assumed, continuity re 
quires that the integral of the total perturbation velocity 
ur vanish when taken over a period in the y’-direction. 
Therefore, a constant equal to the average value of u7 
in Eq. (19) should be subtracted from the right side of 
that equation. However, the next step in the analysis is 
to expand the right side of Eq. (19) in a Fourier series, 
and the constant term in this expansion will be just the 
constant that should be subtracted. Thus, continuity 
will be satisfied if the constant term is simply dropped. 
This is the procedure to be followed. 

Eq. (19) can be transformed into a complex Fourier 
series by means of the Poisson Summation Formula,’ 
which states that for suitable f(7) 

~ fn= > | f(sje"°""* dz — (20) 
lo apply this to Eq. (19) it is necessary to evaluate the 


integral 


% 2 

exp [—A°(y’/d,; — 2)?]e~“"""" dz 
« a 

where 


K? = - cos? a,(d,/Y)? (21) 
After introduction of a new variable, § = (y’/d, — 2)K, 
and completion of the square on the exponent, the inte- 
gral is 
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, 


l ( Te ) 
> exp —Zaim — : 
Kk“ d, K: 
| =| -(¢-F) | 
“7 — te | dg 
« K oa 
VT y" 7m 
= exp | —Zein — ZZ 
K ad k? 


With this result, the application of the Poisson Summa 


tion Formula to Eq. (19) gives 


- , y’ rm 
= . 4, exp —LTIm — 
U, K m d AK 


where ® denotes the real part. The constant term in 


Eq. (23) is the average value of u;/u., which must be 
dropped to satisfy continuity, as explained above, so 


that u7/u, is to be written 


9 
uy 2VIrs 4 WH 

= —— >> exp 1 —Zrtum — (24 
U, K ” l ad 


The real part is implied from here on 

The time dependence of the y’-coordinate of a fixed 
point x, of the rotor blade can be obtained from Fig. 2. 
The coordinate x, 1s measured rearward from the center 
of the rotor blade, and the time origin, ¢ 0, is taken in 


the position shown in Fig. 2-—1.e., 


when the center 
of the rotor blade lies on the extended chord line of the 
zero-th stator blade. It is seen from the figure that 
y’ = U(x,/V, — t) (25 
LU is the downward velocity of the rotor blade, and J’, is 
the velocity relative to it. 
The velocity ur is in the direction of x 
an angle 8 with the rotor blade, so that the component 
Therefore, 


, Which makes 


giving upwash at the blade is —w, sin 8. 
the upwash at the rotor blade, v”, is, from Eqs. (24 


and (25) 


9 ” 7 
a A Vr _ 7-H x,/I 
v’ = —u, sin B . > xpi e , 
K meal K 


(26 


Here vy, is the circular frequency with which the rotor 
blades pass the stator blades, defined by v, = 2rl’/ds. 
Besides the explicit dependence on x, exhibited in the 
last exponential of Eq. (26), there is a dependence 
through u, and A, both of which depend on x’ and there- 


; 


fore on x,. The relation between x’ and x, is x’ = 


bsec a, — 0.7¢, + x,V./V,. With this relation, one can 
write 
od Gls x. ¥, " 
= sec a, + * — O.; (27) 
Cs Cc ( Cr | » 


It is this expression for x’/c, which must be used in u, 
and A and which gives them their dependence on x,. 
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Now that the nonsteady upwash 7’ has been ob- 
tained, the nonsteady effects caused by the viscous stator 
wakes at the rotor blade can be found by applying non- 
steady thin-airfoil theory to an airfoil with this upwash. 
However, the dependence of v” on x, just discussed 
makes this process complicated. Since the theory out- 
lined so far is an approximate one, another approxima- 
tion, which makes the application of nonsteady airfoil 
theory to wv simple, will be introduced. Namely, it is 
assumed that the variation of u, and K with x, does not 
have a large effect on v”’, so that in these two quantities 
v,/c, may be taken fixed at some appropriate value be- 
tween 1 and 1. In other words, in calculating the 
variation of v” along the rotor-blade chord, one as 
sumes some mean values for the quantities u, and A. 
The consequences of this assumption are investigated in 
the numerical calculations reported in the next section. 

If the fixed value of x, in u, aad A is denoted by x,» 
and the corresponding value of x’ from Eq. (27) is xo’, 
then use of Eqs. (21), (16), and (17) in Eq. (26) yields 


the following expression for v’": 


| — a? 
Dr m=1 
where 
: V, 242 VCp sin@B 0.686, 
" =4dr - : : a x 
V, Xo /Cs + 0.3 V2 cos a. 
Cains’ { 0.680, \2 Cp°xo’ 
exp | — a0 (29) 
Cs V2 cos a, / Cs 


Here o, denotes the stator row solidity 2c,/d,, and 


vo’/cs is given by Eq. (27) with x, = x,. 


Nonsteady Lift, Moment, and Wake of Rotor 


The calculation of the nonsteady aerodynamic quan- 
tities induced at a thin airfoil by a relative upwash of the 
form of v” can be easily obtained from the results of a 
paper by Sears.‘ There it is shown that if a thin airfoil 
experiences at nonsteady upwash of the form 
iis (30) 


ty 
1 => We 
Vv = Uo 


the lift and the moment about the mid-chord are given 
by 


/ 


Lit) = 2rcVuS(w)e™ 
M(t) = L(t)/2c (31) 


Here v is a constant independent of x«, which may be 
complex, |’ is the stream velocity, c is the half-chord 
(not the chord as in reference 4), x is measured rearward 
from mid-chord, and S(w) is a complex function of the 
reduced frequency w = vc/V, which was shown in 
reference 5 to be 
S(w) = jtw[Ko(iw) + Ki(iw)]} (32) 
where Ay and A, are modified Bessel functions of the 
second kind. 
The application of this isolated-airfoil result to a cas- 
cade blade is consistent with the idea of a “‘first-order”’ 
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theory, as previously mentioned. It is assumed that 
the effects of passage through viscous wakes are smal] 
compared to steady effects. Thus, in calculating nop 
steady quantities, the induced effects of other blades 0; 
the rotor cascade and of the stator cascade ar: neg 
lected. This approximation, of course, limits th 
validity of the results to cascades of relatively oy 
solidity. | 
Since each teri of v”’ is of exactly the same form as 
vin Eq. (30) and since the problem is a linear one due t 
the restriction to thin airfoils, the result of Eq. (31) can 
be applied term-by-term to v’ to find the lift and 
moment of the rotor blade. If lift and moment co 


efficients defined by 
C,(a) == [L’"(t) (1/2) pV ,?2¢ 
Cy,""(t) M(t) (1/2) pV,7(2 


are used, the lift and moment at the rotor due to v” arg 


where c, is the rotor half-chord and w, is the rotor re 
duced frequency, v,c,/1’,. 

Thus, the lift and moment coefficients are obtained as 
sums of harmonic components. In flutter and vibration 
analyses, for example, the magnitudes of the harmonies 
themselves are of primary importance, rather than their 
sum. 

Other nonsteady aerodynamic quantities associated 
with 7’, such as the nonsteady circulation and non- 
steady vortex wake shed by the blade, can be found by 
standard methods of nonsteady thin-airfoil theory (see, 


for instance, references 6 and 1). 


NUMERICAL CALCULATIONS AND CONCLUSIONS 


A few sample calculations of C,’" have been made in 
order (1) to obtain the order of magnitude of the lift on 
a rotor blade due to its passage through the viscous 
wakes shed by the stator blades, (2) to evaluate the 
effect on this lift of Cy’, the stator profile drag co- 
efficient, and (3) to investigate the sensitivity of the 
results to the choice of x,o/c,, the point on the rotor blade 
at which the constant wake width 2Y and velocity at 
the wake center 7, are evaluated. 

The geometric parameters taken are thought to be 
typical of present-day axial-flow compressors and aré 
the same as those used in the earlier calculation of non 
steady interference effects in a turbomachine stage 


They are 


Os = | 
CIC, = I 
a = 45 
B = OO 
b’/c, = 0.2 
where b’ = b — c, cos a, — c, cos (8 — a;,), the axial 


distance from stator trailing edge to rotor leading edge 
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TABLE | 
{gSOLUTE VALUES OF First Two HARMONICS OF ROTOR LIF! 
COEFFICIENTS 

. Cp* Cri" C1,’ 
«i /2 0.0] 0.029 0.020 
0.02 0.057 0.037 
1/2 0 Ol] 0.023 0.014 
0 O2 0.045 0.026 


[wo typical values were taken for Cp’—namely, Cy’ = 
(.01, 0.02. For each, the first two ternis of the series for 
C,"—.e., the first two harmonic components——-were 
obtained. 

To test the importance of the choice of x,», 
—c,/2—are used in every case. 


two val- 
ues—namely, c,/2 and 
This means that, in calculating the variation of the up- 
wash velocity over the rotor-blade chord, the quantities 
u.and AK have been evaluated at the */s-chord and the 
The results are shown 
of the first 


chord stations, respectively. 
in Table 1 where the absolute value |C,,,,”" 
two harmonics, m = 1, 2, of C,’" are tabulated. 

The convergence of the series in C,’"—-that is, the 
rapidity of decrease of successive harmonics—is not as 
rapid as might appear from the fact that the con- 
vergence factor in G,,”" [Eq. (29)] is exp (—m?r?/K?). 
Although for large m this factor will decrease rapidly, for 
small m it is extremely close to unity because 7*/K? is 
In fact, for the 
cases calculated, L107. 
the second harmonic is smaller than the first by a factor 


proportional to Cp*, which is small. 
mw’ /K? is of order Therefore, 
of only about two. 

The effect of the C,p* term in the exponent of C_,,”" is 
seen to be fairly small, because the results for the two 
different values of Cp‘ are in about the same ratio as the 
respective drag coefficients—namely, 1/2. It seems, 
therefore, that the lift produced at the rotor due to its 
passage through the viscous wakes shed by the stator 
blades is nearly proportional to the profile drag co- 
efficient of the stator blades. 


DU 


E 
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The large change in x,./C, from 1/2 to 1/2 pro 
duced changes of the order of 25 per cent in |C_,,,/" 


Thus, the approximation involved in taking constant 
values of u, and A seems to be justified for practical 
purposes. 

The results given in Table | may be compared with 
similar results given in reference | for the nonsteady 
effects at a rotor blade due to the relative motion of the 
circulation carried by the stator blades. Since typical 
compressor blades operate at steady lift coefficients of 
about unity, the values of |C_,, 
with the magnitudes of the ratio of nonsteady to steady 


may be compared 


lift on a rotor blade given in Fig. 7 of reference 1. If 
this is done, it is seen that the first harmonics (m l 

are of about the same size for C,° = 0.02, while the 
second harmonics (m = 2) of C,” are larger than those 
This is because the convergence 
It therefore 


given in reference 1. 
is slower here than was the case there. 
appears that, at the higher stator profile-drag coeffi 
cients, the viscous-wake effects on a rotor blade are of 
about equal importance with the circulation-induced 


effects discussed in reference 1. 
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The Design and Calibration of a Variable 
Mach Number Nozzle 


JOHN ROSEN* 


Division of Steam Engineering at The Royal Institute of Technology 
Stockholm, Sweden 


SUMMARY 


A method of designing a symmetric Mach Number nozzle 
with adjustable flexible walls has been developed. The 
tintious variation of the Mach Number is controlled by one 


single screw mechanism, which in a defined manner simultane- 


con- 


ously affects the coordinates and derivatives of the inlet and 


throat curvatures up to the inflection point. The flexible part 


of the curvature between the inflection point and the test sec- 


tion consists of a spring, with linear thickness distribution, 


rigidly clamped at the level of the test section (150 X 150 mm.? = 
5.9 in. X 5.9in.). Due to small curvatures and a suitable choice 
of the linear thickness distribution of the flexible part of the 


wall, a good Mach Number distribution (+0.015 or less) has been 
realized within the whole Mach Number range from 1.5 to 2.5. 


INTRODUCTION 


— EVER INCREASING DEMAND for a more advanced 
wind tunnel technique explains the keen interest 
recently shown in the development of devices for con- 
tinuous variation of the Mach Number. Some of 
these tunnels fitted with such devices can be char- 
acterized by their working principles, as for instance, 
the sliding-block type, the flexible-wall nozzles and 
the Drougge wedge. 

The investigation described below concerns a small 
flexible-wall tunnel (150 X 150 mm.” at the test sec- 
tion). To avoid the complicated mechanism generally 
characteristic of these constructions, a new method 
was found for affecting the flexible walls in a well- 
defined motion by means of one single screw mecha- 
nism. The purpose of this tunnel is to study the Mach 
Number effect in some processes which previously were 
very difficult, or even impossible, to study in a setup of 
conventional fixed-wall nozzles. Asa comparison, some 
prominent features of the variable nozzle design will 
be mentioned: 

(1) All measurements can be recorded continuously 
as a function of the Mach Number. 

(2) The time-wasting nozzle changes are eliminated. 

(3) The starting of the tunnel can be performed 
smoothly. 

(4) Instructive films can be taken over complicated 
processes, visualizing the influence of Mach Number. 

Some important investigations which could easily be 
carried out with such a device are: 

(A) The pressure recovery for*supersonic diffusors 
running at different Mach Numbers and the cause of 

Received June 1, 1954. 
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lower 


instabilities which generally occurs at mass 
flow ratios. 
(B) Some hysteresis effects associated with th 


starting conditions for supersonic diffusors. 
(C) The regulation of an air intake when diverging 
from the design Mach Number. 


SYMBOLS 


coordinates according to Fig. | 


2.9 “ 
Y, }, Z = coordinates according to Fig. 7 
M = Mach Number 
a = load 
o = stresses 
Subscripts 
t = test section value according to the Foelsch method 
1 = test section value according to the spring calculatior 


(see further Figs. 1-7) 


NOZZLE DESIGN 


A Mach Number range from 1.5 to 3.07 was selected 
as being the most interesting range for future needs 
A flexible symmetrical nozzle was chosen, the following 
conditions being imposed for the contour: 

(1) At Mach Number 2.5 the nozzle profile approxi- 
mately corresponds to a Foelsch profile. 

(2) In the whole Mach Number range only small 
curvatures are allowed. 

(3) To get a rigid system the spring has to be short 
and strong. 

(4) In order to avoid the pressure changing sign along 
the flexible part of the wall, communication is estab- 
lished between the test section and the space outside 
the springs. 

(5) The springs are rigidly clamped at the level of 
the test section. 

(6) The inlet and throat contours IOBA up to the 
inflection point (A) at the supersonic side are made 
out of one solid block. 

(7) Corrections for the boundary-layer increment 


have been estimated at 0.60°. 


Theoretical Treatment 

The Foelsch method was applied for designing a Mach 
Number nozzle the which should corre- 
spond to the elastic line of a spring, flexible downstream 
of the inflection point (A) and rigid upstream of it 
The Mach Number (1/,) at which this should happen 


contour of 


+ The Mach Number range 2.5-3.0 is to be examined at 
Svenska Flygmotor AB, Trollhattan, Sweden. 





equd 


whet! 


Be 


I 


ulatior 


‘lected 
needs 


owing 
proxi- 
small 
short 
along 
»stab- 
itside 
vel of 


o the 


made 


ment 


Mach 
orre- 
ream 
of it 
ppen 


dq at 





was chosen as 2.500, when the maximum slope angle 
¢, was to be 7.0°. (See Fig. 3. The radius (R) at 
the throat is 400 mm. and the part between B and A is 
straight The data given above are sufficient for deter 
mining the nozzle profile. 


Fs Vt — 
G oY. 14 tT, = _ = 2.6364 
} iV, 
M, = 2.500 
r\ A, = 
Q 0 oe 2.14 mT = = | = 2.0355 
r A, 
My = 2217 
75) mim. 
28.44 mm. 
sin ¢; A, al sae 
= Pa” = .«.¢0 mm. 


¥1 A 
The coordinates of the contour are given by the 


equation 


a r cos (8, — 90) — r, cos g; + 
1 cos (O, —- 90+ 8B (1 
y = rsin (0, —90) + 1 sin (0, —-0+ 6 (2) 
where 

ry = ywysin g = 475.91 mm. 

= /™% = 232.82 mm. 
v) = 7, COS g, = 470.38 mm. 

M-r(6 — @,) 


The coordinates between A and F are presented in 
Table 1. Due to the fact that yo and y, are known at 
a given Mach Number the following values can be 
determined 

k = (vy) — ¥o)/sin g, = 240.51 mm. 
Because the radius (R) at the throat is 400 mm.: 
1=« — Rig(¢,/2) = 216.04 mm. 

The .-coordinate at the throat is 
vo = [(vo + R) cos g; — R] ‘sin g; = 207.21 mm. 
The coordinates of the point B are 


Xp =X + Rsin ¢ = 255.96 mm. 
)» 
) 


Ve = vo + R(1 — cos g) = 31.43 mm. 


TABLE | 
The Nozzle Shape According to Foelsch 
gi = 7° and AM, = 2.500 

Point é x’, mm y, mm 
oO —0.835 000 28.44 
B —(). 681 48.7 31.42 
A 0.0 263 .2 57.75 
0.035 274.3 59.09 
0.091 292.0 61.07 
0.180 320.0 €3.88 
0.307 359 8 67 .24 
0.442 102.3 70.10 
0.587 448.0 72.41 
0.744 497 .6 74.04 
I 1.000 578.5 75.00 


The total length 7’ from the throat to the test section 
then is 


= X, — Xo = 3/8.52 mm 


The part IOBA of the contour is a solid block, which 
is fixed in relation to the tangent through the inflection 
point. If now the flexible part of the wall between A 
and F is considered as the elastic line of a beam with 
linear thickness distribution acted upon by a point 
load at A and with its end clamped at the level of the 
test section, the following conditions must be satisfied 
in order to have the elastic line osculating to the Foelsch 
A and F 


must be equivalent. If now the thickness distribution 


profile. The coordinates and derivatives at 


1S 
he = hol + 8 3 


the differential equation of the elastic line of the beam 
(case A) will be 


_, Ebhy' é 
° — os | 
12PL (1 + Bé)8 
Ebho® é | ’ 
= 4. ¢ 5 
12PL? 28(1 + Bé 287(1 + Bé 
Ebhy® E 
0° -_ = — : — 
[2PL* 287 (1 + Bé 
_In (1 + BE) + Co& + 6 


where 
b = 15.0 cm. 
L = $1.50 em. (the projection of the Foelsch pro 


file between A and F on the centerline 


The boundary conditions are 


6’(1) = 0 “C= —- = f(s 
60) =0 <*CG,=0 


By comparison with the data from the Foelsch profile, 
the values of 8 and P/h,* can be calculated. 

12PL? | 

Ebh,* 2 (1 + 8)? 


we a |; ba a ol 2+ 38 | 
aii) = n ( ae _ = 
kbh,* LB? 287 (1 + £B)? 


(0) = = ig 7.0° = 0.12278 7 


1.725 em. (S 


By dividing Eq. (7) by Eq. (8), P/A,' is eliminated, 
and 8 becomes 2.46. In order to limit the stresses in 


the beam (Case A) 
6PLE 


2 = - < 3,000 kp./em.2 (9 
bho?(1 + Bé)° 


(oO) me 


. GPL (1 — £8) , ; 
o(g) = =U gives £&= 


. bho?(1 + Bé)? B 


= 0.4065 


From Eqs. (7) and (9) 
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BESiGn OF A VARIABLE 
TABLI 
The Nozzle Shape According 
7 “i 

a ‘ ind 
mm ¥¢ 
) R35 {) 0) 
)_ 6S] 1S. 75 7 
( 0) 963 .2 7 
O5 hws 278.9 6 
0.1 3] 204 7 61.5! a) 
Q.9 62.0 326.2 64.81 5 
03 Q4 357.7 67.57 j 
126.0 389 2 69 80 3 
6:5 157 .5 120.7 71.55 2 
rs) $52 2 72 29 2 
0.7 183 .7 73.86 l 
Pee) 515.2 74 51 2 
0.9 46.7 74.88 U0 
1.0 578.2 75.00 () 

0.32 cm. 


242.5 kp 
Further spring data are given in Table 2. 

Phe expansion curve between A and F has also to be 
that 


creasing pressure distribution exists along the spring 


into account. Assuming a linearly de 


taxell 


between A and F, the load distribution can be written 
g(é Po (1 é 10 

ud thus the elastic line becomes 

bho! | - | 

) B+ BE+ SBE + BE] In (1 + BE 

6P L 

3(2 ByE 68 13 (6 + 98 + 36*)-In (1 4 

é (4+ 58+ @*) In (1 3)-— (44+ 38) 8)]8 (11 


In the above equation the following values have been 


used 

¢ 7.00 

D 15 em. 

h 0.52 cm 

3 ? 46 

L 31.50 em. 

P por b = 0.6-15 9 kp./cm. 
At M, 2.500 this gives a maximum deflection of 
0.21 mm. (See further Table 2.) The spaces outside 


the springs are evacuated by means of the jet in the test 
section but not to such an extent that the static pres 
sure value there is reached. Therefore the deflection 
will than indicated 
with the Mach Number 


Within the approximation made in obtaining the 


be less above and will decrease 


differential equation for the elastic line of the spring, 
we make use of another outstanding feature. This is 


the node, defined as the center of convergence of the 


tangent from point A when this is deflected by a point 
oad at ¥ 

; ‘ ; 
The general expression for the &-coordinate of the 


node is thus (see Fig. 8 


2.46, & 0.44647. Because 


MACH NUMBER NOZZLE 187 
» 
to the opring Calculation 
I; » 500 
I o 7 TT, 

é n kp./mm ky nn kp./mn 

{) 
00 0 
(Ht) ‘) t) ‘ + th) 
S4 
16 0O.0 19. 3¢ —() 6] 18.75 
10) 0.12 97 01 0.23 27.24 
52 0.18 29. S6 0.54 10.20 
61 0.?] }) 5D ) BO 0 S84 
SO 0.2] 0. 22 0.18 30.40 
OS 0] Ig 42 0.08 29 50 
17 0.12 8. 40 0.01 28.41 
Q? 0 06 97 . 30 —0. 05 27.25 
13 0). OV 26, 19 —() 09 ?26.10 
OO O.0 95.1] 0.13 4 OS 


of this property a supporting arm with a bearing at the 
node and clamped at the point A will entirely control 
the system upstream of the point A. Because the point 
A moves in a circle and the spring cannot extend, a 
compensation spring iS placed at the level of the test 
See The 

i 


#) of this spring is 0.25 mm. 


section. Figs. 5 and 6 maximum de 
flection 

These preliminary calculations first aimed at deter 
mining a theoretical nozzle profile for a certain Mach 
Number according to Foelsch, and then accommodating 
it to the elastic line of a spring. 


The next step is to check the spring for another 


single load P—e.g., that corresponding to ¢; 1.0 
With the same curvature IOBA as before, the new 
elastic line between A’ and F can be calculated. The 


new test section Mach Number (.1/,) becomes 1.S77 if 


the boundary-layer correction is excluded Upon this 
new IOBA’ profile it is possible to erect a Foelsch pro 
file A’F’. 
produces theoretically a Mach Number 0.05 lower than 
See Tables 


Finally the discrepancy between the 


This nozzle however is 10.33 em. shorter and 
the corresponding spring nozzle 3 and 
tf and Fig. 4. 
flexible-wall nozzle and the Foelsch profile has been 
checked at a 


given in Tables 5 and 6 and Figs 


value of ¢ 8.60 The results are 


3 and 4. 


DESCRIPTION 


The entire nozzle contour is made out of one steel 


spring (material: Bofors S 145) with a minimum 


PABLI 
The Nozzle Shape According to the Spring Calculation 
tt) 
Point : Yr, mm v4, mm ¢ 
oO —0.773 () 19.12 0) 
B —{() S84 97.9 50.10 $00 
0) 243.4 65.17 $00 
Q.050 259.2 3.9] 
0. 100 974.9 6 31 Ha 
0 POO) i) 4 69.19 14 
0 200 327.9 70). 7é 2 5S 
0 400 69 4 73.0 > in 
O 500 100.9 73.03 1.60 
) HOO) $32 73.80 1.140 
0. 7OO 163.9 74.35 0 83 
QO SO0 $95 44.72 0.52 
0 900 526.9 74.938 0.25 
k 1 O00 558 . 4 75. OOO 
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a; = 120 kp.’mm.? (see Fig. 5). This spring is 
clamped at the test section and at the part upstream 
of A, joined by bolts to a supporting arm pivoted 
about the node. By this means a screw mechanism 
acting upon the supporting arm, via an articulation, 
entirely defines the curve IOBA within the whole 
Mach Number range. The supporting arm also has to 
prevent the spring from cracking as could happen if 
the pressure on the inlet surface were transmitted by 
the spring. 

Due to the fact that the node is fixed and the lengths 
of the supporting arm and the spring are determined, 
a special device (see Figs. 5 and 6) is fitted in the test 
section to conserve the condition that the end of the 
spring is clamped though it changes axial position at 
different Mach Numbers. This device consists of a 
frame-spring mounted in such a way that the nozzle 
wall is pre-inclined 0.60°, approximately corresponding 
to the correction for the boundary-layer increment. 

The compensation spring is further preset to give 
higher axial tension to the nozzle spring at lower Mach 
Numbers, thus giving negligible influence on the elastic 


TABLE 4 
The Nozzle Shape According to Foelsch: IOBA’ as in Table 3 





Point é x’, mm y, mm 
O —0.773 0 19.12 
B —(). 684 27 9 50.10 
A 0 243.4 65.17 
0.044 257.4 66.10 
0.226 313.5 69.14 
0.318 343.6 70.31 
0.419 375.4 71.26 
0.525 108.8 71.82 
F’ 0.672 $55.2 72.15 
TABLE 5 
The Nozzle Shape According to the Spring Calculation 
eo, = 86 
Point é x’, mm va, mm gs 
(a) —() S68 0) 16.93 0 
B —().678 59.8 21.43 8.60 
A” 0 273.4 53.43 8.60 
0.05 289 2 8.40 
0.1 304.9 58.36 7.95 
0.2 3336.4 62.44 6.77 
ts 367.9 65.84 5.56 
0.4 599 4 68.58 $. 44 
0.5 130.9 40.75 >. 45 
0.6 162.4 72.40 2.54 
O:7Z $93.9 73.60 1.80 
0.8 525.4 74.40 1.13 
0.9 556.9 74.86 0.58 
Ik 1.8 588.4 75.00 0.0 
TABLE 6 


The Nozzle Shape According to Foelsch: TOBA” as in Table 5 


Point g x’, mm y, mm 
O —(0.868 0 16.938 
B —().678 59.8 21.45 
A” 0) 273 .4 53.73 

0.069 295.1 56.88 
0.217 341.7 62.82 
0.392 296.8 68 . 52 
0.595 460.8 73.64 
0.835 536.6 77.84 
1.122 626.7 80.58 
F” 1.35 700.1 81.23 


line of the flexible wall at higher Mach Numbers. Dy 
to the thickness distribution of the spring the elastj 
line has to be inclined 0.72° + 0.60 Ll .32 

The screw mechanism consists of a symmietrica]h 
working gear system, which transmits the control 
forces via an articulation to the supporting arm at th 
level of the throat. The height of the throat can by 
read on a circular scale placed upon one gear, so that 
after calibration, the test section Mach Number ea; 
be arbitrarily adjusted from the outside during the ru: 
Another device for Mach Number indication on phot 
graphs or films is mounted in the light beam of thy 
schlieren system. At the inlet, a thin (1 mm.) spring 
is fixed with one end upon the nozzle block and th 
other on rollers running in tracks in the floor and bot 
tom. The spring osculates to two contour blocks, on 
moving with the nozzle block and the other fixed to the 
floor and bottom. (See Figs. 5 and 9. 


SCHLIEREN SYSTEM 


For flow visualization in the test section two separat 
schlieren systems are used, which work perpendicular 
to one another asin Fig. 10. 

The principal of these two systems can be fed alter- 
natively by continuous light or by a spark (1~—2usec). 
The change of light source is performed by a _pneu- 
matically controlled mirror to allow both continuous 
observation during the run and instantaneous photog- 
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OF A VARIABLE MACH NUMBEI NOZZLE 





5 


6S ve 





ht beam passes via 
900 


raphy. From the first slit the lig 
and 
After that it 


another pneumatically 


two spherical mirrors (250 mm. diameter 


mm.) to the second slit or cut-off edge. 


The linear thickness distribution goes to a large screen or via 


hp = h (1+ 8 controlled mirror to a camera. 

In order to take 
dispersion prism can be 
slit. The schlieren edge is then replaced by a slit large 


enough to let one color through when the 


phe rt graphs, a 
after the 


coloured schlieren 








attached immediately 


pressure is 
constant in the test section. 

The vertical parallel light beam 1s 
spark in the SO 
} f = 1,000 mm. 

and fitted with pneumatically controlled cutoff edges. 
(See Fig. 10. 


produced by a 


focus of a lens of min. diameter and 


The camera side is optically identical 


RESULTS AND DISCUSSION 


The complete results of the tunnel calibrations ar¢ 
Fig. 11. A compilation of the cali 
bration results is given Table 7, Mach 
Number divergences are shown for all Mach Numbers. 
Mach Number was de 


06 A 0.7 


nozzle 
the 


shown in 


of the node 


€-coordinat¢ 


rsus the . 
in where 


The accuracy of measuring the 
termined from the scatter of test points obtained from 


several tests made under the same conditions. The 
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Summary of Nozzle Calibration Results 








( Because a fix diffuser was used during the calibration, the flow could sometimes separate before the exit of the nozz!l | 
why some values in Table 7 are excluded 
Vertical rake: O, O, Z 
V vA 
theor mim, — 62 12 22 —2 S 40) 1S WV 
1.602 (] 522 | 832 1. 534 1.536 1.535 1.534 | .53l 1.533 £0) O03 
1.833 l 5 1.759 1.757 | 75a 1.753 |. 756 l.soe 1.756 £0 O08 
2 O65 2. 029 2 O38 2 041 2 O41 2 O36 2 027 » O35 +1) OOS 
2.160 2 2.150 2.147 2.137 ?. 134 >. 141 } 
2.391 (2.296 2.299 4 2.320 2 316 2.302 2.304 2.307 0.013 
2.450 2.344 2.339 2.339 2.342 2.340 EO. O04 
2.533 2.429 2.408 2.411 2.426 2.419 
Horizontal rake: O, Y, O 
WV, } 
theor min. —60 1) —2() 0 10) 30 () VW 
2.065 2.028 >. 042 2.042 2.041 
2.261 (2.215) 2. 229 y. 230 , . 229 2.231 2.224 2.229 t) O05 
2.391 2.317 ? 320 > 390) 2 3220 2.220 2 390 +() O05 





error in the Mach Number measurement in the whole 
+ 0.007. 
The two springs are symmetrically mounted with 


region was thus less than 


an accuracy of +0.02 mm., which explains the good 
horizontal Mach Number distribution (+0.005). The 
vertical Mach Number distribution is also good, es- 
pecially for Mach Numbers lower than 2.00 (+0.008). 
This is obvious, due to the fact that the curvatures 
of the nozzle in this region are very small. Above 
Mach Number 2.00, however, the divergences become 
greater (+0.013) until the design Mach Number is 
reached ( ~2.39, theoretically 2.50). For even higher 
Mach Numbers the divergences increase once again 
and the Mach Number profile also changes character. 
(See Fig. 11.) 


CONCLUSIONS 


A method of designing a symmetrical flexible-wall 
nozzle for continuous variation of the Mach Number 


has been developed by using small wall curvatures in 
the lower Mach Number region and by applying the 
Foelsch theory for a higher Mach Number (2.50). 
The results show the variation in Mach Number to be 
Mach 


1.5 to 2.5 was obtained by one single screw mechanism 


+0.015 or less. The range of Numbers from 

The length of the nozzle from the minimum section 
to the beginning of the experimental test section was 
approximately constant and measured 4.5 test section 
heights. 
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a weak foundation upon which to build an approxi- 
mate method of solution. 
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On the Buckling of Stringer Panels 


Including Forced Crippling 


P. P. BIJLAARD? 


Bell Aircraft Corporation 


ABSTRACI A, i 

The analysis considers the effective moment of inertia of the 
ners when located on one side of the plate and the effective B 
idth of the plate at the unloaded edges It is shown that the C 
rigid beam webs supporting the unloaded edges of the plate D 
greatly increase the effective width of the plate acting with the E 
first stiflener G 
Also considered in the analysis are the effects of shear deforma ] 
tion of the stiffener web and the flexibility of the attached flange ] 


of the stiffener; each may cause an approximate 10 per cent reduc- 


in the effective moment of inertia in computing the critical 


stress for the case in which the column deflection of the stiffeners y” 


predominates 


An explicit formula is presented for the critical stress for buck 





1 the forced crippling mode, 
ny number of stiffeners 


which is valid for panels having 
A rapid method for computing the ] 


buckling stress is given. 
kK 


NOMENCLATURI 


half wave length of buckles in the direction L 
of the compressive force (Y-direction R 
= width ofa bay 1 
= distance between the rivet line and the 
middle plane of the stiffener web Y 
= spring constant in case of buckling in the y V:Z 
forced crippling mode in reference 6 a 
= distance of middle plane of skin to axis of o 
gravity of beam webs that bound a 
stringer panel Y 
= effective flange width of stiffener ; 
k = [(b*t)/(#*D)]o, buckling stress coefficient : 
Y = number of bays in a stringer panel i 
= pitch of rivets in attached stringer flange . 
= parameters in expression for ky’ in Eq. (61 v 
= radius of gyration 4 
= plate thickness . 
= displacement in the Y-direction , 
= displacement in the }-direction et 
= deflection (in the Z-direction ) 
= coordinates Subscripts 
= distance of center of gravity of effective h 
stiffener cross section to middle plane of : 
skin 
= distance of center of gravity of stiffener ; 
proper to middle plane of skin 
A, Ci, Co, C3, Cy = constants fy 
Received January 21, 1954. Revised and received November 
29, 1954 
rhis is part of a study undertaken by the author at the ” 
request of the Bell Aircraft Corporation. The author wishes to ’ 
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by Messrs. Arthur Schnitt, Alexander Krivetsky, and George S$ 
Johnston 
+ Consultant Also, Professor of Structural Engineering, uu 
Cornell I niversity ti 
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bendin 


critical 


* 


effective area of cross section of bounding 
beam webs and area of cross section of 
stiffeners proper, respectively 

nb, width of stringer panel 

spring constant or foundation modulus 

flexural rigidity of plate 

modulus of elasticity 

modulus of rigidity 

effective moment of inertia of stiffeners 

effective moment of inertia of stiffeners, 
taking account of flexibility of attached 
stiffener flange 

effective moment of inertia of stiffeners, 

taking account of shear deformation of 

stiffener web and of flexibility of at 
tached stiffener flange 

moment of inertia of stiffeners proper 

buckling stress coeflicient for forced crip 
pling applying for an infinite number of 
panels 

length of stringer panel 

C,.,63)/('D), 
4 
buckling in the forced crippling mod 

excess tension in the skin in case of up 
ward buckling of stringer panel 

shear force in stiffener 


parameter in theory of 


transvers¢ 
coordinate axes 
k.'/k 


a/b, ratio between half wave length and 


in Eq. (60 
width of a bay 

coefficient of interaction in method of split 
rigidities 

A,./bt 

A bt 

quantity that is small with respect to unity 

7T/a 

Poisson’s ratio 

normal stress 

shear stress 

coefficient of interaction in method of split 
rigidities 

Airy’s stress function 


g 


effective 


equiva 
caused 


forced 


lent 
by bending of attached stiffener flange 


crippling 


web beam 


infinite 


maxim 


um deflection, amplitude 


refers to “outer’’ as additional subscript to 6, or 


fa 


stiffene 


b,, or B,, where it refers to 


“+r, except as additional subscript to k, 


‘simply supported”’ 


or as subscript to C, where it refers to “‘shear”’ 


upper 


web 
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1, 2, 38 = refers to the number of stiffeners in the stringer 
panel 


INTRODUCTION 


i bes PROBLEM OF BUCKLING of stringer panels into 
the conventional mode has been treated in several 
publications." * * However, at the writing of the 
subject report, none of the references offered any sug- 
gestions for the evaluation of the effective moment of 
inertia for an arbitrary number of stiffeners located 
on one side of a panel. This problem was brought to 
the attention of the author whose investigation resulted 
in a report on the effective moment of inertia of stiffen- 
ers in stiffened flat plate panels with 1, 2, 3, and an in- 
Recently, somewhat simi- 
Refer- 
ence 5 assumes that the unloaded edges y = 0 and y = 
nb of the skin are not restrained in the Y-direction (Fig. 
1), so that the shear stresses, 7,,, are zero here. On the 


finite number of stiffeners.‘ 
lar calculations were published in reference 5. 


other hand, reference 4 takes into account, in a conserv- 
ative way, the attachment between the skin and the 
beam webs. It also considers the flexibility of the 
attached stringer flange, which results in a decrease of 
the effective moment of inertia of the stiffeners. 

In the present paper, formulas are derived for the ef- 
fective moment of inertia of the stiffeners, taking into 
account the compatibility of the strains in the skin and 
beam webs OA and CD (Fig. 1) and the reduction of the 
effective moment of inertia of the stiffeners by the 
shear deformation of the stiffener web and by the 
flexibility of the attached stiffener flange. General 
formulas are given for buckling into the forced crip- 
pling mode, extending reference 6, which presupposes 
panels with a large number of stiffeners. 


DETERMINATION OF THE EFFECTIVE MOMENT OF 
INERTIA OF THE STIFFENERS 


The stringer panels considered are assumed to con- 
sist of a skin supported by stiffeners of equal cross 
section, located at equal distances b, and to be bounded 
by webs (Fig. 1). The panel is loaded by compressive 
stresses in the direction of the stiffeners. 

The determination of the effective moment of inertia 
of the stiffeners which has to be taken into account in 
computing the buckling stress of the stringer panel re- 
duces to that of determining the effective width of the 
skin which can be taken into account as flange of the 
stiffeners. This effective width will be computed for 
a plate (skin) with one stiffener (Fig. 2) and for a plate 
with an infinite number of stiffeners. From these cases 
those for plates with 2 and with 3 stiffeners will be 
derived. The effective moments of inertia of the 
stiffeners so determined will be reduced to take account 
of shear deformation of the stiffener web and of the 
flexibility of the attached stiffener flange. 


Plate with One Stiffener 


In this case (Fig. 2) the number of bays = 2. Itis 
assumed that plate and stiffener buckle in half waves 
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a which vary sinusoidally in the X-direction, which 
assumption will subsequently be proved to be exact 
Hence the deflection of the stiffener is 


W, = Wem SIN AX 
where 


X= 7/a ) 


Assuming that in the half wave length as presented 
in Fig. 2 the stiffener is buckling upwards, at incipient 
buckling according to Fig. 2b excess tensile stresses ¢ 
will be introduced in the skin (plate). These stresses 
og, may vary in the J}-direction as shown in Fig. 2d 
The entire state of stress in the skin will have to 


satisfy the differential equation (reference 7, page 25 
(O'b Ox) + {2[O'b (Ox2dyv2)]} + (O'S dy") = 0 (3 
where ® is Airy’s stress function and 


oy = 0°@ Oy? 
= 0°} Ox" | 
—O°P (Oxdy) 


o Y 


Try 
Eq. (3) is satisfied by 
® = J) sinAv ) 


where } is a function of v alone and X is given by Eq. 
(2). Inserting Eq. (5) into Eq. (3), its general solution 
is (reference 7, page 45) 
@ = (C, cosh Ay + Cy sinh Ay + 

C;y cosh Ay + Cyy sinh Ay) sin Ax (6 


The effective width b,, of a bay of width 6 (Fig. 2d 
is determined by the condition that, assuming the ex 
cess tensile force 7 in it to be equally distributed over 
the width 6b,,, the strain in the Y-direction should be 


equal to the actual strain at y = O (Fig. 2a), so that 
T/(b.tE) = (or — vo,)ymo/E i 


The excess tension 7’ in a bay, being introduced into 


it by the shear stresses 7,, at v = 0, may be expressed as 


T=! [ (Try) yxo AN S 


From Eqs. (4) and (6) 7,, is proportional to cos Ax. Let 


(ticlent = A COS AX 


Then from Eqs. (7), (8), and (9) 

Dep = (A/X) [sin Av ‘(o; — vo,) y= 10 
Using Eqs. (4) and (6), this becomes 

bo = A tr2[Ci + v)A + 2C4)} (11 


From Eq. (11) the effective width b., is independent of «, 
so that the stiffener (stringer) will behave as if it had 
a flange of constant width 26,, along which the stress is 
Hence the stringer behaves as tl 


equally distributed. 
The assumed 


it has a constant moment of inertia J. 
sinusoidal distribution of the skin deflection in the X- 
direction satisfies the differential equation for its de- 


flection (reference 1, page 337) 
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which DV*w + to,(0°w/Ox?) = 0 (12) Ov Ox = —(1/G) [0° (Oxdv)] — (Ou Oy 

* €Xact i - | . 

From symmetry only half of the stiffener can be con- O*7 Ox? = —(1/G) [0° (Ox*dy) | 
sidered as support of the edge of the bay to the right O°" (Ovd% 15 
Fig. 2), so that from Eq. (r) on page 346 of reference ; ; 
: ee ; ; In here Ou, Ox may be expressed in terms of the strain 
one of the boundary conditions of the skin at y = Ois : he 
: €,, SO that, using Eqs. (4), 
pA 1/9 RI oe 4 a i Ginn 13) 
2 FI (O° Ox*) I 1(O “ O71 a Ou /Ox = 6. = ( E (ts = ve 

Se ud : [3p -2Ayy7))]' — ») ~ Qap) -2) 2) ; . 

eR 2 — v) [d*w/(Ow*dy)]} — (1/2)A,o,(0°w/dx*) (13 (1/E) [(02 /dy2) — v(O2 dx 16 

cipient . 

sses ¢ where A, is the cross section of the stringer proper. After inserting this into Eq. (15), Eq. (14) can be 

tresses Since this boundary condition contains only even deriv- written as 

ig, 2d atives with respect to x, it is satisfied by a sinusoidal dis- 

ive to tribution in the X-direction of the deflections of stringer ((1'G) [0%@/(ox*oy)] + 

e 25 and skin. The other boundary conditions (Fig. 2), (1/E) {(0®@/oy*) — v[0*®@/(Ox*dy) J} ) O (le 
w/oy = Oaty = é w= f *w/Oy"? = ‘ 71 : 

0 (2 ou ‘ - J 0 and bes 0 ind 0 ; oy U at rhe outer edges y = +0 of the plate have a zero deflec- 
¥. f — independent os the distribution “ the de- tion w. They are assumed to be connected to edge 
flection in the X “direction. _Hence the sinusoidal members, running in the X-direction, with effective 
stringer deflection assumed in Eq. (1) is indeed correct. cross sections A. (Figs. 2a and 2c). Here “effective 

The constants C; and C; in Eq. (11) have to be deter- cross section’’ means the cross section of a bar, concen- 
t . - . . . - 

mined from the boundary conditions at y = 0 and y = trated at y = bor y = —), and in the plane of the plate, 

One of the conditions at y = 0 is given by Eq. (9). which has the same rigidity against elongation by a 

The other condition is, from symmetry, that the dis- force at y = +0 and in the plane of the plate as the 

placements v in the }-direction are zero, or rather, as actual member. (The way to calculate A, is shown in 

was used in reference 8 for determining the effective the Appendix.) Hence from Fig. 2e one of the bound 

y Eq. width of the flanges of J beams, that the curvature of ary conditions at y = bis 

ution the edge vy = O in the plane of the plate is zero: 

_ : tr,,dx — A, do, = 0 
(0°7/Ox"),—0 = O (14) or 
6 ro express Eq. (14) in terms of ® one obtains from Eqs. [t7,, — A,(do, dx 0 (1S 

Id } " 

<A Further one has 

Pex- | ; . 

Try = G[(Ou Oy) + (Ov/Ox)] = —O°@/(Oxdy) fia — 19 
over yi y= : 
d be where u is the displacement in the Y-direction and G is From the four boundary conditions, Eqs. (%), (14 

the modulus of rigidity. Hence (18), and (19), using Eqs. (4) and (6), one finds 
‘ 
into 
das 
C (1 — v)Ab6, sinh 2\6 + (1 + v)A2b2(1 + 26.) + (1 — v) sinh? Ad 
b= : * - 4 
S 2r*b + (2db6, coth Ab + 1)A? sinh 2A 
C. = —[(1 — v)/2d2]A 
Let . | 20) 
C; = —[(1 + »)/2a]A 
g 
; l—yp I+ p Ci 
C, = —— coth Ab JA - coth Ab 
2°b 2r b 
10 where 
5. = A./(bt) 2] 

1 Insertion of C, and C; into Eq. (11) vields 

: : $}6.(1 + cosh 2\d) + [(sinh 2\d)/(2A5)] + 1} él 

rT %e = RT PD = - Ps : h 2» 

2(1 + v)*A7b7(26, + 1) + (5 — 2v4+ v7?) + (1 $+ v) (8 — v) (26,X0 sinh 246 + cosh 2db 

18 

if where 

d 

: dO = x/B) 

; 23 


_ B =a bt 


lfA, = 6, = 0, Eq. (22) reduces to 
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: 2 sinh 2Ad + 2rd se 
) _ f) Pe | 
Nb (1 + v) (3 v) cosh 2Ab + 2(1 + v)*A7b? + 5 ye 


is represented graphically in Fig. 3 for A, (bt 


Faq 2? 
hq. (=- 


between zero and two. 


Plate with an Infinite Number of Stiffeners 


Eqs. (1) and (6) also apply in this case (Fig. 4 
Similarly to Eq. (7) the effective width 0,; is deter- 
mined by the equation (Fig. 4 

T/(b.¢E) = (62 — vey) ymn/2/E 25) 


where from Fig. 4 the tensile force 7 in a subpanel of 


width Dis 


eX 
1 =— 2 ce ny AX (26 
Assuming 
Try)y=p/2 = —A COS AX (27) 
from Eqs. (25), (26), and (27), neglecting va, with re- 
spect to o,, 
b.;= (2A/nX) [sin Ax/(or)y=0 72] (28) 
[he reason for neglecting the influence of o, is the fol- 
lowing. If at x = Oa half wave, buckling upwards, is 
idjacent to a half wave that buckles downwards, from 
consideration of antisymmetry at x = O the displace- 
ment v in the }-direction will be zero there. In other 
cases v at x = O will practically remain zero owing to 
restraint by the connection of the skin to a rib. With 
in infinite number of panels of width } a line y = 6/2 


will remain straight, so that then everywhere at y = 
)2 the displacement v is zero. This means that the 
contraction in the )-direction is restrained, so that for 
large a/b ratios the apparent modulus of elasticity for 
extension in the -direction is increased from / to k + 
) as a consequence of the stresses o, = vo, 
the Y-direction. Since the effective width 

actually is used for cases where v is finite and the 
contraction in the Y-direction will be fully prevented 
only if nb is large with respect to a, lest b,; is overesti- 


l-*% gen- 


erated in 


mated, the influence of a, was neglected in Eq. (28). 
Only for small 8 values, where o, is opposite in sign to 
7,, the influence of o, decreases b,;, but the difference 
is negligible here. 

From symmetry in the present case Eq. (6) reduces 
to 


(29) 


® = (C, cosh Ay + Cyy sinh Ay) sin Ax 


Using Eqs. (4) and (29), Eq. (28) becomes 


= 2A /(\*4/ CA cosh (AB/2) + C;[2 cosh (Ab/2) + 
(Ab, 2) sinh (Ab/2)]}) (80) 
The boundary conditions at y = 6/2 are given by Eq. 
7) and by the condition that the edge y = 5/2 re- 


mains straight, as discussed above. Analogous to Eq. 


17) this condition is 


((1/G) [0%@ (Ox20y)] + (LE) { (0% dy 
v[O*@ ‘(Ox*dy) }} ) O (51 

Using Eqs. (4) and (29) the conditions (27) and (31 

vield 

C {{(1 — v) sinh (Ab 2 1+ v) (Ab/2 x | 


cosh (Ab 
[2\ sinh (Ab 


2)]/ [2A sinh? (Ab 


2))\A | 


Insertion of C; and C; into Eq. (30) gives 


C, = (1 + pv 


({4(cosh Nb — 1)] 


b.; = V{A((B + v) X 
')h 33 


sinh Ab — (1 + v)Ad]} 


Eq. (33) is presented graphically in 
10] 


where Ab = 7 8B. 


Fig. 5. It was also given as Eq. in reference 6 


Plate with an Arbitrary Number of Stiffeners 


For a plate with an arbitrary number of stiffeners 
the boundary conditions are less simple. For example, 
for the plate with two stiffeners in Fig. 6a, assuming it 
}-direction (Fig. 
6b), the stresses or displacements at the edge y 5/2 
rhis edge 


to buckle into one half wave in the 


of the center bay are not directly known 
neither remains straight, nor are the stresses a, equal 
to zero. It may be shown, however, that the stresses 
g, at the edges do not appreciably influence the effective 
width &,. 
an infinite number of stiffeners, dealt with in the pre 
ceding section, where the edge y b 2 
a single bay reinforced by 
+ b/2 


This may be done by comparing the case of 


remained 
straight, with the case of 
stiffeners at both edges, where at the edges y 
the stresses o, vanish (Fig. 6d). Also for this case 
Eqs. (29) and (30) may be used. The boundary con- 
ditions at y = } 2 are given by Eq. (27) and by 


(g)enuw2 = 0 (34 


y/y=b/e 


Using Eqs. (4) and (29), Eqs. (27) and (34) yield 


C, = —$[b sinh (Ab/2)]/[A(sinh AB + AD) IIA | 
C, = {[2 cosh (Ab/2)]/[A(sinh AB + Ad) I} A f 


in reference 8, 


and D 


referring to box and U beams under transverse loading, 


These constants are similar to A, 


where the same boundary conditions [Eqs. (27) and 
(34)] apply, if in there 2b, A,, and a, are replaced 
by 6, A, and X, respectively. Insertion of C; and C, into 
Ea. (30) gives the effective width 

(56 


b., = {(sinh Ab + Ad) [Ab(cosh AB + 1) ] fd 


7/8. Eq. (36) is given in Fig. 5 by the 


Comparison with the solid curve for 


where Ab = 
dashed curve. 
b,; shows indeed that the influence upon the effective 
width of the stresses o, at the edges y = +0 2 is ex- 
tremely little. This is partly due to the neglect of 
a, in Eq. (25), of which the reason was explained in the 


preceding section. Hence, in the case of an arbitrary 
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number of panels and assuming buckling into one half 
wave in the }’-direction, it will be appropriate to assume 
the effective width of the internal bays equal to b,; from 
Eq. (33) and that of the external (outer) bays equal 


to b,, from Eq. (22). 


Effective Width of Stiffener Flanges 

(A) Plate with One Stiffener (n = 2)—From the 
above section of the same title it is evident that the 
effective width of the stiffener flange is 


fa = 2b (37) 


eo 


where b,, is given by Eq. (22) and Fig. 3. The way to 
determine A, is shown in the Appendix. 

(B) Plate with Two Stiffeners (n = 3) 
with one half wave in the }-direction (Fig. 6b) it fol- 
lows from symmetry that the effective width of the 
flange of each stiffener is equal to the effective width 
b.o of an outer bay plus half the effective width 6,; of 


For buckling 


the center panel, so that 


fea = Deg + (1/2) be (38) 


e2 


where b,, is given by Eq. (22) and Fig. 3 and 6,; is given 
by Eq. (33) and by the solid curve in Fig. 5. From 
reference 3 buckling into two half waves in the J- 


direction is practically not governing. Since reference 


1955 


SCIENCES UL, 


3 deals with stiffeners with given flexural rigidities, this 
applies if for both modes (one and two half waves in the 
\’-direction) the effective width of the stiffener flange 
would be the same. Actually in the case of two half 
waves (Fig. 6c) from antisymmetry with respect to the 
line y = O the latter will not change in length, so that 
for each half wave it will act as a member with an in 
finite cross section A,. Therefore, in the case of two 
half waves the effective width of the flange of each 
so that buck 
ling into two half waves in the }-direction is still Jess 


stiffener is larger than that from Eq. (38), 


likely to be governing. 
(C) Plate with Three Stiffeners (n = 4) 
reasons as in the former case, also here only buckling 


For the same 
with one half wave in the Y-direction is considered 
For small ratios 8 = a/b the excess normal stresses g, in 
the stiffener flanges (the plate) will not influence each 
other much, so that the effective flange width of th 


outer stiffeners is about 
feo = Deg + (1/2) Bi: (39 
while that of the center stiffener is 
Fae = Des (40 


From Figs. 3 and 5 it may be seen that for small ratios 
8 the flange widths /;, and f;, will not differ much, so 
that also the effective moments of inertia of the stiffen 
ers are practically equal. In that case, from reference 
9, the deflection varies sinusoidally in the }-direction, 
so that the deflections of the outer stiffeners are (1 /+/2 
times the deflection of the center one. Since the inter- 
nal work stored in the stiffeners is proportional to the 
square of their deflections, the influence of the moment 
of inertia of the outer stiffeners on the total internal 
work to be done at buckling is about (1/72)? = | 
times that of the moment of inertia of the center stiff- 
ener. Hence, the moment of inertia of the two outer 
stiffeners together has about the same influence as 
that of the center stiffener alone. Therefore, the 
effective flange width of all stiffeners may be assumed 
to be the average of f;, and f;, from Eqs. (39) and (40), 
so that 


f.3 = (1/2)beo + (3/4)0; (41 


J €o 0 


where b,, and b,; are given in Figs. 3 and 5. UH, for 
example, A, = 0, for large a/b ratios, from Figs. 5 and 
5, boo = be; = b. From Eq. (41) this leads to a total 
flange width for the three stiffeners of 3.75 6, which 1s 
somewhat smaller than the actual total effective width 
of 40. 

(D) Plate with an Infinite Number of Stiffeners 
From the above section with the same title the effec- 
tive flange width of all stiffeners is 

f.; = de, (42 


where b,; is given by Eq. (30) and Fig. 5. 


EFFECTIVE MOMENT OF INERTIA OF A STIFFENER 


The flexural rigidity of the skin (plate) is fully 
accounted for in the action of the skin as a plate. 
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BULCELING OF 


[herefore, in computing the effective moment of in- 
ertia of a stiffener the rigidity of the plate against ex 
tension alone is taken into account. If A,, /;, and z 
we the cross section and proper moment of inertia of 
the stiffener and the distance of its center of gravity 
irom the middle plane of the plate, respectively, the 
enter of gravity of the effective stiffener cross section 
is determined by the distance z (Fig. 7) 


g = ({(Az,)/(A, 


+ fic) (43) 
where follows from Eqs. (37), (38), (41), or (42). 
Hence, the effective moment of inertia is 


T=1,+A,(s, — 3s)? + f.tz? (44 


EFFECTIVE MOMENT OF INERTIA 
SHEAR DEFORMATION OF THE 
ATTACHED 


REDUCTION OF THE 
FF A STIFFENER BY 
STIFFENER WEB AND FLEXIBILITY OF THE 
STIFFENER FLANGE 
In the case of deformation by bending alone the re- 
straining transverse force R acting on an element of 
the stiffener of length dx is (Fig. 8) 
R, = —dV = EI(d‘*w,/dx*) dx (45) 
or, using Eq. (1) 


R, = METw, dx (46) 


support of a stiffener, without any 


by an elastic foundation with a modu- 


In the case of 
flexural rigidity, 
lus C,, the resisting transverse force per length dx would 
be R, = C,w, dx, so that, 

16), the flexural rigidity of the stiffener is equivalent 


by comparison with Eq. 


toa foundation modulus 
G = ME (47) 


The shear deformation of the plate was fully taken 
into account in calculating its effective width. The 
shear deformation of the stiffener web causes an angle 


rotation. 


dx =-1.11V/4,6)1 (48) 


yy = dw, 


Where |’ is the shear force (Fig. 8), A, is the cross sec- 


tion of the stiffener web, and the factor 1.1 takes ac- 
count of the variation of the shear stresses in the web in 
the vertical direction. From Eq. (48) and Fig. 8 the 
restraining force acting on an element of length dx, if 
only shear deformation of the web occurs, is 


R, = —dV -(A,G/1.1) (d?w,/dx?) dx = 
\2(A,G/1.1)w, dx (49 


Denoting 


R, would be C.w,dx, so that 


ot shear deformation by C,, 
from Eq. (49 


C, = A*(AgG/1.1) (50 


Hence, the plate may be assumed to be supported by 


two spring systems, built in series, 
Stants C, and C 


flarnt; e e ° 
flection by a unit transverse load per unit length 1s 


Therefore, 


the imaginary foundation modulus in the case 


with spring con- 
the total stiffener de- 


sT 
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w= (1/G) + (1/C,) = (GQ + C.)/(GC,) = 1/G, 
so that the equivalent spring constant 1s 
C,, = (C,C,)/(CG, + C,) (51 


, with G = E-/(2 + 2p 
shows that the 


Using Eqs. (47), (50), and (51 
and v = 0.3, comparison with Eq. (47 
effective moment of inertia of the stiffener is reduced by 


shear deformation of the web to 


I’ = [1/(LAD IG, = [(@WA,)/(287 + a?A,) | (52 


The reduction of J by shear may be as large as about 
(52) it will be 


If in Fig. 9 curve k, 


10 per cent. From Eq. smaller, the 
larger the half wave length a 1s. 


in the equation 


Scr = Rer[(w*D) / (b7t 


represents k, 
) } (53 


for the critical stress, where D is the flexural rigidity of 
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the skin, without taking account of the reduction of 
I due to shear, curve k,,’ is that taking account of J’ 
from Eq. (52). 

If the skin were supported by a stiffener that is in- 
finitely rigid against bending and shear, but with a thin 
flexible attached flange, its buckling stress coefficient 
would be given by curve ky in Fig. 9. The actual 
buckling stress, taking account of all deformations, will 
, resulting from interaction of the 


" 


be given by curve k,,’ 
curves k,and k,,’. This interaction curve has a similar 
form as the curve for o,,, resulting from interaction of 
the curves o; and o», for plate and column buckling, 
respectively, in Fig. | of reference 10. Curve k,,” in 
the present Fig. 9 has a minimum k,, at a small a/6 ratio. 
This is the buckling stress for forced crippling con- 
sidered in reference 6. The right part of k,,”, for larger 
ab ratios, is slightly below the curve for k,,’.. This 
reduction of the long wave buckling stress by flexi- 
bility of the attached flange has to be considered here. 

The equivalent spring constant C; as a consequence 
of the flexibility of the flange is calculated as shown 
in reference 6 for c,,, but without taking account of the 
deformation of the stiffener web in its own plane (Case 
c), which influence may be assumed to be taken into 
account in C, and C,. 
can occur in the entire structure, it will be safer to 
assume Cy, equal to c,, referring to the upward half 
waves in reference 6. Hence, the total deflection by a 
unit transverse load per unit length, taking account of 


C,, Cy and C,, 1s 
w= (1/C,) + (1/C,) + (1/Cy) = 1/Cisy 


If only one or a few half waves 


so that the equivalent spring constant is 
Coss - (C,C,C;) (C,C, = CAL, i C,C,) (54) 


Using Eqs. (47), (50), and (54) in comparison with Eq. 
(47) leads to a reduced effective moment of inertia of 
the stiffener 


I” = (1/EM)Cisy = { (4A Cy) /[An(1OOET + 
atC,) + 28a7JC,|}I (55) 


Also, the reduction of J by flexibility of the attached 
flange may amount to about 10 per cent. 

In order to avoid complicated calculations, C; may 
be assumed to be equal to c,,, as calculated for ob- 
taining the graphs in Figs. 23 and 24 of reference 6, 
where panels with a large number of bays are con- 
sidered and 

R = (€eqb*)/(r*D) (56) 

so that 
Cy = Ceq = (41D/b*)R (57) 
The pertinent values of R are given in Fig. 10 for vari- 
ous ratios p/by, (bg/t,)/(b/t) and (b,/t,)/(b/t), from 
which C; can be calculated by Eq. (57). The distances 
b, t, ba, by, and t, are indicated in Fig. A3. 
is the pitch of the rivets in the attached stringer flange. 
An example of such a calculation is given in the Appen- 
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BUCKLING IN THE FORCED CRIPPLING Mop 


In reference 9 a direct formula was derived for the 
buckling stress of stringer panels with any number of 
stiffeners by the method of split rigidities, using the 
more refined method of comparing the work done by 
the deflecting and restraining forces, in a way similar to 
that done in references 11 and 12. This formula yields 
results that are in practically complete accordance with 
those from the graphs of reference 3. 

From Eq. (47) an elastic support with a foundation 
modulus ¢,, is equivalent to support by a stiffener 
with a flexural rigidity. 


CED) 55 = Gea X* (58 


d 


and a cross section A, = 0. Hence, the above men- 
tioned formula for the buckling stress of stringer panels 
also applies for forced crippling if the flexural rigidity 
of the stiffeners is expressed by Eq. (58) and A, = 0 
This results in the following buckling stress coefficient 
k., in Eq. (53) [Eq. (93) of reference 9]. 


ke = ko + {(1 + @)/[201 — y6)]} X 
[1 — (1 — {[401 — yO)a]/(1 + a@)?})""] ky (59 


In Eq. (59) 
a oe: ky (60 
whereas, from Eqs. (34) and (92) of reference 9, 


ky = [(B/n?) + (1/8)? 
k, = Re? (61 
ko’ = s — (2/n?) + [q -— (1 n') a2) 


9 


where R is given by Eq. (56) and s, g, and 7@ depend on 
the number of bays, , and are given in Table 1. 
If in Eq. (59) the form within the square brackets is 
written conservatively as 
1—(1—6” =1 — [1 — (1/2)e] = (1/2)« 
it reduces to 


ae = Ro + | (RiRe’) (Ry + ko’) | (62 


With the present notations this formula is equivalent 
to Eq. (57) of reference 6, which, as was stated in its 
derivation on pages 25-29 of Part I of reference 6, 1s 
indeed conservative. For the case calculated on page 
34 of Part II of reference 6, it was found that R = 3.07 
and 8 = 0.86. Further n is infinite, so that from Table 
1 of the present paper s = 2.5, q = 5, and 76 = 0.707. 
Hence, from Eq. (61) 


ko = 1/87 = 1/0/86? = 1.35 
by = KB = 3107 0:86)" = 2.27 
ko’ = s + gB? = 2.5 + 5(0.86)? = 6.2 


so that from Eq. (59), k,, = 3.12. This is only 2.5 per 
TABLE 1 
n § 7 79 
2 2.27 2.45 0.570 
3 2.36 3.58 0.616 
4 2.37 4.20 0.639 
co 2.5 5 0.707 
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BUCKLING OF 


cent more than k,, = 3.04, found on page 34 of Part I 
of reference 6 from Eq. (57) of that reference. 

An error of not more than | per cent is made if it is 
sumed that for any number of bays the minimum 
buckling stress occurs for the same half wave length as 
for an infinite number of bays. Hence, values R in the 
present Eq. (57) can be read from Fig. 10, while the 
half wave length a = 8b follows from the equation for 
3in the second column on page 27 of Part II of refer- 
ence 0 


8 = [(K.,/2) — 0.0134K,,24]7~ “ (63) 


where A,, is the buckling coefficient for an infinite num- 
ber of bays and can be read from Figs. 23 and 24 of 
reference 6 as a function of p/bg, (b,/t,)/(b/t) and 
h/t.) (b/t).* As stated in reference 6, the graphs 
in these figures were calculated for ratios b/t = 25 and 
t = 0.5 but can be used sufficiently accurately for 
1) > b/t > 20 and 0.6 > ?¢,/t > 0.4. 
number of bays n, the buckling coefficient k,, is given 


Hence, for any 


by Eq. (59), where @ is given by Eq. (60), ko, ki, and 
b.’ by Eq. (61), and s, g, and y@ are given in Table 1. 
An example of calculation is given in the Appendix. 


PRACTICAL DETERMINATION OF BUCKLING STRESS OF A 
STRINGER PANEL 


To determine the buckling stress for long wave buck- 
ling (in contrast to forced crippling) of a stringer 
panel, one has first to estimate the half wave length 
for which it becomes a minimum. After estimating 
the effective flange width of the stiffeners and thus 
their EJ, this can be done by examining, in the graphs 
of reference 3, the curves for constant //J/dD (in the 
present notations /://bD) situated to the left of the 
dashed curves, where the abscissas represent the ratio 
1b (for curves to the right of the dashed curves the 
abscissas represent 2a /b, 3a/b, etc., since in the present 
paper a denotes the half wave length). The actual 
minimum of o,, will occur with a ratio a/b which is a 
little smaller than that for which the curves in reference 
3 show a minimum, since the effective moment of 
inertia increases with the ratio a/b. This is due to 
the larger effective flange width for larger ratios a/b 
as well as to the smaller reduction of the effective mo- 
ment of inertia of the stiffeners by shear deformation 
of the stiffener web and by flexibility of the attached 
flange for larger a/b ratios. 

For the ratio a/b so chosen, one calculates the effective 
flange width of the stiffeners from Eqs. (37), (38), (41), 
or (42), where b,, and b,; are found from Figs. 3 and 5. 
The area A, in Fig. 3 is calculated from Eq. (A3) as 
shown in the Appendix. Then the effective moment 
of inertia J of the stiffeners is computed from Eq. 
+4). This moment of inertia is reduced to 7” from 
Eq. (55) to take account of shear deformation of the 
stiffener web and flexibility of the attached flange. 

“As follows from Eq. (180) of Part II, in the conservative 


equation for ¢.. on page 2 of Part I of reference 6, the middle 
solidus should be omitted, so that ¢,, = 0.60 E (t,/b,) (t/b). 
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Finally k,, is scaled from the graphs of reference 3, or 
computed from Eq. (32) of reference 9, in both cases 
using /” from Eq. (55) instead of J. If one is not sure 
that the chosen ratio a/b leads to the smallest possible 
critical stress, one has to consider another ratio a/b. 
In most cases this will not be necessary, since only 
half wave lengths that are an integer part of the 
stringer length have to be considered. 

If k,, is higher than k, = 4 for buckling of the bays as 
simply supported plates, indicated in Fig. 9, the latter 
value is governing. 

Finally, one has to investigate whether or not a 
smaller buckling coefficient k,. results from buckling 
into the forced crippling mode (Fig. 9). This compu 
tation has been worked out in reference 6 for a plate 
with an infinite number of stiffeners. In the preceding 
section this has been completed in order to use it also 
for cases with 1, 2, and 3 stiffeners. A numerical 
example of this calculation is given in the Appendix. 
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x Y def 
A, = I,/(1,? 2? 
O 3 5 v > > P=| : Io/(r0* + *) AB the 
yZ bD 
ty where 
A D ™ 0.1 
B Ki/2 Bes 7” = f,/A, (A4 and 
FIG. Al If more girders are present (Fig. A2), at buckling, fron a 
antisymmetry in C, the superimposed elongation jn will 
the Y-direction of the skin is zero, so that A, is infinit 
——— However, in O it is given by Eq. (A3), so that also 
. Tv P= =T v here A, should be assumed to be given by Eq. (A3 (B) 
throughout. W 
A D 
+ 
FIG A2 EXAMPLE OF THE NUMERICAL COMPUTATION B= 
A) Esti ° ' the Opti Half Wav ther 

° Bijlaard, P. P., and Fisher, G. P., Interaction of Column and (A) Estimation of the Optimum Half ave Length cur 

— Buckling in Compression Members, NACA TN 2640, March, The data of the stringer panel considered are given the 
952 : ~: 9 Q9-: I , 
: 19, Ad. > Tey = 72. Ls = ).: ‘ = = 

1! Bijlaard, P. P., Analysis of the Elastic and Plastic Stability = Pag i memer, 4 ts ™ Cane te, 4 : wid 
of Sandwich Plates by the Method of Split Rigidities—IIT, Journal 0.102 in., by = 0.80 1n., ts = 0.040 in., by = 0.73 in., b = 1 
of the Aeronautical Sciences, Vol. 19, No. 7, pp. 502-508, July, = ©1355 im., p = 0.75 1n., ad, = 5/32 in. The panel web 
1952. length is 12.5 in. The number of bays n = 2. The tanec 

12 Bijlaard, P. P., and Fisher, G. P., Column Strength of H-Sec cross section of the stringer is found to be A, 0.0754 pres 
tions and Square Tubes in Postbuckling Range of Component sq.in. Its proper moment of inertia J. = 0.008184 in.‘ sa 
Plates, NACA TN 2994, August, 1953. j q: Mite ie — = piggiet Lane nes ; ~ mm 

Estimating the effective flange width /, of the stringer as (Fig 
; , - oe Si » flange section is f,t = 5 2) = 
Appendix Te ; 2b > in., the flange section is fel (0.102 | 1.64 
0.51 sq.in. The distance z, from the center of gravity 
COMPUTATION OF THE EFFECTIVE CROSS SECTION A, st wpe ee bas ae eee pene yn oe pate (Fig. 
(n tan Wane is 0.391 in. With these values one finds from Eq. ( oy tt 
z = [0.0754(0.391)]/[0.0754 + 0.51] = 0.0502 in : 

If panel OC buckles (Fig. Al, a), the webs OA and Then from Eq. (44) oe 
CD act as edge members. The effective edge member 
through O and running in the X-direction, should I = 0.008184 + 0.0754(0.391 — 0.0502)* + (C) | 
have the same rigidity against elongation in the Y- 0.51(0.0502)? = 0.018214 in. W 
direction, by a force through O in the X-direction, as aie ic 

. . . ° ° 4 2 O 

the girder OA. Point O is located in the middle plane Further a 
e ° _ ° ° ° ° x \ 
of the skin. To determine the width of the tension D = Ef®/[12(1 — v?)] = 0.0000972E —_ 

- : . ‘ - . sit ara | ere 
skin AD that can be considered to belong to the girders, | @a3 
it is observed that it has the same boundary conditions so that bD) 
as the plate considered in Fig. 6d, so that its effective es 

° . . " 92 “¢ ° | Ave 
width B,, is given by Eq. (36), if there b is replaced by B=5 i | ae 
B (Fig. Al, a). Hence, the width of the tension skin a es 

, ; tsalgh . bs=0.73 ! hc 
belonging to the girders OA and CD is (1/2)B,, (Fig. b=25 Cis 
Al, b). B,,/B can be obtained from the dotted curve 075 > f*b,= 0.355 with 
in Fig. 5, if in Fig. 5 the notations a/b and b,,/b are re- ' ‘ ad le 0.375 245 
placed by a/B and B,,/B. The axis of gravity of the : [t=Q102-y —}, bt, 
cross section so obtained can be determined. Let it be ' j A seal 

3/4x3/4x3/16 3/32R ap finds 
located at a distance e below the middle plane of the L t,=040 57 
oi 


compression skin (Fig. Al, b). The (fictitious) stress _— 


in such a section at z = O (Fig. Al, a), caused by a unit 0.125 WEB by 80 | paie 
force P along the X-axis (Fig. Al, c), is noo ase the < 
o = (1/A,) + (Me/I,) (Al) 62.8 


where A, and /, are the cross section and the moment of 








inertia with respect to },Y; (Fig. Al, b) of the girder. - “—_ (D) 
The bending moment A/ is equal to Pe = e. Inserting t=0102 W. 
this into Eq. (A1) and requiring that the elongation, and ALL RIVETS ARE 5/32" AT p= 3/4" SPACING. 62.8 
thus the stress, atz = 0, is equal to that of a member DISTANCE BETWEEN RIBS IS 12.5 INCHES. the g 
with effective cross section A, through O, one obtains — 
ALL MATERIAL IS 75 ST(SKIN & STRINGERS). r 
reler 


o = (1/A,) + (e?/J,) = 1/A. (A2) Fic. A3 





(A4 


zy, fron 
tion in 
nfinite 
at also 


|. (A3 


panel 
The 
).O754 
4 in. 


ger as 


1 in.! 


S). 





FI (bD) = 0.018214E [2.5(0.0000972)] = 75 


Assuming a reduction of the moment of inertia by shear 
deformation of the stringer web and by flexibility of 


the attached flange of about 15 per cent, one has J] + 


bD 64. Further, 6 = A,/(6t) = 0.0754/(2.5 X 
0.102 0.296. Since, with a panel length of 12.5 in. 
and 6 = 2.5 in., the ratio 8 = a/b can be only 12.5/2.5 


= 5, 5/2, 5/3, ete., it follows from the graphs on pages 


12 and 13 of reference 3 that the smallest buckling stress 


will occur with a = 12.51n. and @ = 5. 


(B) Calculation of the Effective Cross Section A, of the 
Webs 


The width of the lower (tension) skin (Fig. A3) is 
B=5in. The effective width is given by Eq. (36) if in 
there bis replaced by B. It is also given by the dashed 
curve in Fig.5. Witha/B = 12.5/5 = 2.5 from Fig. 5, 
the effective width B,, = 0.78B. Hence, the effective 
width of the tension skin belonging to a web is 0.39 (5) 
= 1.95in. It is found that the axis of gravity of the 
web section so obtained (Fig. A4) is situated at a dis- 
tance e = 4.543 in. below the middle plane of the com- 
pression skin. Its cross section is Ag = 1.645 in. and 
its moment of inertia /, with respect to the axis },}; 
(Fig. A4) is 14.472 in.4.. Hence, 7,2 = J,/A, = 14.472 + 
1.645 = > 


8.8 sq.in. Then from Eq. (A3) 


A, = 14.472/(S.8 + 4.5437) = 0.4925 sq.in. 


so that A,/(bt) = 0.4925/0.255 = 1.93 and a/b = B = 
5. With these values from Fig. 3, b,, = 0. 


(C) Computation of the Reduced EI of the Stiffener 


With b,, = 6 the effective width of the stiffener flange 
from Eq. (37) is fu = 2 6.. = 2b = 5in. This is the 
same effective width f, as was assumed above (if this 
were not so, it would have to be checked whether with 
the newly found value /,; and the resulting ratio EJ + 
bD) the same optimum half wave length a is found as 
above, which will generally be the case). Hence, 
I = 0.018214 in.4as above. Further, in Eq. (55) A, = 
bt; = 0.8 & (0.04) = 0.032 sq.in., a = 
C; is given by Eq. (57), where R is found from Fig. 10 
with p/b, = 0.75/0.355 = 2.11, b/t = 2.5/0.102 = 
24.5, bg/t; = 0.355/0.040 = 8.87, (b,/t;)/(b/t) = 0.362, 
by. ts = 0.80/0.040 = 20, (b,./t,)/(b/t) = 0.817. 
finds from Fig. 10 that R = 1.30, so that from Eq. 
7) Cy = (0.0000972F 4 2.5*) & 1.3 = 0.00079F. 
Insertion of a, Ay, 7, and C; into Eq. (55) gives a re- 
duced moment of inertia /” = 0.8377. Hence, from 
the above, with EJ /bD = 75, EI" /bD = 0.837 (75) = 


62.8 


12:5 in., and 


One 


(D) Computation of the Buckling Coefficient k., 


that is, EJ’/bD = 
(bt) = 0.296, B = a/b = 5, n = 2 
the graphs on pages 12 and 13 of reference 3, one finds 
ke, = 3.045 


reference Q. 


With the values found above 


es a 
02.8,6 = A, from 


The same value is found from Eq. (38) of 
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7.625x 0125 
3/4x 3/44 x 3/16 
1.95x 0.102 
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(E) Buckling in the Forced Crippling Mode (Simplified 

Method) 

For the ratios p/b, = 2.11, (bg/ts) /(b/t 0.362, and 
(b,,./t,)/(b/t) = 0.817 as calculated above, as found 
above from Fig. 10, R = 1.30. 
of reference 6, with the same ratios, k,, for an infinite 


From Figs. 23 and 24 


number of bays is 2.08. Hence, from Eq. (63), 


8 = [1.04 — 0.0134(2.08)?"] = 1.02 
From Eqs. (61) and Table 1, with” = 2, 
ky = [(1.02/4) + (1/1.02)}? = 1.528 
k, = 1.3(1.02)? = 1.35 
ko’ = 2.27 — (2/4) + [2.45 — (1/16)]1.02? = 4.25 


9 = 0.570 


L325 = 3.15. 


From Eq. (60) @ = 4.25 Hence, from 

Eq. (59), 

ko, = 1.528 + (C+ 3.15)/[2(1 — 0.57) ]} X 
fi — (1 — $[4C) — 0.57)3.15)/(1 + 3. 


1.35 = 2.67 


(F) More Accurate Calculation of the Buckling Stress for 

Forced Crippling 

An accurate calculation according to reference 6 leads 
to R = 1.505, 1.475, and 1.44 for B = 0.98, 
1.005, and 1.08, respectively. With these values one 
finds from Eq. (59) k,, = 2.78, 2.77, and 2.78, respec- 
tively, so that k,, = 2.77 for the optimum value 6 = 
1.005. 
above is 4 per cent conservative. 


a/b = 


Hence, the approximate value k,, = 2.67 found 


(G) Minimum Buckling Stress 


The minimum buckling coefficient is that found for 
forced crippling. Using the approximate value k,, = 
2.67, which is calculated easily, the elastic buckling 
stress is 
Ser = er [(w2D) /(b*t)] = 2.67 1(0.0000972E x?) 


[2.57(0.0102) ]} = 0.00401 = 43,000 Ibs. per sq.in 


Actually for the material used (75S-T) this stress is 


With b/(t Vk.,) = 2.5 
1S of reference 6, 


within the plastic range. 
(0.102 79267) = 15.0, from Fig. 
ocr = 36,000 lbs. per sq.in. 
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The Effect of the Temperature Dependence of To apply Eq. (1) to turbulence measurement, the usual assumy 


King’s ‘‘Constant’’ A on the Hot-Wire 
Sensitivity Coefficient 


J. R. Ruetenik 
Senior Engineer, Aeroelastic and Structures Research Laboratory, 
Department of Aeronautical Engineering, MIT 


January 12, 1955 
N A RECENT SERIES! of turbulence-level measurements (1’) 
using a hot-wire anemometer the constant-current? and con 

stant-resistance® methods of calibration gave sensitivity coeffi- 

cients differing by 25 per cent over most of the velocity range (3 

to 16 ft 

essentially constant current 


per sec.). In both cases the wires were employed at 
The discrepancy was traced to the dependence of the coefficient 

A in King’s Equation 
I?R = |A(T) + BV U| (T — T.) (1 
on the wire temperature, 7; /?R is the rate of heat loss, l’ the 
air speed, and 7, the unheated wire temperature, taken as con 
stant. This A(T) effect, as hereafter called, had been identified 
and experimentally evaluated by King but has been omitted in 


a number of interviewing investigations.* % > 6 
King empirically obtained in air 
AzA,{1+c¢T — T,)] tz 
where A, is a constant, and c = 0.00114 at 7, = 17°C To con 


firm Eq. (2) a 0.00015-in. diameter by 0.040-in. long platinum 


a, = 0.4, 0.7, and 


wire was calibrated at three ‘‘overheats’’: 
(R — R.)/R,. 


per sec., and two wire orientations were tested: 


1.0, where ay’ The air speed was 15 to 35 ft 
normal and 
inclined 40° to the stream. The comparison to Eq. (2) is shown 
in Fig. 1; B was assumed constant (justified by King and con 
firmed in the present test). The agreement is good 

The A(T) effect will now be incorporated in the turbulence 
measurement. Since the completion of this investigation,! the 
author’s attention has been directed to the elegant analysis by 
Betchov® ® incorporating (1) the A(7) effect as well as (2) the 
end losses. He applied Eq. (1) locally along the wire and de 
rived corrections to the “perfect-wire’’ equations in terms of the 
“‘basic’’ data (wire dimensions and properties, A,, etc.). The 
present report obtains a correction to the measurement of turbu 
lence level in terms of quantities somewhat more readily meas 
ured in the laboratory by considering first-order differences from 
the mean operating conditions. The analysis includes (1) and (3 


the second-degree term in the resistance law 
R = R,{1 + a(T — T,) + BT — T,)? (3 
Item (2) is omitted, so the analysis evidently applies to wires 


with length-diameter ratios of 200 or greater (see reference 8). 


It is convenient to adopt the symbols? 
2 = VU/{{A(T)/B] + VU} 
dy (R — R,)/Ry 
a, =(R. — R,)/R, 


502 


tion of a quasi-equilibrium state of the flow field will be admitted 
Eq. (1) is expanded in a Taylor’s series keeping only first-order 
terms, in the manner of reference 3, giving 


AE/E = —}(ay’s)/[2(1 + ®)]} (AU/U (4 


where FE 1s the voltage drop across the hot wire rhe terms 
dy’, z, and ® correspond to the mean conditions / and [ Eq 
(4) is identical to Eq. (12) of reference 3, for a constant-current 


wire, when ® = 0. # incorporates (1) and (3 For typica 


operation ® can be approximated by 


j 


Eq. (5) is plotted in Fig. 2 for platinum and tungsten (a and 3 
from reference 10) In reference 1 the ® correction brought the 
two measurements within 5 to 10 per cent agreement 

The following remarks are in order 

(1) The ® correction is generally larger for platinum due t 

(2) For tungsten (3) terids to offset (1), but for platinum they 
are additive. 

(3) Roughly @ ~ a,’ for ay’ < 0.5. Since the signal is also 
proportional to ay’, it is not necessarily feasible to reduce a 
to reduce ®. 

The writer would like to acknowledge the helpful discussions 
with Drs. Stanley Corrsin and Leslie S$. G. Kovasznay 
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Optimum Stage for Axial Flow Turbines 


N. Van Le 


Aerodynamics Group, AiResearch Manufacturing Company, 
Los Angeles, Calif 
anuary 10,1955 


ime PROBLEM OF FINDING the stages of an axial flow turbine 
that would run at the highest efficiency under prescribed 


pressure ratio, diameter, speed and mass rate of flow has been 
approached on a trial basis, guided by experience 

A solution to this problem is presented, aiming at a systematic 
investigation, eliminating the cut and tried process and intro 
ducing possible use of computing equipment. Despite the three 
dimensional flow pattern in turbines, their efficiency is calculated 
with mean values. On this basis, a turbine stage characterized 
by a velocity diagram and its h-s representation 1s considered in 
Fig. 1, with the definitions shown 


The stage efficiency ns7 can be calculated as 


NST = W/W, =h h h h es 2¢ 1) 


The ratio 7 7 can generally be taken as equal to 1, or estimated 
from the degree of reaction of the turbine stage. The loss co 
efficients ¢., ¢r, derived from test data, generally, form an in 
tegral system—for example, that of reference 1. The loss co 
efficients depend mainly on the turning angle, solidity, clearance, 
aspect ratio, Reynolds Number, ete. If a type of velocity dia 


gram defined by its degree of reaction is selected, it is possible to 
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Fic. 1 
W = hy hys, Shaft work 
Wy, = ho, — hess — C.*/2¢, reversible work when the exit velocity 
C» is not assumed to be lost 
t; = Ah,/(C\?/2g), loss coefficient in stator 
tr = Ahpr/(W2?/2¢), loss coefficient in rotor 
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Fic. 2. Characteristics of impulse stage. 
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Fic. 38. Characteristics of impulse stage 


represent the characteristics of different combinations of this 
diagram. Fig. 2 shows a velocity diagram chart of the impulse 
type (zero reaction). For each point on this figure there is a 
corresponding velocity triangle which yields complete informa- 
tion for the calculation of ¢., fg. Application of Eq. (2) permits 
the establishment of a map of lines of constant nsr, Fig. 3 
Now, every point on Fig. 2 yields alsoa shaft work W’. From the 
given pressure ratio, the efficiency nsr can be calculated with 
Eq. (1) independently of ¢,, &e. Lines of constant nsr are then 
calculated and shown on Fig. 2 

Upon combining Figs. 2 and 3 and taking the intersects of 
lines of the same nsz, a turbine that satisfies the prescribed con- 
ditions is found. It operates at the considered value of nsr 
The intersects constitute an operating line along which the 
efficiency is known and, therefore, the highest one can be selected, 
Fig. 3. No other combinations of velocity diagrams exist 
that would satisfy the problem 

By changing the degree of reaction, a new operating line is 
established. The combination of all these operating lines will 
permit the establishment of the absolute optimum turbine stage, 
defined by its degree of reaction and the parameter C,/1 

It is readily seen that this method is adaptable to modern com 
puting machines as any judgment based on experience is not 
involved. It is conceivable that different ‘‘unit’’ turbines cor- 
responding to different sets of prescribed conditions can be in 


vestigated to establish relations between these units. 
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A Similarity Rule for the Interaction Between 
a Conical Field and a Plane Shock* 


Y. H. Kuo 

Graduate School of Aeronautical Engineering, Cornell University 
Ithaca, N.Y 

January 14, 1955 


N THREE-DIMENSIONAL SUPERSONIC steady flow past an ob- 
I stacle, it frequently happens that a conical field touches, and 
thereby interacts with, a plane shock, resulting ii a phenomenon 
known as a “triple point.”” In the limit of infinitesimal shock 
strength, problems of this kind have been treated by linear 


in which the "’ is disguised as a logarithmic 


theory, “triple point 
singularity superposed on an algebraic one at the Mach cone 


In reality, of course, the shock cannot change discontinuously but 


* This study was sponsored by the Mechanics Branch of the Office of 


Naval Research 
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decays to zero strength under the 
through a definite region whose 
strength of the shock. 
this phenomenon necessitates a re-examination of the approxi 
mation assumed in the linear theory 

In the case of a weak shock, the flow behind it ma 


V stall be 


considered as isentropic and irrotational. Let the velocity po 


tential be ¢ in the cylindrical coordinate system (p, 6, in which 


is in the free-stream direction If the flow is conical, ¢ = 
Uz{1 + f(r, 6 csand8 = VY V.2—1; land 
are, respectively, the free-stream velocity and Mach Number 





] where r = Bp 


The differential equation satisfied by the function f(r, 4), 
cording to Lighthill,? is 


fh— Pr — (7 +1)M orf — rfp + (1/2) (Ff —2 i 
(y —-D)M.2f — rf, (1/2) Uf - 2 

1/2)(r? + B2) (f.2 + 12h? 2M ao%*fd1 +f — rf-) — 
M 287,21 l—(y¥-—1)Ma%f—? 

(1-9, + ¢-“Heg) — 2M At Yo — Fee + C= 0 


Here C and y stand for, respectively, the cubic terms of the 
and the ratio of the specifi 
When the 
disturbances are small and no terms other than linear ones are 


exact equation, not containing f,,, 


heats. Subscripts denote partial differentiation. 


kept, the coefficient of f,, is approximated by 1 — r? and the 
equation becomes linear. For such cases, ry = 1 defines the 
singular Mach cone. It is no surprise that linear theory break 
down in the neighborhood of r = 1, inasmuch as the terms 


neglected are of the same order as 1 — r. To remove the ambi 


guities inherent in the linear solution in that neighborhood, non 
linear terms are indispensable 

If the small parameter of the flow is a, such as the angle of 
attack of a wing, the shock strength, and so the velocity dis 
turbances, will be of the order a. In the region of interaction let 
it be assumed that 
3 atk” F(E,n 7 
1—r=(ak)"— > 2 
@-—& = (ak)'s } 


where 4 is a constant giving the location of the interaction zone 
and k is another constant to be determined. It is now assumed, 
on the basis of the fact that the interaction takes place far from 
the solid boundary, that the character of interaction—namely, 


the exponents /, m, and n—depends solely on the local conditions 
and the differential equation. 
velocity disturbance ¢, — U = l(f — rf,) are both of order a, 
If, in the coefficient of f,,, 1 — r? is of the same order 


These two conditions uniquely 


Since the shock strength and the 


l—m=1. 
asf — rf,, then m must be unity 
define] tobe 2. The first term of Eq. (1) is therefore O(a). Now 
if fag is also of the first order, m must be 1/2. Therefore, by re- 


taining only the first-order terms, Eq. (1) simplifies to 
vith 


This shows that the effect of the interaction would confine itself 
to a region that extends in the radial direction a distance of 
O(a) and subtends an angle of 0(./a) and that two flow fields 
described by fi(r, 0) and f2(r, @) are similar if the product ka is the 
same for both cases. It is of interest to point out that, by a differ 
ent interpretation of the variables, this similarity rule holds 
equally true for unsteady two-dimensional flow, such as the 
diffraction of a plane shock by a corner. 

One of the boundary conditions that the solution to Eq. (3 
icross the shock 


must satisfy is the shock condition—namely, < 


(a) f is continuous, and, (b) 
Af, = (2K¥ + 1)] (vy -)DU—-Sf-) + 2Mao7f - 
2f, + 2M. ~1B%r — 1) — M,,~48°(dr/dd)?)] 4 ; 


where Af, signifies the jump in f,. From the similarity trams 


formation, f is O(a?) and therefore the shock condition reduces to 


o 


F: = 4& + 2(dt/dn)? 


action of the conical fielq 
dimensions depend on thy 


Consequently, a correct description of 
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if the disturbances in front of shock are zero 
If there is no disturbance outside the Mach cone and if both 
” ind 1 — r near the Mach cone are of the order a’, the 
came argument leads to the following similarity transformation 

1 /B\ka'F(é / 
=? ss Ra)*s +) 

v — dy = kan \ 


iffecting the validity of both Eq. (3) and Eq. (5 





The similarity considerations above show that in the (r, 3 
plane a three-dimensional purely supersonic conical flow in the 
neighborhood of both the “triple point’’ and the Mach cone r = 1 
is governed by not a linear but a nonlinear differential equation 
of the mixed type, with the transitional line Fe = 2 separating 
the elliptic and hyperbolic regions. These similarity rules lead 
us to conclude, without solving the differential Eq. (3), that the 

rturbed shock strength of a first-order plane shock can be ex 


pressed in the form 
Ap/p.~. = —\(y.,,2a 
with a shock curve of the form 
l—*y? Rak ” 7.1 
ind that the shock strength of the Mach cone without disturb 
inces in front of it has the form 
Ap/t =— WV. *ka?)/3) Fes, S 
witha shock curve 
l—r = (ka)*k(n 8.1 
Two special solutions of Eq. (3) have been studied They will 
be briefly discussed in the following 


A) Suppose: 


} H(u —_ Qult + of 1 — Lo n + mo)? T 2(n Q 
w= E+ & + o(n + m)* 
vhere g() is a polynomial of degree 4 and £, 70, and o are con 


st 


ants. The differential equation for //(«) is, from Eq. (3 


2u — H')H”" — (1 — Pod’ — 4e(1 — Qo) nu = 0 10 
where HH dH/dyu. This equation can be explicitly integrated 


with the result for 60 =~ 1 


H' —4opn)|H’ — (1 — 2o)p 1+ 20)/[2(1 2o)) = const. (11 
When bo 1, the solution can similarly be found and contains 
i logarithmic term 

If the constant o lies in —1/2 o < 1/2, the integral curve 
H', w) will have a modal point at «4 = O. For small yu, one 


branch is positive and the other is negative and intersects the 


singular line H’ = 2yu at uw < 0 with infinite slope. In the special 
case where there is no disturbance outside the Mach cone, o = 0, 
then F; = H'(u 


the flow across the Mach cone is continuous 


Thus, on the expansion side where H’ > 0, 
On the other hand, 
on the compression side where //’ < 0, the flow across the Mach 
cone must be discontinuous with a shock strength from Eq. (8 
proportional to —ka?3~'!H'(uo) and shock cone 1 — r = 
4~"(Ra)?H'(uo) according to Eq. (5 These results seem to con 
firm Lighthill's earlier conclusions regarding the character of the 
shock in the case of a triangular wing with angle of attack, lying 
inside the Mach cone 


B) By the similarity rule, assume 


F = (n+ no Gl? 12 


where 


With this form of solution, it is evident from the shock condition 
(Eq. (5)] that an admissible shock curve is ¢ = const. = ¢ In 
the case of a conical field interacting with a shock, such as the 
tip Mach cone from a rectangular wing at angle of attack, the 
velocity disturbance F: and the coordinates (£, 7) of the “‘triple 
point’ are, respectively, —1 and (—1/2, 1/\/2) from the shock 
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condition. With these values, it is easy to deduce G'(fo) = 2fo 
The shock condition |Eq. (5)] then determines f) = —(1/4), from 
which it follows that m = 1/(yV2 One can conclude finally 
that the perturbed shock curve is & = —(1/4)(n + 1/V/2)? on 


which the shock strength decays according to (y + 1/y 2 
The numerical integration of the resulting equation for G 


with the initial values & = —(1/4),G = (1/24), and G’ = —(1/2 
shows, however, that the integral curve (G’, ¢) forms a cusp 


after hitting the singular line G 2 + 4¢2 at & = 0. The 


solution in this neighborhood is of the form $~ C,, ¢)” There 
= 

is, at the moment, no evidence to support the existence of this 
singularity in the interaction zone, even though other properties 
of the solution are acceptable Further exploration along this 
line is certainly desirable 

It might be mentioned, in passing, that in trying to solve the 
same problem, Legras* chose a similar shock curve as the initial 
line to construct a solution in power series. Unfortunately, his 
process has not been pursued far enough to determine whether 
this form of solution is actually possible. In the light of the new 
results, it seems to be a matter of fundamental importance that 


Legras’ solution be reexamined 
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Effect of High Advancing-Tip Mach Number 
on Helicopter Blade Flapping 


Kenneth B. Amer and H. Hirsch 
Aircraft Division, Hughes Tool Company, Culver City, Calif 
January 14, 1955 


R™ RENCES 1 AND 2 discuss severe retreating-blade stall as 
the cause of helicopter rotor blades striking the fuselage 
or tail boom. The purpose of this note is to propose an addi 
tional, perhaps even more important, cause of blade striking, 
especially appropriate for modern helicopters 
Indications are that the blade-striking problem is becoming 
more severe with modern helicopters, and it is hypothesized that 
the blade-striking problem is accentuated by the current tendency 
toward operation at higher tip and forward Mach Numbers 
Recent NACA data indicate that airfoils such as the commonly 
used NACA 0012 experience a forward shift in center of pressure 
with increasing Mach Number before the well-known large rear 


ward shift occurs. It is believed that this forward shift in center 
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Fic. 2. Spanwise variation in angle of attack increment per 


unit root angle increment for a typical rotor blade. Advancing- 


tip Mach Number = 0.7. 


of pressure, combined with the higher dynamic pressures asso 
ciated with operation at higher advancing-tip Mach Numbers, 
results in a pitch-increasing twist of the blade on the advancing 
side and leads to a condition in which blade-striking occurs. 

In order to indicate the manner in which the forward shift in 
airfoil center of pressure with increasing Mach Number twists 
the blade and leads to the blade-striking problem, a preliminary 
quasi-static study on a sample rotor blade was conducted. In 
this study, the effect on the blade of an increased in blade pitch 
on the advancing side (90° azimuth), such as would be caused by 
upward collective and/or rearward cyclic control motion, was 
computed. Figs. | and 2 present the results of this study in the 
form of the spanwise variation in angle of attack increment per 
unit root angle increment. The outboard angle of attack incre- 
ment is greater than the root increment because of blade twist 
due to the lift increments near the tip acting forward of the center 
of gravity. 

Fig. 1, for an advancing tip Mach Number of 0.6, shows a 
small increment in angle of attack change due to blade twist 
Fig. 2, however, for an advancing tip Mach Number of 0.7, shows 
a large angle of attack increment due to blade twist. The large 
increase in blade twist indicated in Fig. 2, as compared to Fig. 1, 
occurs due to the forward shift in chordwise center of pressure 
associated with the increase in Mach Number and to the in- 
creased dynamic pressure. Also, it reflects the fact that at the 
higher advancing tip Mach Number, there is a significant de- 
stabilizing effect in that the blade twisting causes additional lift 
which, in turn, causes additional twisting. 

Figs. 1 and 2 indicate that upward collective and/or rearward 
cyclic motion, which increases the blade angle of attack on the 
advancing side, may cause a large angle of attack increase at that 
position per unit control displacement under conditions of high 
advancing-tip Mach Number. This large increase in angle of 
attack would cause the blade to rise considerably in the forward 
position. Because the blade is operating in a centrifugal field, 
this large amount of flapping represents a large amount of po- 
tential energy. As the blade continues to rotate the potential 
energy is converted to downward kinetic energy and again to 
potential energy represented by an unusually large negative 
blade-flapping angle in the rearward position, where striking may 
occur. (This situation is aggravated by retreating-blade stall- 
ing, as discussed in references 1 and 2.) Thus, magnitudes of 
control motion that are safe at lower tip and forward speeds may 
become unsafe when higher tip and forward speeds are reached 

Assuming that blade twisting due to forward shift in center of 
pressure at high advancing-tip Mach Numbers is an important 
cause of the striking problem, it would appear desirable to design 
the blades of high-speed helicopters with the chordwise center of 
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gravity moved forward about 3 per cent of chord from the usya 
position, at least over the outer 20 per cent of the span. This 
forward shift in blade chordwise center of gravity would elimina, 
the unstable cg-cp situation at high advancing-tip Mach nyy 
bers; thus, a fundamental cause of the blade-striking problen 
would be eliminated. At lower tip and forward Mach numbers 
the cg-cp relation would be stable, which is acceptable from sue} 
standpoints as blade dynamic loads and flutter. It is recognized 
that the proposed shift in chordwise center of gravity would 
effect control loads at lower tip and forward speeds becaus, of 
the separation of the center of gravity from the center of pressure 
However, it can be shown that this effect is minor in relation to 
other sources of control loads 

The above hypothesis arose out of work on high tip-speed 


rotors conducted under contract for the United States Air Fore: 
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Analysis of Thermal Stresses in Conical Shells 


Marvin Stern 

Design Specialist, Convair, San Diego, A Division of General 
Dynamics Corporation, San Diego, Calif 

January 18,1955 


SUMMARY 


The differential equations for the problem of induced thermal stresses in 
conical shells under axially symmetric temperature distribution are derived 
Contrary to past analysis of this problem, the development presented herein 
allows for variations in temperature profile, Young’s modulus, and skin 


thickness. The results appear as only two second-order linear differential 


equations that are relatively simple and allow for solution by machine meth 


ods for wide ranges of the above parameters 


INTRODUCTION 


I A RECENT PAPER! an idealized case of the subject problem has 
been analyzed. This past analysis will be extended here to 
include allowance for the following effects: (1) temperature 
profile to be arbitrary, (2) Young’s modulus to vary as a function 
of temperature, and (3) skin thickness to vary arbitrarily 

The development in this analysis parallels that of Tsien in his 


treatment of the problem for heated flat plates.* 
ANALYSIS 
(1) Displacements and Compatibility 
For displacements u in r direction and w in s direction (inward 
normal to surface) the strains at a distance 2 away from the 
median surface are: 


e- = (du/dr) — 2|(d?w)/(dr?)| | 1 


€p9 = [(u/r) — (w/tan y)] — 2(1/r) (dw/d? 4 


Let E(r, z) be the variable Young's modulus, v the constant 


Poisson's ratio, and a the coefficient of thermal expansion. Then 








Fic. 1. 
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(1) Static Equilibrium 


he usua the stress-strain relations are 
1 , Static equilibrium of an element requires that 
Thi € ) 1 E(r ' Or — voy al(r 1 | 
liminat ; (d/dr) (rN,) — Ny = 0 
h nun € 1 / [EG eee alr { 
: (d/dr) (rO + Ny/tan > = () 2.1 
roblen ld . \ 
nverting Vicids 
umbers : (d/dr) (rM,) — Msg — rQ, 0 
1 r 9\t 
m suc , — p2)} 4 a _ 
ne SUCe LAr, dis i. slags , he normal shear term Q, is eliminated in the following way 
Enized i |E(r, v)yal(r, | From the first two equations above we find 
Would . 2\t = (1.2 
) E(r, ay Ne si QO, tan y + A const 
Use of YT Rr {a ' T(r x 
' “Aly ~ i (ls 
ressure Since NV, = Q, = Oat the boundaries, it follows that 
tion to In order to work with sectional forces and moments, set \ —(, tan 
-speed \ Se d Therefore, our equilibrium equations become simply 
Force 
Ns So d d/dr rN Ny | on 
- 1.3 (d/dr) (rM, WV» — rN, cot { 
Vi = Se d 
B P -_ : , — 
V I . 1 (Il) Reduction to Two Equations in Two Unknowns 
rt ' ] TyH~- € 
The equilibrium Eqs. (2.2), together with Eqs. (1.6) which 
19 . : . . 
e, 908 define forces in terms of displacements, will now be reduced to 


Integrating Eqs. (1.2) and substituting Eqs and (1.3 


two second-order differential equations in terms of two dependent 
variables 


du (" <i ) Introduce & = dw/dr as one of the dependent variables and sub 
\ D T » — - 
hells La r  rtan 


stitute the definitions of Eqs. (1.6) into the second of Eqs. (2.2 


‘ai bits d/dr)Def [r(dé/dr)| + vt} + (d/dr)r My 
D, ( T )-a Dz} (1 /r)é t v(dt/dr)|{ 4 Wr+rN, cot 


7 Expanding and collecting terms, the first of the required equations 


it “ du 
N D ( ‘ _ = is obtained 
r r tan dr 
j j dg (! d | D )s (" l ¢, D 
1 du Ww T T og —_ —_ og , 
D ( : + as *) _ V; dr? , dr . dr r yr dr . 


sses in 4 dr lr 
rived ry 1.4 cot \ 1 di; _ 
herein - dD net . — os : D D. dr 7 
1 skin : : d? . r tan > 
ential In order to obtain the second equation needed, we first invert 
meth du 1 dw \ ~ 
D bon ) — iM the first two equations of Eqs. (1.4) to solve for u and du/d? 
dr? r dr 


Vo D (" 7 W ) ; cd - w cot 7 + } 1+ y)DIli}Nr (3.2 
1 has r r tan > d? 











a - du/dr il 1 — v2)Dolj (N, — vNy) -— j1 l v)Dol} N 
€ to 1 du dx ‘ ’ : ; 
ture D j re? Ip? — M The first of Eqs. (2.2) gives a substitution for Ny in terms of 
yr dr dr 7 ‘ . 
‘tion N, Inserting this substitution and eliminating u from Eqs 
3.2) by differentiation, we find 
: where 
1 his 
d , d d , 
| ] ° rN, — vN + £ cot > t \ 7 
D&W : E(r, d dr(1 — v*?)Do\d dr(1 + »)D 
l—» 
. ] ad , l 
| . V}—- > rN, \ 
DO : E(r, d (1 — v*)D dr 1+ v)D 
ard l i ; 
the . Expanding and collecting terms, the second of the required equa 
dD ; E(r, z)z? dz (15 tions is obtained 
, 
4 a d?N (° d ] dN l—vj)d Do 
l , + — log D — og 
\ Vil? EG T(r, d dr? r dr ; dr r dr 
| i 
ant ° (l—». (1 — pv) d 2 aN 7 | 
: l : ‘ cot 7 —& = log Do - \ — (3.3 
1€1 Vr? E(r, s)T(r, s)s dz Dor? > dr , dr 
I ah 
mal IV) Collection of Equations 
Bi choose our median surface for s = 0, dependent upon Our two second-order differential equations in the two de 
Ar, 2), sothat D, = 0 Eqs. (1.5) reduce to pendent variables &, N, are 
\ Dy \(du/dr tT vi(u/r) — (w/r tan y) } — Ny (d*t/dr?) + fir) (dt/dr + folr)éE + flr) N, = flr i 
4.1 
_ Vi Doy \(ae/r — (w/r tan 4 + v(du/dr) } — N\ 1 (d?.N,/dr? T Jal’ (dN,/dr t Jel? N, + IWnE = Js f 
1 1.6 
f D»\(d?w/dr*) + [v(1/r) (dw/dr)|{ — My where the f;(r) are known functions of r given in Eqs. (3.1) and 
Wy —Ds\(1/r) (dw /dr + [y(d?w/dr?)]} — My (3.3 
The boundary conditions are simply: atv = 1, fe, 
Eqs. (1.6) above give the representation of the net sectional forces . at 
ind moments in terms of the two displacement variables u, w 
ind the given thermal force and moment MV, = 0 = —Dz.\(dt/dr) + |(v/r)z}f — Ma (4.2 
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}) Simplification for Special Problem 
It will be shown here that a simplification of Eqs. (4.1) results 
from application to one of the more typical of the subject prob 
lems 
Consider the cone subjected to a temperature distribution re 
sulting from heat transfer through a laminar boundary layer 
Further, assume that the cone has been optimally designed in 
the sense that the skin thickness variation has been proportioned 
so that the resulting heat flux distribution permits the outer and 
inner skin surfaces to be isotherms. To within a first-order ap- 
proximation then, it is seen that the isotherm pattern is the same 
at all cross sections—i.e., the geometric picture of the isotherm 
pattern is merely scaled by the skin thickness h(r 
Using this scalar magnification property, we set 


The plate stiffness parameters become 


so that S dr = | 
Dir) = f E(r,s) dz = fE[T(r)|h(r) dr = 
Werf ia (5.1 
Dir) = f E(r, zs)? ds = f E[T(r)|h\(r)r* dr = 
hr) f° E(r)r? dr = exh 


The logarithmic derivatives become 





(d/dr) log Do = (d/dr) log h(r) | 
(d/dr) log Dz = 3(d/dr) log h(r f 7 


Further, since we have at our disposal the choice of temper 
ature reference, we choose this reference so that the thermal nor 


mal force vanishes everywhere—i.e 


o> 
} . 
ah(r : a 
= E(r)T(r i=0| 
l1— A 
We are left with only the thermal moment, which is now 
ah*(r) ( _. - - 
Mr(r) = 1 E(r)T(r)7r dr = e3h*(r o.4 
, oe 


Using Eqs. (5.1) through (5.4), the differential Eqs. (4.1) can 
be written in the following form: 


dt (. a ; i) (" =i d i) 
4 + ¢ ge H — — ocraiict- 
dr? r ar mai dr > r dr ' 


cot 4 2c, 1 d 
= og 
coh N, 2 hdr ae 
ro 
d2N, (- d ' : \ l—v)d , 1 
a — og nl - og ft. 
dr? r dr ” dr r dr . 
| v* cot 4 
& = () 
cohr? 


Eqs. (5.5) can be simplified still further by assuming that the 
optimally designed skin thickness variation can be approxi- 
mated by a curve of the form 


The resulting equations now become 


dt '. dé 1 
-+4+—(1 — 3n)— ——(1 + 3m) — 
dr? r dr r? 
cot y . _. 2nc 
rN. = , 
bic, bce 
.6 
d?N, 1 dN, ] 
—+—(3 +n) +-—(1 — v)nN, 4 
dr? r dr r? 
1 — v*) cot 4 
, t=0 
bco 
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he solutions to Eqs. (5.6) above can be easily approximat. 


by series expansion 


REFERENCES 


' Huth, J. H., Thermal Stresses in ( ical She Journal of the Aeron 
tical Sciences, Vol. 20, No. 9, pp. 613-616, September, 1953 

2 Tsien, H. S., Similarity Laws Siressing Heated Wing ] nal of t 
Aeronautical Sciences, Vol. 20, No. 1, pp. 1-11, January, 1953 


Transient Stresses in Airfoils Under Blast 
Loading* 


H. F. Ludloff and F. J. Marshall 


New York University and U.S. NOTS, China Lake, Ca 
January 24, 1955 


Bhar NOTE DEALS with the distribution of deflection, sh 
and moment in a simplified airfoil model (beam), if 





latter is struck by a strong, plane blast front. The airloads 
be considered as known theoretically or experimentally The 
type and intensity of the vibrations induced in the beam are 
be derived 

Among the various possible geometrical arrangements, two 


may be pointed out as of particular theoretical and pr 
interest (see Fig. 1 


(1) Strong blast advancing almost perpendicularly t 


surface [at angle (7/2) — e Case I: The shock front 
split into two parts. The part moving along the uppet 
will advance, as if against a shallow compression corner and ex- 


perience a Mach reflection;7 the shock moving along the lower 





surface will produce a weak expansion domain, as in 


his problem originated from a discussion while the senior a 












summer consultant at U.S. NOTS his preliminary investigation is t 
supported, in part, by the Mechanics Branch of ONR 
f See Ting, L., and Ludloff, H. I i imi iB Journal of the 
Aeronautical Sciences, Vol. 19, No. 5, pp. 317-328, May, 1952 
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CASE I 


Fic. 1. Two cases of blast loads. 
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irts of the shock front will end perpendicularly on the 


very smail distance AB apart 





face, on places which are 
other Since either shock diffraction produces only 





small deviations from the undisturbed pressures behind the 
shock front, it is seen that a strong concentrated force (fp; + € — 
f along the sur 


»),) from above, will advance with shock speed | 


II) Strong blast front advancing almost parallel to airfoil 
Case II: In this case, the lower part of the 


shock front will experience a regular reflection producing an ex 


irface (at angle « 


nely high pressure load, Pe, on a restricted part AP of the 
lower surface. The upper part of the shock front will be diffracted 


become quite weak near the upper surface, E Since 





P and E are increasingly far apart from each other, the airload 


g from below is spread over an appreciable area AP. Now, 





interesting feature is the fact that the end point of the load 
istribution, ?, advances over the lower surface with a speed 
sin €—i.e., at a speed which may be well in excess of the 
speed of propagation of the elastic waves produced in the air- 
beam in this case a discontinuous shear distribution will 


ivance into the undisturbed part of the beam, similarly as a 


In either case, the flexural vibrations of the airfoil (beam) will 

the predominant role. Since high frequencies are excited 
blast loads, the rotatory inertia and shear have to be taken 
In other words, the forced vibrations of ‘‘Timo- 
As known, the velocity 


0 account 
henko-beams’’ have to be computed 
f elastic waves in such beams tends with decreasing wave 
ength toward a finite limit 

hus, the analytic statement of the problem is to find solutions 


of the following differential equations 


ett-2-k, = —q; 5, 4 2Dr-ke — 7-2A-s = O; 
heres = V + 7-2\-wandk = M+7-2d-y l 
= shear, J = bending moment, w = angular velocity, y = 
transverse deflection. £& +7 = dimensionless axial and time co- 
ordinates. \ = L/2r, is the beam length in radii of gyration 


For our loading 
= P for 0 E< uct 
= Ufor E> u-t 
a : . . : : 
where the speed of the load (or load end point P respectively ), u, 
™m  « . . c : , 
may be larger or smaller than the limit speed of the elastic waves 
whine b4 
which is set equal to one 
Furthermore, there are two boundary conditions on either end 


of the beam, and four initial conditions 





5009 


FORUM 


The case of a concentrated moving load is easily deduced from 
the preceding case of distributed loading 

Solutions have been obtained in any one of the following three 
Wavs 

1) by the classical method of expansion 1n eigenfunctions 


Due to the discontinuity of the loading function, small wave 


lengths are excited with great intensity Hence, the slow con 
vergence of the series makes an evaluation of results prohibitive 

2) by the method of characteristics which is applicable, 
because differential equations (1) are hyperbolic rhis repre 
sentation is very appropriate for numerical evaluation at early 


times 


3) by application of the Laplace transform. In this repre 
sentation, it turns out that the solutions can be represented by 


a finite number of terms, so that an analytic discussion as well 


as a numerical evaluation could be successfully carried out 


rhe deflection curves y(£, 7) as well as the distribution of shear, 


j 


d rotation have been derived for u l 


bending moment ar 


This is comparable to supersonic, sonic, and subsonic conditions 


The analytic solutions obtained after (3) are of the form 


| rt Se — Ve] a - rT, : T d 2a 
ZAM « 
ind 

P ‘ 

Zr J 
1) ba 

T= Pr: We: — r,¢ r f/u »b 
> 
Here 
@ = rw =f r-J,{2Nt - — I 
AT. D(X V 7? — &2]; U(r, 6) = Vor 1 for tr 


Discussion of these expressions is greatly simplified, if char 


acteristic coordinates are employed 
2 and 3 show the early development in time of the dis 
9% = y El-P 106 and shear |’ = 


Figs 
tribution of deflection 
V-P-L 10 It is remarkable that for « > 1, the deflection 
is alwavs downward, while for x 1, there are two places (near 
the end of the elastic waves, and at one end of the wing where 
the deflection is upward. It can be shown from Eqs. (2a) and 
2b) that this is true for all times 


= ], the shear (and moment) distribution is discontinuous <¢ 


Furthermore note that at 
it 
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7 = & Comparing deflection curves at equal load length, the 


g 
deflection is larger for u < 1, as it should be. 
The absorption of energy by the beam as well as the energy 


spectrum, as functions of u and 7 are being investigated. 


On a Stress Tensor for Turbulence 


Willis H. Braun 
Lewis Flight Propulsion Laboratory, NACA, Cleveland, Ohio 
January 24, 1955 


r j SHE REYNOLDS EQUATIONS for steady turbulent shear flows of 
an incompressible fluid are 


+ (—puju;), ; (1 


pU;U;, 3 = BU, 34.35 — P, 
The l’; are components of the mean velocity; the u; are fluctua 
tion velocity components. For a solution of these equations to 
be obtained it is necessary first to find an independent relation 
between the mean flow quantities and the Reynolds stresses 7, 
puju There are several properties required of such a rela 
tion: 

(1) The formulation of the stresses must be a tensor. Hence, 
it should be composed of sums of products of tensor quantities 
such as the mean velocity components and their gradients and 
gradients of the mean pressure. No quotients or reciprocals of 
tensor components are allowed for they have no meaning in a 
tensor formulation. Reciprocals of invariants formed from 
the tensor quantities are permitted, but, because of the regular 
nature of the mean flow, it is assumed such functions are ex 
panded in power series of the tensor components. The net re 
sult is to bring the expression for the stress back to the stated 
form. 

(2) It must be consistent with the boundary conditions 
There are two boundary conditions on the mean velocity 
requiring that it vanish at solid boundaries and attain a maxi 
mum at some other surface. Therefore, the equation of motion 
should be of second order. This, in turn, requires that no de 
rivative of the mean velocity higher than the first appear in the 
stress tensor. 

(3) The magnitude of a stress component should remain in 
variant while its sign undergoes a change due to the reflection of 
a coordinate that corresponds to one of the directions associated 
with that component. If a coordinate is reflected which does 
not correspond to either direction associated with the stress, then 
the component should chaage neither in magnitude nor in sign 

In addition to satisfying (1), (2), and (3), it is desirable that 
part of the stress tensor should represent a dissipative process 
That part is expected to take the form of an irreversible process 
in a fluid—namely, the sum of products of diffusion coefficients 
and powers of the gradient of the velocity. In general, there 
should also be a portion of the stress which performs reversible 
work. 

The effects of requirements (1), (2), and (8) can be simply 
illustrated by examining the fully-developed turbulent flow be 
tween two infinite, parallel flat plates. Let x, be the Cartesian 
coordinate parallel to the plates and in the direction of the mean 
velocity U;, and let x2. be perpendicular to the plates. The com 
ponent of turbulent shear stress which acts to retard the flow is 
T», and according to (1) and (2) it should be of the form 


Ty = AmeUiUi, 2 + Aan iC Ui, 2)? +... + 

AgneU?U), 2 + Asn 12(Ui, 2)? + .. (2) 
Because the pressure gradients are either extremely small or con 
stant, their contribution to the stress is assumed to be negligible 


The 4A;; are tensor coefficients of proper dimensions, some 
of which are found to vanish by (3). The shear stress is then 
T» aed (Aone? + Asin 14 — ‘+ Ui, 2 -+- 


(Asuna U1? + Asien U1 eo ais ) (U4, 2) TT = (3) 








28 





ul +575 logy, +55 
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Eq. (3) resembles the general formulation of an irreversib 
process. The eddy diffusion coefficients, which are the terms j 
parentheses, are associated with the relative mean velocity by 
tween the fluid and the solid boundary. Their dependence on 
velocity difference rather than on a velocity gradient corresponds 
to the production of the stresses by the turbulent motions whic! 
are of such scale that the mean velocity changes appreciably over 
one eddy. The classical form of a dissipative process involves 
only terms linear in the velocity gradient. According to refer 
ence 1 the classical generalization of Eq. (3) for the dissipativ: 


stresses should be a deviator tensor 


(Tij)p = (ajjetrsU.Ui T Qijkl U._.U.U, U, 4+ [ : _ 
(1 3 )65;(@ppkirs U uy T @ppkly l UiU,U, T l ’ | 


The portion of the stress tensor which does reversible work 
has normal components only and is associated with the turbulent 
pressure 


(Tir = —(1/3)punttm 5; — #6 


Because the mean velocity gradient is associated with the dissi 
pative portion of the stress, the turbulent pressure will be repre 


sented as a function of the mean velocity components alone 
we = B;U;U; + BignUiUjU.Ui 4 6 


The total turbulent stress is the combination of Eqs. (4) and (4), 


Tu = (Tuo ~ vhs (7 


The coefficients a;, and 8;;... must in general depend upor 
the dimensionless parameters of the flow, on whether the bound 
aries are solid or fluid and on the shape of the boundaries if 
solid. These are the factors that determine the history of the 
flow and, hence, the turbulent structure at a point. At present 
it appears that the coefficients and their parameters’ dependence 
can only be ascertained from experimental data. Consequently 
the stress formulation |Eqs. (4), (6), and (7)] seems to have no 
immediate applicability in the sense of predicting mean velocity 
profiles in flows not already experimentally investigated. Never 
theless, it may be of value in suggesting correlations of exper! 
mental data. 

In this regard the application of Eq. (4) to fully-developed pipe 
flow yields a new expression for the universal velocity profile neat 
the wall. The primary dimensionless parameter of the flow 1s 
the Reynolds Number based on shear stress velocity and radius 
Re = pU*Ly - 


tained from Eq. (4) by performing a dimensional analysis of the 


The shear stress (using the mixed tensor) is ob 


a;j, ... and introducing dimensionless variables. The resulting 
form is 
T,' = p(l*)? > f,( Re) (,*)?"1, 2* 5 
n 1 
where 1,* U,/U* and the f,(Re) are unspecified functions 


The integrated equation of motion is then 


l 


1 on + 1° 


l - 
u,* + u,*Re f,(Re) (u,*)2* = 5 Re - 7 3 


n= 


Comparison with experiment? indicates that the leading terms 
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in the sum in Eq. (9) are integral powers of (1*)* Re Then 
ear the Eq. (9) reduces to 

uy* + a(m*)? = y (10 
where 3 Re(1 — r) and r is the dimensionless radial coordinate 
of the pip Eq. (10) for the mean velocity profile near the wall 


is compared with the logarithmic law (with constants as given in 
3) in Fig. 1. It varies smoothly from the linear law, 
laminar region to a nearly logarithmic curve where 


reference 


y, in the 


the turbulent stress predominates. 
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The Obukhoff Spectrum of Homogeneous 
Isotropic Turbulence 


Knox Millsaps 
Chief, System Dynamics Research Branch, Wright Air Development 
Center, ARDC, Wright-Patterson AFB, Ohio 


January 31, 1955 


A M. OBUKHOFF! HAS SUGGESTED that the inertial transfer 
¢ term in the spectral equation for homogeneous isotropic 
turbulence may be given in terms of the spectral distribution 
by postulating that the microcomponents of the turbulence 
constitute a Reynolds stress which produces a mean rate of shear 
in the macrocomponents of the turbulence. It is naturally as- 
sumed that the Reynolds stress is proportional to the energy 
in the microcomponents and that the mean rate of shear in the 
macrocomponents is effectively measured by the root mean 
square of the rate of strain of the macrocomponents. The 
Obukhoff postulate when applied to the universal equilibrium 


situation leads to the following spectral equation: 


c= f° F(R) dk + V2 f, F(k) dk lf * REF (B) ak | 


wherein the customary notation? has been used. This integral 
equation is usually solved by numerical methods, and a closed 
solution does not appear to have been obtained. Consequently, 
it may be interesting to note that the substitution used by J 
Bass* in his treatment of the Heisenberg spectrum may be em- 
ploved and reduces the integral equation to a cubic algebraic 


equation. The resulting closed solution is 


l 
Fk) = — [A+B ria oe 
3k? 


for 


ind 
F(k) = 0 for (ey?/4v3)'/* < k < 

where 

bs 

- ‘ (27<*y4 + 18ey?v3kt + Qy5k ) 

1324 
B= eks 27e*y? + 4tevkt 

Sy 3 


and the primes denote differentiation with respect to k 
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On the Normal Force Produced by Multifinned 
Tail Assemblies in Supersonic Flow 


Shapiro 
Redstone 


J). Leith Potter, William D. Murphree, and Norman M 

Research Division, Ordnance Missile Laboratories, 
Arsenal, Huntsville, Ala 

February 7, 1955 


A’ THOUGH REASONABLY SATISFACTORY METHODS have been 
devised for the calculation of normal forces on tail assem 
blies having two equally loaded fins 180° apart, very little prog- 
ress has been made toward solution of the corresponding multi 
fin problem where four or more fins may be acting. Two re 
ports!: ? dealing with the subject have come to the attention of 
the authors. The earliest of these is by Luther and Jaeger 
and contains experimental results for several families of four 

six-, and eight-fin tails. The second report, by Miles, is a theo 
retical analysis using slender body theory. The latter report in 
cludes generalized results and apparently represents the only 
study of its type. However, the limitations of slender body 
theory are such that the writers, while appreciating the value of 
this work, wished to study the problem further. For this reason, 
it was decided that a brief but systematic supersonic tunnel 
program would be conducted for the purpose of obtaining addi 
tional experimental data on the lift of multifinned tail assemblies 
and that this would be followed by an attempt to generalize the 
data in a form suitable for design purposes. This note presents 
the data obtained thus far and a tentative proposal for an inter 


ference parameter which successfully correlates all the experi 
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Fic. 2. Comparison between predictions of slender body theory 


and some experimental data for clipped-delta plan forms 

















A-B IS THE TRACE IN THE BASE PLANE OF A MACH 
CONE HAVING ITS APEX AT A ROOT CHORD LEADING EDGE 


Fic. 2. 


mental data in reference 1 as well as that secured by the authors’ 
experiments 

The models tested in the Redstone 7- by 7-in. tunnel are de- 
scribed on Fig. 1. A trapezoidal or clipped-delta plan form was 
chosen in order to provide for subsonic and supersonic leading 
edges with adequate fin area for easily measured forces. The 
body was made as long as possible considering test section size, 
and calculations by the method of characteristics indicated that 
the nose should have negligible influence on the flow at the fin 
locations when angles of attack are small. In the test program 
already completed, all models were mounted with two fins in the 
i.e., the four-, six-, and eight-fin configurations had 
Measurements 


pitch plane 
two, four, and six effective fins, respectively. 
of normal! force and pitching moment were made at angles of 
attack equal to 0, 1, +2, +4, and +6 degrees with numerous 
repetitions to verify accuracy. Reynolds Number based on free- 
stream conditions was close to 1.94 X 10° per in. during all tests 
and Mach Number was 2.86. <A description of the NOTS ex 
periments with fins of rectangular plan form is given in reference 
1 and will not be repeated here. 

The performance of finned bodies will be discussed in terms of 
the quantity Cyg which is the slope of the normal force coeffi- 
cient curve at zero angle of attack. The writers determined 
values of Cyg by plotting (Cy/a@) as a function of the angle of 
attack, a, and fairing a curve through a = ().. After this, the 
finless or bare body value, Cygg, was subtracted from Cyg to 
give Cy7rg, which represents the total contribution of the fins 
including all interference effects. Having these data from ex 
periments, the next step is to compare results with Miles’ theory.” 


This is done on Fig. 2, where 
€, = Cynrq (n fins) /Cnrq (4 fins) 
and 
a, = tanw/tan pw 


the angles being defined on Fig. 3. 
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Referring to Fig. 2, it is evident that the theory may be yse 
for predicting approximate levels of the ratio, ¢,, but there ar 
features of the data not recognized by slender body theory For 
example, the detrimental effects of fin-on-fin interference may } 
reduced when the overlap of pressure fields is lessened. (Co 
siderations such as this led the senior author to formulate thy 
interference parameter described below 

The intention in this case was merely to search for any signif 
cant trend in the data when plotted as a function of som 


parameter which reflects the degree of mutual interference } 
tween fins. Asa simple first step, the ratio of a particular Mact 
cone radius, y, to the length between fin roots, /, was selected 
(See the notation shown on Fig. 3.) This ratio is 

I = (y/)) = c,/[BD sin (a/n 
where 


B= V M-—-1]1 


n > 2, total number of fins 


Presumably the body-on-fins and fins-on-body interference com 
ponents differ depending on the number of fins, but it will be seen 
later that the parameter, /, at least partially compensates thes 
effects 
In addition, Cy 7q of the six- and eight-finned tails was ad juste 

for the reduced angle of attack and projected area of the fins not 
positioned in the horizontal plane. If a, represents the geometric 
angle of attack of a rolled fin with zero incidence, it may be show: 


that 


a, = sin”! (sin a cos @ 
where 
a@ = angle of attack of nonrolled fins 
@ = angle of roll from horizontal 


Effective projected area of a rolled fin was obtained by multi 
plying panel area by the cosine of ¢@. These adjustments are 
meant to remove the penalties simply due to multifinned tails 
not having all fins operating at maximum angle of attack and not 
having all individual fin normal forces acting in the total normal 
force direction. The remaining differences between tails similar 
except for number of fins should be chargeable to interference 
effects 


fins on each tail are ineffective, 


) 


To summarize, in the present case where a — () and 2 















Me | 4 fins) = a 
5S, (4 fins) = 2S, 
a, (6 fins) = 0.866e 
S. (6 fins) = 3.46.8 
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be used a, (8 fins) = 0.805a lO 
é 
















































































lere are S. (8 fins) = 4.838 
\ For . 
saneil All values of Cy7e used hereafter are based on a, and .S,, the ’ 
lay D afiea : - ‘ . : 
Ces latter symbol representing total effective fin area. Sy is the area u 
ite the of cne fin O8 U 
The clipped-delta data and typical examples of the NOTS ‘ . 9 
signif rectangular-fin data are shown on Fig. 4. Considering the 
aa normal scatter in experimental data of this type, the parameter, 
som sis : ice ; 
nee | is very effective. On the strength of evidence shown on Fig ’ 
r Mad {it appears that relatively accurate values of Cy7q for six- and ral 
act , , P : : 
eight-finned tails can be obtained merely by plotting available O06 . ~~ 
lecter igh . I 6 © oO oe 
four-fin data as a function of 7. Pending further analysis, a ° 
stronger conclusion will not be suggested, but the apparent lack Pl 
f better methods for dealing with the problem seems to warrant u'Vv’ fe | 
this preliminary report ——_ PA 
O4LU ‘ 
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showt 
Note on Velocity Fluctuations in Equilibrium O 2 4 6 
= . ; ; ° 
Turbulent Shear Flow Along Solid 
Boundaries* 
Fic. 2. Comparison of w’/l, (experimental points) with 
Heinrich B. Helmbold V u'v'/U, (curve) as deduced from faired experimental curves 
chool of Engineering, University of Wichita. Wichita, Kan of u’/U, and v'/U> for plane boundary layer flow with no pres 
c © 140cc sure gradient (plane wall, see reference 2 
February 8, 1955 
nulti 
'S ar | , 3 
hess ECENTLY TWO CAREFULLY PERFORMED experimental inves an . , 
tails . 7 s (a) The root-mean-square values of velocity-fluctuation com 
| tigations have been published by J. Laufer! and P. §$ é : i 
d not = * : ihe eye ponents are interconnected by the relationship (Figs. 1 and 2 
Klebanoff? which reveal common features of fully-developed 
ormal . : : 
ae | turbulent flow in pipes and plane boundary layers with zero we = Vi l 
milar ‘ : 
pressure gradient It may be expected that these features are ee 
rence within the ranges of 
19 haracteristic not only of these special types of flow but of 
ind 2 . ‘ ° “4 : > » «2 — ' 7 . 
equilibrium turbulent shear flow along solid boundaries in gen 0.1 n S 1.0 in the pipe (independent of Reynolds Number 
eral and 
*This study was performed under Office of Naval Research Contract O.1<yn 0.5 in the plane boundary layer (intermittency factor 
Nonr 201/01 y « l for n>O0O 33) 

Here » denotes the nondimensional distances from the wall, re 
zz T ferred to the pipe radius KX and to that distance 6) from the flat 
| | | : . ; ‘ 

— 7 t | plate where the shear-correlation coefficient vanishes (6 
} | 1.26 in Klebanoff’s notation), respectively. Whether the range 
csi atecall : mre ; 
of validity for the plane boundary layer actually is restricted by 
| 
= ae intermittency in the outer region of the layer, or if it merely 
seems so because of a falsifying effect of the intermittency on the 
measurements, remains open to discussion 
I (b) The ratios of the root-mean-square velocity-fluctuation 
os components in a plane parallel to the mean velocity and per 
pendicular to the walls of the pipe and flat plate are a ‘‘uni 
versal’ function of the nondimensional distance from the walls 
a (Fig. 3) in the form 
r 4 u’ =k + 0.379 (2 
ee oe |} | The constants and ranges of validity are dependent on Reyn 
| | } | ie | : 
: F: 6 olds Number; 
ptt | | | | | | bb UR s0.25xi0%_ 
k = 0.438 (pipe), 0.05 < < 0.5, U,sR/v = 0.25 X 108 
| | 6 > (pipe), ) n , 1 “ 
= |___|e R/v =0.025xI0 nbs a . 
i | oa Up R/ —— k = 0.45 (pipe), 0.15 < » < 0.6, UpR/v = 0.025 X 108 
——}_ } | J | | | ve | oi 2 a k = 0.44 + 0.01 (plane wall), 0.08 < » < 0.6, U) 6o/v = 
| ; 0.1 X 108 
—t + + 4 i + ee | _+—_— —— ae a 
1) = | | ;} | U4 The significant point about Eq. (2) is not the special function 
wil + oh ail 1 = = 7 : 
7 5 1O 1.5 20 chosen for the representation of the experimental values, but 
I Fic. 1. Comparison of w’/U; (experimental points) with the existence of a “‘universal”’ function 0/1’ = f(n For greater 
the ata ‘ : ; Sa aiisaia > Og ‘ ti aiutiinn oe eee 
Vu'v'/U; (curve) as deduced from faired experimental curves distances, 7 0.5, the function deviates from Eq. 2 and tend 
ou /U;, and v'/U; for pipe flow (see reference 1). to f(1) = 1 in the center of the pipe. The apparent difference 
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Fic. 3. The “universal” function v’/u’ = f(n 
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Fic. 4. Lateral distributioa of 
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Fic. 5. Lateral distribution of turbulence scales in pipe (L/R) 
and plane wall (L/6;) flows. 
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between the functions for pipe and plane boundary layers in thj. 


region is subject to the same considerations about intermittency 
effects mentioned above 

(c) To complete the description of the field of velocity fly 
tuations in that region where these two cases are strictly cop 
parable—i.e., where intermittency does not interfere (0.1 < » 


0.33)—it also may be stated that the longitudinal fluctuations 


u’'/U, are approximately the same with the pipe flow and th 
plane-wall boundary layer at sufficiently high Reynolds Numbe 
(Fig. 4). Together with the statements (a) and (b) this means 
that in an important region of both flows the kinematic proper 
ties of the fluctuation fields are essentially the same; the only 
major difference is the magnitude of the shear correlation eo 
efficients, about 0.42 with the pipe flow and 0.50 with the bound. 
ary layer at high Reynolds Numbers. In other words, the fix 
tuation fields agree essentially with respect to the intensities 
However, they are much different with respect to the scales as 
can be seen by comparing the nondimensional quantities L/R = 
[uv/u'(dU/dn)] for the pipe and L/5) = [uv/u'(dU/dn)| for the 
plane wall (see Fig. 5). The length Z can be interpreted as a 
measure in the y-dimension of the large eddies in turbulent mo 
tion. 


REFERENCES 
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NACA TN 2954, June, 1953 
? Klebanoff, P. S., Characteristics of Turbulence in a Boundary Layer with 
Zero Pressure Gradient, NACA TN 3178, July, 1954 


Limitations of Linear Theory in Predicting the 
Pressure Distribution on Triangular Wings 


Elliott D. Katzen 
Aeronautical Research Scientist, NACA, Ames Aeronautical 
Laboratory, Moffett Field, Calif. 


February 11,1955 


7 PURPOSE OF THIS NOTE is to call attention to the fact 
that even for relatively thin, low aspect ratio wings the 
local surface slopes may be sufficiently large so that linear theory 
cannot be used for calculating the lift distribution, even at angles 
of attack as small as 3°. Experimental measurements that 
illustrate this were made on a triangular wing 5 per cent thick 
with an aspect ratio of 4, at a Mach Number of 3.36. The wing 
was tested in the 1- by 3-ft. supersonic wind tunnel at the NACA 
Ames Laboratory. The semispan wing was a reflection-plane 
model mounted on a boundary-layer bypass plate. The wing 
section, plan form, and positions of the pressure orifices are 
The wing was composed of similar sections of 
Each section 


shown in Fig. 1. 
5 per cent thickness ratio in the stream direction 
was a circular arc to the maximum thickness point at mid-chord 
and a parabolic are from the mid-chord point to the blunt trailing 
edge. The trailing-edge thickness was 50 per cent of the maxi 
mum thickness. Asa result of this wing geometry, the thickness 
distribution was not conical with respect to the apex in that the 
ratio of local thickness to chord was not constant along rays from 
Therefore, the pressures due to wing thickness at 
zero angle of attack were not constant along rays from the apex 

If linear theory is applicable, the pressure distribution due to 
angle of attack (lifting pressure) is independent of that due to 
On this basis, the wing 


the apex. 


wing thickness at zero angle of attack 
pressures can be considered as the sum of the pressures due to wing 
thickness at zero angle of attack (thickness pressures) and those 
on a flat-plate wing of the same plan form at angle of attack 
Since the flat-plate pressures due to angle of attack should be 
constant along rays from the apex, it follows that if linear theory 
is applicable for the complete wing, the pressure due to angle ol 
attack should be constant along rays from the apex regardless of 
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whether or not the thickness pressures are conical with respect to 
this point 

The lifting pressures are shown as a function of the direction of 
rays from the apex in Fig. 2 It can be seen that there is poor 
wreement between linear theory and experiment; the lifting 
pressures are not at all constant along these rays 

A more complete understanding of the nonlinear effects on the 
outboard portion can be obtained by reconsidering the geometry 
of the wing and the pressures. The wing is a cone with respect 
to the tip; the geometry can be described by directions from this 
point since there is no characteristic length in the problem. 
Therefore, the wing pressures should be constant along rays from 
the tip for that part of the wing forward of the region of influence 
of the apex. If the thickness and lifting pressure distributions 
ire plotted as functions of the direction of rays from the tip as 
in Fig. 3, it can be seen that the thickness and lifting pressures 
ire nearly constant along such rays. Comparison of the results 
of linear theory and experiment shows good agreement for the 
thickness case. For the lifting case, however, linear theory 
indicates a constant value of lift regardless of ray direction, 
whereas measured values vary systematically with change in ray 
direction. Apparently the local surface slopes are sufficiently 
high for this Mach Number to give rise to a nonlinear component 
of lift in much the same way that nonlinearities are produced on 
unswept two-dimensional wings of the same section at the same 
angle of attack. The linear theory predicts constant lifting 
pressures for the unswept two-dimensional wing In contrast, 
the shock-expansion theory accounts for the finite thickness of 
the section and predicts increased lift at the leading edge and de 
creased lift at the trailing edge, resulting in a shift in the center 
of pressure toward the leading edge. It is evident that, as for 
the unswept, two-dimensional wing, nonlinear solutions that 
properly account for thickness effects are necessary for accu- 
rately predicting the lifting pressures on the triangular wing at a 
Mach Number of 3.36 

For that part of the wing forward of the region of influence of 
the apex, a nonlinear solution can be obtained by extension of 
the shock-expansion method. This has been done by Vincenti, 
of the NACA Ames Laboratory, who considers the curved sur- 
face of the wing as the limit of an infinite number of planar sur 
faces intersecting along rays passing through the wing tip. The 
changes in flow across each element, consisting of an intersection 
ind subsequent flat surface, are then found from considerations 
of geometry and simple-sweep theory. The result, in the limit 
of the curved surface, is a pair of simultaneous ordinary differen- 
tial equations whose solution is readily obtained by numerical 
means. Theoretical curves obtained by this method for the 
Present wing at 0° and 3° angle of attack are shown in Fig. 3 
It can 


ve seen that the variation of lifting pressure with ray 
direction, previously noted, is estimated more accurately by the 
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shock-expansion theory. The differences in absolute value of 
lifting pressures as estimated and measured may be attributed to 
small deviations in the model ordinates from the design values; 
theory indicates qualitatively that this is the case but does not 
provide a quantitative correction. Despite these differences, it 
is apparent that even for relatively thin low aspect ratio wings 
the surface slope may be sufficiently large that linear theory 
cannot be used for calculating lift distribution and that this is 
the case for the present 5 per cent thick wing at Mach Number 
of 3.6 and angle of attack of 3 

Deviations from linear theory also occur inside the region of 
influence of the apex. For this region, further development of 
nonlinear theory is necessary in order to calculate the pressure 
distribution 

These results have shown an example of the errors involved in 
superposing thickness and lifting solutions for a triangular wing 
in the manner assumed to be valid in linear theory. Wings of 
other plan forms have similar regions for which nonlinear solu 


tions are also needed 
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Fic. 2. Lifting pressures on a triangular wing with an aspect 
ratio of 4 at 3 degrees angle of attack 
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Fic. 3. Thickness and lifting pressures for that part of the wing 


forward of the region of influence of the apex 





A New Method of Solution of Equation of Compressible Flow 


Ezic Betti 
Laboratorio di Aerodinamica, Scuola Applicazione Aeronautica, Florence, Italy 


February 28, 1955 


ise 
V \ E CONSIDER THE TWO-DIMENSIONAL equation of compressible flow in which the parameter m 1s introduced: 


Oy Fa Oe ¢ Og V, 
l1— m 


») 


¢ is the potential function 

Vz = O¢/Ox Vi, = O¢/oy a* 
dy is the sound velocity at rest. For m = 0 the equation is that of incompressible flow and for m = 1 the equation is that of compressible 
flow. We try to determine the compressible flow, subsonic or supersonic, around a given airfoil by utilizing the solution of incompressible 
flow around the same airfoil and proceeding in the following manner: 

We subdivide the interval 0, 1 of m parameter in » partial intervals of size 4; the corresponding values of m are indicated with m, 
Mo,... Mi, ... my; (where mi+1 — m; = uw, m, = 1). The generic function, solution of Eq. (1) with m m; is indicated gy; and the 
difference gy; — gi;-1 = wx 

We utilize the equations obtained by neglecting the terms in yp? and higher in the following .* 


O*(gi-1 + wX j 1 mM, O(gi-1 + wNX;) |? , 07( gi—1 + wX;) y, m Ol gi-1 
ax? | u*(m;) Ox f Oy? | a*(m;) oh 


m; O7(0;-1 + wX;) O(¢gi-1 + wX,) OC gi-1 + wXi) Ri O*gi-1 J i is 
~ a*(m;) Oxey Ox 2 Ox? a*(mji-1) Ox? 
m1 gi-1?( . 2mi-1 OP Gi-1 OGi-1 OG; 
aX my-1) dy? f " aX mi-1) OxXOY Ox Oy 
ee ee Sarre 


The simplified equations are indicated with the number (3) between the brackets and the above equations with (2). The difference method 
will be used to determine X, of the first of linear partial differential equations (3), go teing the solution of incompressible flow Then from 


the knowledge of gy; = yg + wX; using the second of (3) will be determined XY» and so continuing X;, X,, . : The solution of the 
Now when pu ~ 0, that one will be as the solution of equations of com- 


last equation (3) willresult: gy, = go + wXi + wpXo +... wX,,. 
pressible flow because the difference between (2) and (3) resides in the terms of order yw? and higher neglected in the equations (3). 


The importance of the terms will be made extremely small decreasing the size of intervals—namely, for Au > 0. 


* a2(mj) ao? — [Cy 1)/2] [(Ogs/Ox)? + (Ogi /Oy)?] 


On Use of the General Bending Formula where 


A =SfdA, 0, =fydA, Q, =SxdA, L. =SydA, I, = 


W. P. Rodden 
Missiles Division, The RAND Corporation, Santa Monica, Calif SxdA, Try = SxydA 


February 14, 1955 The simultaneous solution of Eqs. (5), (6) and (7) yields the 
following expressions for the coefficients: 
HE GENERAL UNSYMMETRICAL bending formula for plane : i 
, . p B = (1/A) [M.(0.0, [.,A) + M, (I,A O,”) 4 
stress, o, given by the simultaneous solution of the three cies v p i .)) @ 
‘ P pera ‘ ‘ sas ‘ ‘ (O,T, @) )} (8) 
equations of section equilibrium without the simplification of li . 4 
centroidal coordinates offers a reasonable savings in computing , @ (1/4) [M, (1,4 + M, (Q:Q, 
time and memory space required for stress calculations on a high- 
speed digital computer ; 
Ss ed digi 1 pute = (1/4) | 6A eS WU, 
Writing the equilibrium equations about arbitrary reference 


axes, x and y, we obtain the applied axial load, 
where 


P= dA (1 
Se = AI I,I,) + Q. (Qzl, Ot.) 4 


where dA is an element of section area, the applied bending mo- Oy (Q,1 Q,T zy) 


ment about the x axis, Although these expressions for the coefficients are somewhat 
| Sf oydA 9 more lengthy than the corresponding ones in the usual unsym- 

. metrical bending formula, the computing effort saved by avoid- 

and the applied bending moment about the y axis, ing the calculation of the centroidal coordinates of typical com- 


Me ms ff oxdA . plex aircraft wing and fuselage structures is extremely worth 


while 
Substituting the plane stress relationship, The centroidal moments of inertia needed in deflection cal- 
: culations are obtained from the transfer formulas: 
o0 = Bx + Cy+D 
¥ = (12) 


the equilibrium equations become 


BQ, + CO; + DA =P 
Bl., + Cle + DO, = —M. 
Bl, + Cl, + DQ, M, é where ¢ = Q,/A and ¥ 
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